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1620-1311 Linear Programming System 
(1620-C0-04X) Program Reference Manual 

This is a general-purpose programming system providing 
sophisticated mathematical techniques to determine the most 
efficient use of various resources in carrying out alternative 
activities. A number of programs are stored on the 1311. 
Each routine is called into core storage by procedural con- 
trol cards when it is to be executed. The sequence of con- 
trol cards defines the solution procedure. 

Besides direct optimization, the program provides extensive 
post-optimal analysis to indicate the effect of changes in 
constraints, costs, or technology. Automatic procedures 
have been incorporated for determining the frequency of in- 
version and appropriate data accuracy, but both may be ad- 
justed by the user to accomplish specialized purposes. 
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PROGRAMMING SYSTEM ABSTRACT 



The IBM 1620-1311 Linear Programming System is a general-purpose programming system 
designed to provide the 1620-1311 user with sophisticated mathematical techniques to deter- 
mine the most efficient use of various resources in carrying out alternative activities. The 
system is composed of a number of programs which are stored on the 1311. Each program 
routine is called into core storage by procedure control cards (or agendum cards) when that 
program routine is to be executed. The sequence of the control cards defines the solution 
procedure. 

Besides direct optimization, the program provides extensive post-optimal analysis to indicate 
the effect of changes in constraints, costs, or technology. Automatic procedures have been 
incorporated for determining frequency of inversion and appropriate data accuracy, but both 
may be adjusted by the user to accomplish specialized purposes. 

GENERAL DESCRIPTION OF PROGRAMMING SYSTEM 



Linear programming (LP) is a mathematical technique for determining the optimum use of 
various resources (such as capital, raw materials, manpower, and plant or other facilities) 
to attain a particular objective (such as minimum cost or maximum profit), when there are 
alternate business activities to be carried out. LP can be used to allocate, assign, schedule, 
select or evaluate the uses of limited resources for various jobs: blending, mixing, dis- 
tributing, controlling, budgeting, ordering, bidding, cutting, trimming, pricing, purchasing, 
planning, and transportation of goods from plants to warehouses to outlets. 

The IBM 1620-1311 Linear Programming (LP) System optimizes an objective function subject 
to the constraints of a system of independent linear equations. For the minimum 1620 con- 
figuration, the system will accept problems having up to 100 equations (including alternate 
objective functions) and up to 1, 000 columns (including structural variables and multiple 
right-hand sides). The revised simplex algorithm with special features is used to obtain a 
two-phase solution. Source language for the system is IBM 1620 SPS II-D. 

The IBM 1620 Monitor I System controls the calling in of the LP System, and its I/O opera- 
tions. Upon completion of an LP run, control is returned to Monitor. This means that a 
linear programming problem can be stacked with other jobs run under the 1620 Monitor. 

The IBM 1620-1311 LP System is composed of a number of programs which are stored on the 
disk. Each program is called in only when it is to be executed. The cards used to control 
this procedure are called agendum cards. These cards are prepared by the user and entered 
into the system through the card reader. As the agendum cards are read sequentially by the 
computer, the appropriate agendum routines are called in from the disk and executed to per- 
form the input, solution, output, and utility functions required for an application. 

The input data originates on cards and may be stored and maintained on disk for subsequent 
reprocessing. Data revision is accomplished easily and does not require that all of the data 
files be reloaded onto the disk. Problems can be run using selected subsets of data, with the 
required rows and columns being referred to by name. The data files can include multiple 
objective functions (also called cost rows) and right-hand sides; these are also selected by 
name when a problem is run. 

Processing may be interrupted during calculation, and later continued from the point of 
interruption. 
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FEATURES 



The specification of input data is simple and provides the following flexibility: 

• Variables, right-hand sides, objective functions, and constraints may be selected from 
card or disk input by item name. 

• Data is maintained on disk in compact form. 

• Revision of an entry on the disk does not require that the entire file be re-entered on the 
disk. 

• Several problems may reference different subsets of the same card or disk data files. 

• Slacks are automatically generated and named. 

• A starting basis may be given; if it is not, a slack basis is used. 

• The use of SHARE standard matrix card format provides substantial data compatibility 
with other linear programming systems. 

The following advanced mathematical methods are used to increase processing speed and 
accuracy: 

• Revised simplex (product form of inverse) with a choice of automatic or user control of 
inversion frequency. 

• User may maximize or minimize any objective function. 

• Assignment and variation of floating-point length and tolerances during solution process- 
ing may be under automatic or user control. 

• Minimum floating-point length is five digits; maximum is 18 digits. 

• The bounded variable algorithm is applied to range and <) constraints (and bounded 
variables). In alloy scheduling (blending), for example, the silicon composition of an 
alloy (product) may be constrained to the range of 5% to 7. 5% by one range constraint. 

• Input data files are maintained in compressed floating-point form. This decreases cal- 
culation to about 30% of normal floating-point time. 

Restart is incorporated into the system as follows: 

• Processing may be interrupted; a subsequent run can continue from the point of 
interruption. 

The following extensive checking features are incorporated: 

• Input is checked for duplicate row identifications, split or duplicate columns and right- 
hand sides, and duplicate coefficient entries. 

• Each time a reinversion is to be done during solution, the problem is checked for ac- 
curacy; if error tolerance is exceeded, the mantissa length and tolerances are adjusted 
automatically. 
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Any output may be selective or include all pertinent slacks and variables. All output quanti- 
ties are in fixed-point form, with a specified number of decimals. Output is optional, and 
the output unit may be selected as follows: 

• The user may specify output device for the iteration log as well as iteration frequency. 
During solution a sense switch can be used to inhibit output of the log for certain 
iterations. 

• The output device for other reports, which may also be specified, is independent of that 
used for the iteration log. The output options include: solution amounts for variables 
and slacks, or both solution amounts and optimum price range for variables and slacks; 
opportunity (shadow) prices for nonbasis variables and slacks; sensitivity (restriction 
coefficient shadow price) values for nonbasis variables; and error check for solution 
using rounded output values. 

Post-optimal analysis capabilities include: 

• Marginal value — the change in the objective function per unit increase in the corre- 
sponding right-hand-side constant, assuming that a change of basis is not required to 
maintain feasibility. 

• Reduced cost — the amount by which the cost of a non-basic structural variable would 
have to be reduced before this variable would tie the optimum basis. A reduction in ex- 
cess of this amount would push the variable into the basis. 

• Cost range — the objective function interval of a basic solution structural variable in 
which the basis is optimum; a change in the objective function of the variable which is 
outside that range would force a change in the basis to preserve optimality. 

• Submatrix summary — the partial sum of the coefficients of a matrix, evaluated at the 
current solution; the solution contribution of a subset of the variables on a subset of the 
rows. 

After optimum solution, reoptimization can be efficiently accomplished with the following 
features: 

• Using a new right-hand side (or sides) 

• Changing values in the right-hand side 

• Changing values in the objective function 

• Changing values of coefficients in the matrix 
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MATHEMATICAL METHODS SUMMARY 



The 1620-1311 Linear Programming System uses the bounded variable, product form of the 
inverse, simplex method. The simplex method is based upon the fact that if there are m 
equations (or rows) in the constraint matrix and they are linearly independent, then there is 
a set of m columns (variables or vectors) which are also linearly independent. Hence, any 
right-hand side can be expressed in terms of these m columns (called a basis). The simplex 
method works with these basic solutions, stepping from one to another (by adding one column 
and removing one column from the basis on each step or iteration) , until a solution (called a 
basic feasible solution) is obtained that meets all of the criteria, including the requirement 
that all the column values be non-negative. The bounded variable feature allows the user to 
specify limits on the activity levels for any or all of the variables. Either or both upper and 
lower bounds may be specified. 

After a basic feasible solution is found, the simplex method steps along, examining a series 
of basic feasible solutions, to find one that satisfies the requirement that the value of the 
functional (or objective) row be a maximum or minimum; this is called the optimal solution. 
Not all linear programming problems have an optimal solution. If there is no solution at all 
in non-negative variables, or none which keeps the variables within their specified bounds, 
that LP problem is said to be infeasible . If a feasible solution is found, but the constraint 
rows do not confine the value of the functional row to finite values, the LP problem is said to 
be unbounded . 

In the process of trying to find a feasible solution, row transformations are made on the con- 
straint rows of the matrix. Then an optimal solution is made on the objective function row. 
The transformed coefficients in the objective function row (also called cost row, or functional 
row) are the relative cost factors, or dj's. It has been found advantageous to modify the 
method to obtain solutions faster. The simplex method uses the dj's to determine which 
vector to bring into the basis. The product form of the inverse is used to update this vector 
in terms of the basis. The product form of the inverse is a means of representing the in- 
verse as a product of matrices of a special form. Actually, only one column of each such 
matrix is required to represent the matrix. A new matrix (column), called an eta, is pro- 
duced each iteration. It is used on all succeeding iterations until the matrix is inverted 
again. When reinversion is done, the etas are started again; this reduces the number of etas 
(and also the density of the etas) that must be processed during each iteration. 

In reinversion, the column selected for pivoting is that which "minimizes the number of 
coefficients modified at each step (excluding those which are eliminated at the particular 
step)". * This is essentially a revised tableau simplex iteration, with the pivot choice de- 
termined by the number of ensuing operations. The tableau begins with the current basis 
variables in their original input data form (not updated). The tableau is updated by simplex 
iterations on the remaining variables, and the selected variable is pivoted and entered as an 
eta. This method is particularly suited to a data processing machine with high input/output 
speeds relative to compute speeds . 

MACHINE CONFIGURATION AND PROBLEM SIZE 



The 1620-1311 Linear Programming System will operate on a 1620 Model 1 or Model 2. 



* Bibliography, reference 11. 
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1620 Model 1 



Minimum configuration: 

Core storage of 20,000 digits 
Indirect Addressing (Feature 4650) 
Automatic Divide (Feature 1285) 
Additional Instructions (strip, fill, 

move flag) (Feature 1021) 
1622 Card Read Punch 
1311 Disk Storage Drive 

Optional: 

Automatic Floating-Point Operations (Feature 1288) 
1443 Printer 

Core Storage of 40, 000 or 60, 000 digits 
Additional 1311 Disk Storage Drives 

1620 Model 2 

Minimum configuration: 

Core storage of 20, 000 digits 
1622 Card Read Punch 
1311 Disk Storage Drive 

Optional: 

Automatic Floating-Point Operations (Feature 1289) 
1443 Printer 

Core storage of 40, 000 or 60, 000 digits 
Additional 1311 Disk Storage Drives 

Machine Size Max. No. of Rows Max. No. of Columns 

20K 100 1000 

40K 250 2000 

60K 400 3000 



SYSTEM CONFIGURATION 



The source language used for the 1620-1311 Linear Programming System is IBM 1620 SPS 
II-D. The system is called in under control of the IBM 1620 Monitor I System. (See Bibli- 
ography, reference 12.) The I/O operations in the LP System are under Monitor control. At 
the completion of an LP run, the LP System returns control to Monitor. Therefore, an LP 
problem can be stacked with other jobs to be run under Monitor. 
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GENERAL SYSTEM CHART 




DATA INPUT, AGENDA AND REPORTS 



AGENDA 



User control of the 1620-1311 LP system is exercised through special procedure control 
cards which are known as agendum cards (agenda is the plural) . Each agendum card causes 
a particular program (or set of programs, where program size is a factor) to be read from 
disk storage into main storage. That program is then executed; it will in turn call for addi- 
tional data input if required (either from cards or disk storage) or even call other agendum 
programs. Certain programs will produce selected output on cards for LP solution restart 
or printed reports for user analysis. 

The various agenda accomplish four principal functions; they are explaned in detail in "1620- 
1311 LP Agenda" in the order listed below: 



Data Preparation 



INPUT. Calls the input program to read in data from cards, set up 

the problem on disk, or to reference existing information on 
disks. Problems may contain alternate objective functions 
and multiple right-hand sides. 

REVISE Calls the revise program to change identification information 

or data values already stored on disk. Changes can be made 
to the row and column identification files, the matrix and 
right-hand- side element files, and the stored basis. 



Optimization 



MAX. . 



Calls the solution program (Dual and Primal with use of 
Invert) to obtain the optimal solution to produce the maximum 
value for the selected objective function. 



MIN. . . 



Calls the solution program (Dual and Primal with use of 
Invert) to obtain the optimal solution to produce the minimum 
value for the selected objective function. 



INVERT 



Calls a program which obtains an inverse of the current 
basis; this is seldom done by the user since it is handled 
automatically by 1620-1311 LP. 



Report Preparation 



OUTPUT 



DO. D/J 



Calls the output program to report the optimal solutions on 
cards, typewriter, or printer, and gives basic structural 
and logical variables activity level. 

Calls the programs needed to determine the relative values 
of structural variables in terms of entering the basis, and 
indicates the relative value of changing the constraint limits. 



COST. R 



Calls the programs needed to determine the cost range for 
each structural variable which would not cause a change in 
the solution basis. 
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CHECK. 



Calls the check program to determine solution accuracy, or 
to obtain selected submatrix sums. 



Control and Data Maintenance 



ASSIGN 

END JOB 
ERASEp 



Establishes certain control values for mantissa length, in- 
version frequency, error tolerances, use of output devices, 
and iteration printout frequency. 

Transfers control back to the 1620 Monitor. 

Eliminates the index value which permits the use of a stored 
basis. 



GETOFF 



MAP . . . 



Punches a restart deck when processing needs to be inter- 
rupted prior to obtaining an optimal solution. 

Calls a program which lists the size of the problem (number 
of rows and columns) , and the location of the input data on 
disk. 



PAUSE. 



SAVE.p 



Interrupts proper processing at the end of a particular 
agendum to permit user review and insertion of changes. 

Stores the optimal solution as a basis for later use; the 
current inverse can also be saved. 



DATA INPUT 



There is a variety of data inputs required or useful in solving LP problems. The data input 
cards are the means by which the problem data are placed in disk storage or referenced 
from previous problems already in disk storage. The data inputs always follow an INPUT, 
or a REVISE agendum card. 

There are input indicator cards, input identification cards, and input value cards. They are 
discussed in detail under "Data Preparation" in the order listed below: 



ROW. ID indicator 
Row identification 



Introduces the row identification file. 

Contains the row name and row type (objective function, 
equality, less -than, greater-than, or range constraint) for 
each row in the matrix. Omitted rows are indicated. 



COL. ID indicator 



Introduces the column identification file. If not used, the 
1620-1311 LP System will generate column identification 
from the matrix element cards. 



Column identification 



MATRIX indicator 



Contains the column names and column usage (active or 
omitted) . 

Introduces the matrix element data value file. 
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Matrix elements 



Contains the actual data values and structural variable 
bounds (upper, lower, or both) for each active or potentially 
active matrix element. Both the column name and row name 
are indicated for each value. 



FIRST. B indicator 



Introduces the first or only right-hand-side element con- 
straint value file. 



NEXT. B indicator 



Introduces each additional right-hand-side element constraint 
value file. 



Right -hand-side element 



BASIS, indicator 



VARBLS indicator 



Contains the actual limit values for each constraint row, 
e.g., the constraint to which the row is equal, or the upper 
limit, lower limit or range constraint for an inequality. 

Introduces an advanced initial basis. Usually recycling of 
OUTPUT or GETOFF generated cards. 

Introduces the basic structural variable file. 



Basic structural 
variables identification 



Contains the names of the structural variable (column name) 
to be included in the initial basis together with a basis entry 
type to indicate whether it is to be entered at its upper bound 
or at an intermediate value. 



SLACKS indicator 

Nonbasic logical 
variable identification 



ENDATA or EOF indicator 
REPORTS 



Introduces the nonbasic logical variables files. 



Contains the names of the rows corresponding to the logical 
variables which are to be omitted from the initial basis. A 
basis entry type indicates whether the variable left at its 
lower intermediate or upper bound. 



Concludes the data input files. 



The user can receive various reports from the 1620-1311 LP System that provide useful in- 
formation concerning the solution, the sensitivity of the solution to changes in the values, and 
accuracy of the results. These reports are produced by the OUTPUT, CHECK. , COST. R and 
DO. D/J agenda. Under control of ASSIGN, the report information may be on the typewriter, 
on cards, or on a printer. 

The following list summarizes the reports which can be issued; details are covered in "1620- 
1311 LP Agenda" under "Report Preparation". 



OUTPUT 



Contains as the first card (or a line of print on the typewriter 
or printer) a BASIS, indicator card. This deck can be used 
as an advanced basis in solving subsequent related problems 
when placed in the input deck since the format is compatible. 
For the basic structural variables, the activity levels are 
given. The activity levels are given for basic logical varia- 
bles, and the simplex multipliers for the nonbasic logical 
variables. 
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Gives the row errors, the maximum error (the row error of 
the row having the largest error), and includes the specified 
row limit or limits and the solution row activity level. 

Contains the cost (or profit) interval over which the objective 
function coefficient may range and the structural variable 
still remains in the optimal basis. In addition to the upper 
and lower limit of the range, the report gives the variable 
that would enter if the upper bound were exceeded (by a cost 
change) , and the change in basis that would occur if the cost 
were to go below the lower bound. 

Contains the reduced costs (Dj's for nonbasic structural 
variables and for basic logical variables which are forced 
into the basis because of upper or lower constraint limits. 
This report also gives the current cost and the basis value 
(difference between current cost and reduced cost). Each 
of the rows in this problem is given the type of row (+, -, 
=, RANGE), the row name, and the marginal value. 

Punches certain output cards which can be used as an initial 
basis. 
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USING THE 1620-1311 LINEAR PROGRAMMING SYSTEM 



The basic ideas of the IBM 1620-1311 Linear Programming System are easy to learn. This 
section presents the Linear Programming System as if it were being taught in the classroom. 
Sample problems are used to illustrate the various features of the system; after each prob- 
lem there is a discussion of the routines used therein. 

This section does not discuss all of the features available with the Linear Programming Sys- 
tem. The complete information with all of the options appears later under "Data Prepara- 
tion" and "1620-1311 LP Agenda". You may choose to refer to the appropriate parts of those 
sections as you study this section, or you may give this section a reading for the broad out- 
line before becoming concerned with all the details. The latter course is recommended for 
new linear programming users. The very experienced linear programming user should skim 
through this section for terminology and general format; he should progress quickly to the 
succeeding sections. 

The Linear Programming System is a high-level, problem-oriented system. It is a very 
specialized system in that it will solve only linear programming problems. The use of the 
system involves the same three basic steps as the solution of any problem in data processing. 
These are: input (concerned with placing the data that is to be processed into memory), 
compute or process (concerned with processing the data in predetermined ways to arrive at 
the answers), and output (concerned with recording the answers). 

Control of this system is accomplished by the use of control cards, called agendum cards. 
These cards initiate the execution of routines of the Linear Programming System. 

Five examples in this section explain the program features. These examples are based on 
the sample problem recorded later in the manual. They progressively introduce the various 
agenda with their options, going from simple formats and uses to the more complex. Each 
sample problem is legitimate and complete in itself. Given the initial conditions stated for 
the problem, the problem solution can be keypunched to produce results using the 1620-1311 
Linear Programming System. All the information presented here assumes that the user has 
carefully read the IBM manual "An Introduction to Linear Programming" (E20-8171). This 
manual contains a glossary which should be used to define any unfamiliar terms. 

Before tackling the first sample problem let us review the construction principles underlying 
the 1620-1311 Linear Programming System: 

1. The system has been designed for ease of use and is simple to prepare and modify. 

2. The system is disk-oriented so that moderate-size problems can be handled most ef- 
ficiently; data and solution logic can be saved between runs; reference can be made to 
previously solved problems. 

3. It is efficient to perform multiple runs with slightly modified input data at a fraction of 
the cost of the original run. Special post-optimal analyses can be performed conven- 
iently to determine the extent to which data input can be changed without requiring major 
solution changes. 

4. All operation is under the 1620 Monitor system. 

5. The actual optimization is controlled automatically to produce accurate results efficiently. 
The sophisticated user may override the automatic controls on mantissa length, toler- 
ances, and inversion frequency through use of ASSIGN and INVERT. (These two agenda 
are discussed under " 1620-1311 LP Agenda" and are not covered here). 
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6. A large number of input and computational errors are identified and appropriate mes- 
sages displayed, but the system solution continues with appropriate assumptions unless 
certain major errors occur. 

SAMPLE PROBLEM 1 

The following problem is concerned with blending certain raw materials to produce a desired 
end product. This first sample problem shows how the information can be formulated to 
describe the situation and placed in disk storage for later solution. 

An aluminum alloy smelter wishes to produce a particular alloy at minimum 
cost. The alloy must, however, meet certain chemical constraints. The 
smelter has available various scrap materials and some industrially pure 
materials. There are five scrap materials which are stored in separate 
bins; each scrap has been chemically analyzed. The inventory of each bin 
is known, as is the cost of each scrap material. Two of the bins contain 
powdered scrap; since these have the desirable property of oxidizing rap- 
idly in the furnace, it is required that at least a certain amount from each 
of these bins be used in the solution. 

The information stated above gives a framework in which to associate the actual problem 
data. 

1. The use of each of the seven raw materials is bounded by the available inventory (an 
upper bound) and by the minimum required usage (lower bound). 

• Bin 3 has 800 pounds in inventory and must have 400 pounds used for effective oxida- 
tion. Therefore bin 3 has an upper and lower bound. 

• The aluminum ingot (a pure material) has an unlimited supply available, and no 
minimum requirement. Therefore the pure aluminum usage is unbounded. 

A series of mathematical statements is required to explicitly indicate any applicable 
bounds. It is always implied that the minimum usage is zero unless otherwise stated 
(e. g. , we can't use a negative amount of aluminum). 

400 ^ BIN3 ^ 800 

^ ALUM (need not be stated explicitly since 1620-1311 LP will do this automatically) 

BIN3 and ALUM are the names of the structural variables whose value is the amount of 
usage that particular raw material will get. 

2. The cost of each raw material will be multiplied by the usage of that raw material to 
determine the total cost of the planned mixture. 

• The material in BIN3 costs $. 17 per pound. 

• Aluminum ingot costs $.21 per pound, etc. 
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These values are combined into an objective function statement. 



. 03 (BIN1) + . 08 (BIN2) + . 17 (BIN3) 

+.12(BIN4)+.15(BIN5) + 

. 21 ALUM + . 38 SILCON = VALUE 
VALUE represents the value of the objective function 

3. A material balance equation states that what comes in must go out. Therefore, the sum 
of the weight of the raw materials used must equal the end product weight. 2000 pounds 
of the alloy are required. 

BIN1 + BIN2 + . . . + ALUM + SILCON=2000 

4. There are six elements whose presence in the alloy is to be controlled. Each raw 
material has been analyzed as to its composition (by weight) in terms of these elements. 
The amount of a particular element in the alloy will depend on the usage of each raw 
material and its content of that element. 

• BIN2 has . 04 Iron (FE) 

. 05 Copper (CU) 
. 04 Manganese (MN) 
. 03 Magnesium (MG) 
. 75 Aluminum (AL) 
. 06 Silicon (SI) 

• The total must be less than or equal to 1. 00, since the composition is on a per- 
pound basis. In this case it is less than 1. 00, indicating that other elements are 
present. The amount of such other elements is not of interest in the solution of this 
problem. 

• There must be no more than 3% Iron (FE) in the end product. Therefore the alloy 
cannot contain more than 60 pounds of this element. 

FE ^ 60 

• It is required that all composition values be zero or positive so unless otherwise 
stated there is an implied lower constraint of zero. So iron has an upper con- 
straint limit of 60 and an implied lower constraint limit of zero. 

• There must be at least 75% of the element aluminum in the alloy. 

AL ^ 1500 

• The amount of aluminum actually present is the sum of the amount of the element in 
each raw material multiplied by its aluminum content. 

. 70 BIN1 + . 75BIN2 +...+. 97ALUM ^ 1500 

• A similar constraint statement would be written for each element. Some elements 
might have a range of permissible values, e. g. , both an upper constraint limit and 
a lower constraint limit. 
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To summarize, these are the following explicit statements required to describe the alloy 
smelting problem: 

• Bounds on the structural variables (the raw materials) 

• An objective function (cost of each structural variable) 

• Constraint and balance equations or inequalities (these are collectively called 
constraints) 

Let us see how the user places the alloy blending problem in disk. The first step requires 
the user to write a series of equations and inequalities to express material balance and con- 
straints; he also must write one or more objective functions; he must record the upper and 
lower bounds for each structural variable. This information then should be represented in 
matrix form for clarity and compactness. From the matrix the actual data input cards are 
prepared. 



The general form of the matrix is shown below: 

columns 



objective 
function 



rows 



constraints 



bounds 



structural variables 



Matrix Elements 
(coefficients) 



right- 
hand 
side 



constraint 
limits 



/right- 
jhand- 
jside 
(elements 



To see how this actually works, consider the following matrix, which represents one 
formulation of the alloy blending problem. 
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Bounds 



ROW NAME 


BIN1 


BIN2 


COLUMN NAME 
BIN3 BIN4 BIN5 


ALUM 


SILCON 


RHS NAME 
ALOY1 


VALUE 


. 03 


.08 


. 17 


. 12 


. 15 


. 21 


. 38 




YIELD 


1 


1 


1 


1 


1 


1 


1 


=2000 Limits 


FE (Iron) 


. 15 


.04 


.02 


.04 


.02 


.01 


.03 


^ 60 


CU (Copper) 


„ 03 


.05 


.08 


.02 


.06 


.01 




5 100 


MN (Manganese) 


.02 


.04 


.01 


.02 


.02 






^ 40 


MG (Magnesium) 


.02 


.03 






.01 






^ 30 


AL (Aluminum) 


.70 


.75 


.80 


.75 


.80 


.97 




2 1500 


SI (SI (Silicon) 


.02 


.06 


.08 


.12 


.02 


.01 


.97 


250 to 300 


Upper (Inventory) 
Lower (min req'd) 


200 


2500 


800 
400 


700 
100 


1500 









The matrix shows the coefficients of the variables from the equations and inequalities. For 
instance, look at the row called CU (Copper). These coefficients represent the following 
inequality: 

. 03 BIN1 + . 05 BIN2 + . 08 BIN3 
+. 02 BIN4 + . 06 BIN5 
+. 01 ALUM =£ 100 

The blank under the column SILCON for this row means that there is no measurable amount 
of copper in the pure silicon ingot. 



So the second step in preparing the data input for an LP solution is to prepare a matrix of 
coefficients from the various equations and inequalities. This is the time for names of rows, 
columns, and right-hand sides to be standardized. 



The third step is to prepare the data input forms and the fourth step is to keypunch and 
verify the information from these forms. 



The data input forms on the next pages fully represent the information in the preceding 
matrix. Notice these important conventions: 



1. All information must appear in defined columns. 



2. All names and data must be written exactly as indicated, with correct spelling and 
spacing. 

3. Letters and numbers must be carefully handwritten to avoid confusion (e. g. , the number 
and the letter O, the number 1 and the letter I, the number 2 and the letter Z). In typ- 
ing, be careful to avoid confusion between the number 1 and the letter 1. To avoid mis- 
takes, periods are often used as spacers instead of blanks (e. g. , MIN. . . ) 
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IBM 



1620/ 131 1 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 



JOB NO.. 



JOB TITLE 



.ANALYST. 



DATE. 



PAGE 



0F^_ 



KEYPUNCH- PUNCH DEC. PT. IN COL. 24 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






COMMENTS 


'.iiii 




I3i 


191,,, , 2 . 4 


31 136 |4I 
1 i , i 1 i i i i 1 i 


.46 |5I 


156 |6 1 |66 .73 80 
' « 1 1 1 — 


O IK \m 1 t r\ 
KiPiWh-iXjH 


1 . , , I 


1 i i i i 








1 1 I I .1 


— ' — 1 — 1 1 1 




ViA L 11 F 

VlniLtu iC i 








llll 


_L_J 


o 


YiI.C.L.D. 






# 


111,1 


1 1 1 1 1 


I +. 


nc i i i 






m 


1 1 1 




1 1 1 1 1+ 


CiU, , , , 




( 


t 




' ' ' ' ' ' 


1 . 1 1 


MiN. , , , 

i!L!5j i i i 










llll _J 


1 . . , .*■ 


MiG, , , , 








iii!! 


1 1 1 1 1 




AiL. , , , 

"M*" 1 — 1 1 — 1 — 




1 


1 


i i j i i 


1 1 1 I— 1 


1 1 1 1 l + 


Sil, , , , 




1 ... i'i .... . 


, 


i i i i i 






1,1,, 


( 


, 





i i i i 


C.0,L,.X.D 


1 


III,, 





1 


1 


i i i i i 


i i i — i — i — 


1 


&tl ,K1, 1 , , 







1 ,, , , 


i i i i i 


i i—i — i — i — 




BlT,N,2, , 


, 


, 


, 


i i i i i 






B,I.N,3 , 











I . i ■ 


i i i i i 




Baxf, , 














B.I.NJ, , 














/\|L,U,M, , 














SiI.L.C.M 








i i i i u 


i i i — i — i — 




i , , , , 














I , , , i 














l i i i i 














i , , . , 














i , , , , 








i i i i i_ 


i i i — i — i — 




1 i i i i 














l i i i i 










i. i i — i — i — 




t , , , , 










i i i — i — i — 




l . i i i 








■ i i i i 






1 , i i i 








t i t i i— 



IBM 



1620/ 131 ( LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 



JOB NO. 



JOB TITLE 



.ANALYST. 



DATE. 



PAGE JLoFik 



KEY PUNCH: PUNCH DEC. PT. IN COL. 24 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






COMMENTS 


1 . i i i i 


*l . , . , 


'3| , . . . 


191,,, , 2 . 4 


31 1 I 36 , . , . I 41 , 


I 46 ) 51 , 


|56 .I 61 ,166 _ ,73 _ 80 


MAT.R.I.X 


1 . . . , 


1 , , . , 











1 1 1 1 1 




B|X.N.» , , 


VlA.L.V.E. 


i , . , ...0,3, , , , 








1 1 1 1 1 




BiI.N.l. , 


Yll.E.L.D. 


i , , ,1,- 


1 ,2,0,0,. , , , , , , 


1 .,,,«,,,,, 


1 1 1 1 1 




B.I.N.l, , 


FiC. . , , 


1 , , , , • ,i ,s, , , , 








1 1 1 1 1 






Ciu, , , , 


1 , , , ,.,0,3, , , , 








1 1 1 1 1 




Bil.w.l. , 


M|N, , , , 


1 , , , ,-,0,2, , , , 








1 1 1 1 1 




ftlJU. . 


M|G, , , . 


1 , , , . - ,0.2, , , , 





1 , , , ,. 


1 1 1 1 1 




BiI.N.t. , 


AiL, , , , 


l , , , ,-,7,0, , , , 








1 1 1 1 1 




B|1,N,I, , 


S|I, , , , 


1 , , , ,.,0,2, , , , 


, 


1 


1 1 1 1 1 





i , , . , 


1 i i i , 





. | , , , ,. 





1 1 1 1 1 




ftl.N.2, , 


YlA,L,U,£ 1 


1 0,*, . , . 


1 , 


1 , 


1 1 1 1 1 


i i , . . 


B,l,M,i, , 


f.l.E.L.P, 


1 ■ , .1 


|2,S,0,0, ., 


1 , , , ,., 


1 1 1 1 1 




3\1M,Z, , 


F| E . , , , 


, , , 


1 


1 . , . 


1 1 1 1 1 




M.VJ.2, . 


:.u, , , , 


i . . . ...0,5. , , , 










B,L,M,2. , 


K|H, , , , 


l , , , ,.,0,4, , , , 












Mi6. . . , 


l , , , ...0.3, . , , 










BiI,N,2, , 


A,L, , . , 


I , , , ,.,7,5, , , , 










iiI.N.2. , 


Six, . . . 


i 0,£, , . , 










I , . , . 


i ... ■ 












1 , , , , 


i , , , , 












1 , i i , 


l i . . i 












1 , , , , 


i . . . . 












1 , , , , 


i i . i i 












1 , , , , 


I . . i i 












1 , . , , 


I , , , i 












1 , , , , 


l.i,. 












1 , , i . 


i , i . i 


l , i i , • , i i i i > 










1 , , , . 


1 , i , ■ 


1 , i ■ i_* i i i i i i . 






1 1 1 1 L_ 



1620/131 r LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 

ANALYST DATE PAGeJLoF.^. 



KEY PUNCH- PUNCH DEC. PT. IN COL. 24 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






COMMENTS 


''' ' ' ' 


7 I , . , , 


l3 1 i I I I 


1 i , , i*i , , i i , 


11 11 1 11 11 1 ' 


. . I 46 . . . . J 51 , . , 


|56 |6 1 |66 |73 80 






\/, A 1 It P 


1 i , , i • ,' i ' i i i i 






1 L 1 1 1 




RjT kt a 


v.t c l n 

T 1 1 ,fc, i-i*', 


lii 1 1 1 • . i i i i i 


<( ft /*, 


I i",*^ i^i • , i i i i i 


L J. | 1 1 L_ 






r\E, , , , 


1 i i i i * i^i^i i i i 


1 ' 1 1 '*' 1 1 ' ' 1 




1 L 1 1 L 


11111 




r.i i 

H«i i i i 


1,,, , , , 


1 1 1 1 




L .1 1 1 I 


11111 


AT M a. 


M.LI 


1 , , , i • i^i * i i i i 


1 ' 1 ' '* 


'''''' 


1 1 1 1 L. 




AT Ki a 
Dtl,ri,Oi i 


AiL 


lii 11*1 i i i 


'''''' -* — 


' 


1 1 1 i u 


I I I 1 I 


RiT v\ a 


, , , 


I,.,,* ,o,?, 


1 1,1 


1 ,,,,,,,,, , 


1 1 1 




1 1 1 1 1 


1 i i i 




' ' ' 


' 








1/iA I I1C 
VlftiLiViti 


i , . ,1 z 










HiT M U. 




1 , , ,/ , ., 


1 ,7,0,0, 


I , / ,0,0, 






fl* i»W| # i | 


Fi£, 

rid , , i 


i , , , ,.,0,4 , , 
i i i i i • i", 1 , i i i 










BiI.U.M-, , 


CiU, , , , 


..o,z. , , , 


> 








RiT KI U 


til Id 


i .0,2 

III''* iV^lt-i I i i 










BiI.N.H. . 


AiL, , , , 


i , , , ,.,7,5, , , , 










BiI.N,4.. , 


S.i, . , , 


i /,2, , , , 










i , , . . 


i ., ,. , , 












i , , . > 


i , . , , 












i . . . , 


i , , , , 












i , . . . 


i . . , . 












i , ■ , , 


i , , , , 












i , , , . 


l i , . i 












i , . . , 


i , , , , 












1 i i . . 


i , , , , 












I,,,, 


l i , i i 












I , , , , 


I . . ■ ■ 












i , , , . 


l i i i i 












1 . > i . 


1,1,1 


l i i i i * i i i i i i 














1 i _i i i • i i i i i i 









IBM 

JOB NO JOB TITLE 



1620/ 131 1 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 

ANALYST DATE PAGE OFJ^L 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 








COMMENTS 


1 . i . i i 


7| . . , . 


•3| . . . . 


19, , , , , 24 


31 1 1 i 36 . . . . I 41 . 




|56 16 1 166 .73 80 
1 1 1 1 1 




ftI.W.5. . 


ViA.L.U.E.. 


.... ...1.5, . . , 








1 1 1 1 1 




M.N.S. , 


r,L.£.L,D. 


1 ■ • .»■• 


il ,5.0,0. 


1 .......... . 


1 1 1 1 1 




BiT.N.S, , 


F,C, . . . 


1 . . . ...O.2. , . . 


1 . 





1 1 1 1 1 




BlI.N.5. . 


CiU, , , , 


.,0.6>, , . . 








1 1 1 1 1 




BlI.KJ.5. , 


MiM , , , 


..0.2. . . . 





1 , , , ,. 


1 1 1 1 1 




BiI.N.5. . 


, , . 


1 .0.1. . . . 








1 1 1 1 1 




&1.N.S, , 


AiL. . . . 


■ ro. . . . 


, , , . , 


1 , , , ... 


1 1 1 1 1 




BiI,N.5. . 


Sil. . . . 


1 . . . ...0.2. . . . 





1 , , . , • 1 , 1 , 1 1 


1 1 1 1 1 


i i ■ ■ i 


1 . . . . 


1 . . , , 


■ .... - 





1 1 1 > 1 > 


1 1 1 1 1 




AiL.U.M. . 


ViA.L.UE. 


1 2. 1. . . , 


1 . , . ,., 


1 


1 1 1 1 1 





kii.v.n , 




1 . , J 


1 . 


1 . . , ,. 


1 1 1 I I 




AiL.U.M. , 


F|C. , , . 


0.1. . . . 


1 . 




1 1 1 1 1 




AA.U.M. , 


CI", , , , 


0.1. . . , 


I . 


1 1 , , , * 1 1 1 1 , ■ 


1 1 1 1 1 






A|L, , , , 


1 . . . ...9.7. . . , 




1 ' . 1 1 ■ 1 ■ 1 








S,I. . . . 


0.1. . . . 




I'.''* 








1 . . . . 












1 . . . . 


1 . . . . 












1 . . . . 


1 . . . 1 












1 . , . . 


1 . . . ■ 












1 . . . . 


1 . . . ■ 












1 1 1 1 1 


1 . . . , 












1 . . . . 


1 . . . . 














1 1 1 i 1 . 












1 . 1 . . 


1 1 1 1 1 














1 . . . . 












1 . , . . 


I , . . 1 












1 1 ■ . . 


iii., 












1 . . . . 


1 . 1 1 1 











IBM 

JOB NO JOB TITLE 

KEYPUNCH- PUNCH DEC. PT. IN COL. 24 



1620/1311 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 

ANALYST _ DATE PAGE JLoFil 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






COMMENTS 


1 i i ■ i i 


-J I i I I 


,3 I , , . . 


191,,, , 2 . 4 


31, , _ |36 i f _ _ (41 _ 


I 46 . J 51 
— *■ — 1 — 1 — » — 1 — 1 — 1 — I — < — . 1 


|56 |6 1 _ |66 |73 80 


..... 


Sll.L.C.G.N 


ViA.L.U.E. 


.... ,.,3.9. . . . 












SiI.L.C.p.N 


Y,I.EL,D. 







1 




, , , . . 


S.I.L.C.0.N 


F.E . . 


1 . . . ...03. . . . 


| , , , ,. 






..... 


Sii.L.c.d.k 


S.I. . . . 


.... ...9,7. . . . 









.... i 


III.. 


1 . . . , 










1 1 1 ' ' 


. i . . i 


1 1 1 1 1 


1 . . . . 


, 





1 


1 1 1 1 ! 


. . . i . 


III.. 


1 . . . . 


, , , , , 


', 





■ 1 1 1 1 


. i . . . 


1 I I I I 


1 , , . . 





| , , , ,. 


, 


1 1 1 1 1 





1 . . . . 


1 . . . , 


I , , , , , 


1 1 1 1 1 * 1 1 1 ( 1 1 





1 1 1 1 1 





1 t 1 1 1 


1 . . . . 


, 







11 1 1 1 





1 1 1 1 1 


1 . . . . 


1 , , , , 


1 . 1 1 1 * 1 . , 1 1 1 





1 1 1 1 1 


. . . r . 


1 1 1 1 I 


1 . . . . 




1 1 1 . t • 1 1 . . . . 





■ 1 1 1 1 





—1 1 1 1 1 


1 , . . , 


, , , , , 


lit . 1 * 1 1 1 1 .1 





1 1 1 1 1 





—1 1 1 1 1 


1 . . . . 


1 . . ■ . * 1 . . . ■ . 







1 1 . . • 




1 1 1 1 I 




1 1 . . . . , 










_I 1 1 1 1 




1 . . . . * . 1 1 , , . 


1 1 1 1 1 ' 1 








1 ' 1 ' 1 




1 1 1 1 ... 1 . . 1 . 


1 . 1 1 1 • .. 1 1 1 . 1 








1 . . 1 1 


















■ ■ . 1 ' * ■ 
















1 ' > ' ' • 












1 , . . . • . . . . . . 














1 . . . . ... . 1 . . 














ii..'* 














1 . . 1 . • . . 1 . 1 . 






































1 1 1 1 1 






1 1 . . 1 • 1 1 1 . 1 . 














— 1 1 1 1 1 * 1 1 1 1 1 1— 







IBM 

JOB NO JOB TITLE 

KEY PUNCH: PUNCH DEC. PT. IN COL. 24 



IBM 



1620/ 13 1 1 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 



JOB NO.. 



JOB TITLE 



.ANALYST. 



DATE. 



PAGE_E_OF.iL 



KEYPUNCH: PUNCH DEC. PT. IN COL. 24 



to 
i— 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






COMMENTS 


i I I 1 1 


7 I . i i i 


'31 , . . , 


191,,, , 2 . 4 , , , 


31 1 , . . t 36 , ... I 41 , 


. . I 46 . . . . I"'. . . 


■ 56 .61 .66 ,73 80 


F.l.R.S.T,. 


B|A.L,«J,Y,1 


I.I.I 


1 ■ ■ 1 ■'■ l ■ » 1 . 


i ■ ■ . ■ • 


1 ■ ■ , 


1 1 • J 1 


i i i i i 


1 ■ ■ ■ ■ 


YiI.e.l.o. 


i2,0.0,O. • 





• 


1 1 1 1 1 


i i i i i 


1 i ■ ■ > 


FiE, . . . 


i . ,6.0.., , , , . , 


1 I 1 I I'l 1 I 1 1 I 


I....' 


■ 1 1 1 1 




1 ■ ■ ■ ■ 


C.U. . , . 


i ,1.0.0.. 


* 


1 . . . 


1 1 1 1 1 


i i i it, 




MlN, iii 


i , ,f,o,., , . , , . 


'■■■■'■I'll. 


1 l I I I*. I , I I . 


1 1 1 1 1 


■ i i i i 




i . i 


i . .3,0,. 


!■■'■• 


Ill II* 


1 1 1 1 1 





1 ■ ■ ■ ■ 


AiL, . . . 


ll ,5,0.0. 


1 


1 . I , 


1 1 1 1 1 


■ i i i i 


1 ■ i ■ i 


Sil, . , . 


I .3,0,0 


i .?.S.Q. 


1 I . 1 1*1 I I . 1 I 


1 1 1 1 1 


E.M.P.A.T.A 


1 ■ ' 1 1 


i , . . i 


' ■ ■ ■ ■ 


' ■ > ■ ■ ■ ' 


» 1 I I I 1 ■ 


1 1 1 1 1 


i i i i i 


1 ■ ■ ' ■ 


i . i i . 


1 ■ ■ ■ ■ ■ ' 


1 ■ ■ ■ ■ ■ 


1 ■ ■ ■ I 


- 1 1 1 1 1 





1 ■ ■ ■ ■ 


i i , , . 


1 1 ■ 1 ■•' 1 1 1 1 ■ 


1 


I ■■■■'■■■■ 1 ■ 


1 1 1 — 1 u 


i i i i i 


1 ' 1 ■ 1 


i i i i i 


1 i i ■ i ■ ■ 


1 1 


11111 


1 1 1 1 1_ 


■ i i i i 


11111 


i i i i i 


111 


1 1 1 1 '•' 1 1 1 1 1 




1 1 1 1 u 






i i . . i 










■ i i i i 




i . . , . 















i , . , , 










■ i i i i 




i . i . i 










i i i i i 




i i . . . 














i ... i 














i . i . i 










■ i i i i 














i i i i i 




i . . . . 










> i i i i 




i i . i i 










i i i i i 




i , , , . 














i . . . . 














i , , i i 








1 1 1 1 1 


■ i i i i 




i . i i i 








1 1 1 1 1 






i . , , . 








1 1 1 1 1 



Let us review each type of card shown (the b is used to represent a blank or space): 



ROW. ID This is the ROW. ID indicator card, which precedes the row 

identification cards. 

bbbbbbbbbbbb VALUE 
bbbbbbbbbbbO YIELD 
bbbbbbbbbbb+FE 



These are the row identification cards. One is required for 
each row. In column 12 is a code which indicates the type of 
row: 

b (blank) means objective function 
(zero) means an equality 

+ (plus) means a less-than inequality or a range constraint 
- (minus) means a greater-than inequality 

To determine the row type, look at the right-hand side of the matrix shown earlier. If there 
is no number, then it is an objective function (put a blank in column 12). If there is an equal 
sign, put a zero in column 12. If there is a less-than-or-equal-to sign (<), put a plus (12 
punch) in column 12. If there is a greater- than-or-equal-to sign (>), put a minus (11 punch) 
in column 12. If there is a range of values (A to B), put a plus (12 punch) in column 12. 

The name of the row is also given (e. g. , VALUE, YIELD, etc.). The name is placed in a 
six-position alphameric field with all special symbols permitted except the record mark and 
group mark. 

COL. ID This is the COL. ID indicator card, which precedes the column 

identification cards. If the column identification file is omitted, 
it will be produced automatically by the 1620-1311 LP System. 

bbbbbbbbbbbbBINl 



bbbbbbbbbbb-AL 
bbbbbbbbbbb+SI 



bbbbbbbbbbbbSILCON These are the column identification cards. One is required for 

each column. In the twelfth position on the card may be a code 
(an asterisk) to indicate that this structural variable should not 
be included as an active item in the solution of the problem. The 
name occupies a six-position alphameric field with all special 
symbols permitted except record mark and groupmark. 

MATRIX This is the matrix indicator card, which precedes the matrix 

element cards that contain the data values and the bounds of the 
structural variables. 



22 



bbbbbbBINlbbVALUEb +bbbb. 03bbbb 
bbbbbbBINlbbYIELDb+bbbl. bbbbbb+b200. bbbbbb 



bbbbbbSILCONSIbbbb+bbbb. 97bbbb 

These are the matrix element cards. Each contains the data 
value for one row-column intersection in the matrix. This 
data value is the coefficient for one term in one equation or 
inequality. The matrix element cards must be grouped in 
column sequence; that is, all data values for a given column 
must be together (e. g. , all rows for BIN1, then all rows for 
BIN2, finally all rows for SILCON). The row order within a 
column is not significant and need not be consistent from col- 
umn to column (the row names are VALUE, YIELD, (. . . .), 
SI). Obviously, it is safer to record the columns (and rows 
within the columns) in the same order as in the matrix itself. 
Only nonzero entries need to be recorded; however, if the 
user feels he may need to insert a data value at some later 
time for this row-column intersection, he should prepare a 
matrix element card with a value of zero. The row names 
and column names must correspond exactly to the names used 
in the row identification cards and in the column identification 
cards. If no COL. ID indicator and column identification 
cards were used, then the column names which appear in the 
matrix element cards will be assumed to be the correct 
names, and a column identification file will be constructed 
from these names with all of the structural variables con- 
sidered as active. 

Note the recommended format for the data values: 
ixxxx.xxxxxx. Either blanks or zeros may be used when the 
number does not require all the positions, though blanks are 
more usual. The plus sign is not required for positive num- 
bers, but the minus sign is required for negative numbers. 
Integers do not require a decimal point. In the "Matrix Ele- 
ment" subsection of "Data Preparation", there are more detailed 
rules on additional flexibility permitted in data value format. 

The matrix element cards have one additional function: they 
are used to specify the upper and lower bounds of the struc- 
tural variables. These bounds must both be on the same 
matrix element card and only one matrix element card should 
contain the bounds. Many users will find it convenient to 
consistently use the first (or last) matrix element card for a 
column to specify the appropriate structural bounds. The 
format for showing the bounds is the same as for the data 
values: ±xxxx.xxxxxx with the same construction rules. The 
upper bound, if any, appears in position 31-42 while the 
lower bound is in positions 43 - 54. For instance, in this 
matrix element card: 

bbbbbbBIN3bbYIELDb+bbbl. bbbbbb+b800. bbbbbb+b400. bbbbbb 

BIN3 is the column name, YIELD is the row name, the 
matrix element value is 1. , the upper bound is 800. , and the 
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lower bound is 400. Unless specified, all lower bounds are 
implicity zero. Where proper (because of the nature of the 
equalities or inequalities), either positive or negative bounds 
may be used (1620 - 1311 LP actually transforms these state- 
ments so that the lower bound is effectively zero during cal- 
culation and is adjusted correctly before reports are prepared) . 
The upper bound must be algebraically larger than the lower 
bound. If a negative upper bound is used, a lower bound is 
required. 



FIRST. B ALOY1 This is the indicator card for the first right-hand side. Its 

name is ALOY1 . The FIRST. B indicator card precedes the 
right-hand-side element cards. If there is only one right- 
hand side it may be given a name of all blanks (bbbbb). 

bbbbbbbbbbbbYIELDb+2000. bbbbbb 
bbbbbbbbbbbbFEbbbb+bb60. bbbbbb 



bbbbbbbbbbbbSIbbbb+b250. bbbbbb +b300. bbbbbb 

These are the right-hand-side element cards. There is one 
for each row, except for those rows with a zero limit (or no 
limit, as in an objective function row). The row names are 
YIELD, FE, . . . , SI. If there is only one limit, it is placed 
in positions 19-30. If there is a range limit, the upper con- 
straint is placed in positions 19-30 and the lower constraint 
is placed in positions 31-42. The format for these limits is 
+xxxx.xxxxxx, just like the matrix element data values and 
bounds. The first value may be an equality constraint or a 
lower limit or an upper limit, if there is only one constraint 
limit. The relationship (equality, less than, greater than) is 
specified through the row identification card by the row type 
code. 



ENDATA This is the data termination indicator card; it concludes the 

data input. 

Now that the principal data inputs have been covered, this information can be placed into 
disk storage. The following instructions (agenda and data) store the data ready for problem 
solution: 



**XEQbLP1620 
INPUT. C02000ALLOY 

ROW. ID indicator 
row identification (8 cards) 
COL. ID indicator 
column identification (7 cards) 
/ MATRIX indicator 
matrix elements (48 cards) 
FIRST. B indicator 
right-hand-side elements (7 cards) 
ENDATA 

END JOB 
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Each of the lines on the LP form is keypunched into an individual card. These cards are 
inserted in the card reader of the 1620-1311 system. The start button is depressed and 
three minutes later the eight-row-by-eight-column matrix is stored on disk ready for later 
use. 

Whenever the entire problem data must be read in, a substantial number of cards is 
required. This is why so much emphasis is placed on disk maintenance of problem data. 
Each input card is now examined: 

4=£XEQbLP1620 This is the 1620 Monitor program execute card, which 

instructs the 1620 Monitor to call in from disk storage the 
Linear Programming Initialization routine. This routine 
sets up a common communication area for use by the various 
LP routines. The Initialization routine in turn calls the 
Linear Programming Supervisor routine, which is the con- 
trol or executive program for all other LP routines. The LP 
Supervisor is then ready to read another card. The exact 
form shown must always be used to initiate an LP solution. 
The ++XEQb indicates a program execute instruction (in the 
1620 Monitor language) where LP1620 is the established 
name for the LP program. The 1620 Monitor has a record of 
where the LP program is stored. 

INPUT. C02000ALLOY This is the first LP agendum card. It serves to call in the 

routine which will set up the input data for problem solution. 
INPUT, is the agendum name; all agendum names are six 
characters long, so the period is used to make this name fill 
out the full space allotted. The C indicates that the data for 
this problem is on cards: it is to be stored on disk starting 
at location 02000 (a five-position numeric field). The name 
of the problem stored at 02000 is ALLOY; the problem name 
occupies a six-position alphameric field. 

data input deck These 75 cards contain the entire problem data which is to be 

stored for later solution. 

ENDJOB The final agendum card indicates to the LP Supervisor 

routine that the LP problem has now been solved, that the 
appropriate internal storage areas should be cleared and that 
control should be returned to the 1620 Monitor program. The 
monitor then is ready to initiate processing of whatever job 
is next in the card reader. The problem data, of course, 
remains in disk storage for use at some later time. This 
termination card is used at the end of all LP processing; it 
should not be used between different LP jobs. 

SAMPLE PROBLEM 2 

The LP program has been stored on a 1311 under control of the 1620 Monitor. The data and 
control information for the particular LP problem has also been stored starting at disk loca- 
tion 02000. This problem is concerned with blending certain raw materials to produce a 
desired end product. On disk have been stored (1) the data concerning the available inven- 
tories of the raw materials and their composition, (2) the demand for the end product, and 
its composition constraints, and (3) the cost of each available raw material (ingredient). 
The purpose now is to solve this problem using linear programming: to determine the best 
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(optimum) combination of ingredients to use, to produce the end product, within the stated 
constraints at a minimum cost. 

The following five instructions will accomplish this goal: 

4=4=XEQbLP1620 
INPUT. D02000ALLOY 
MIN. . . 

OUTPUTbbbbbbbbbbbb23 
ENDJOB 

*4=XEQbLP1620 This is the 1620 Monitor program execute card, which brings 

into the main memory the first LP routine. 

INPUT. D02000ALLOY The INPUT, agendum card indicates that the data for this 

problem has already been stored on disk starting at location 02000. The original 
problem name was ALLOY. 

MIN. . . This agendum card calls in the routines which will develop an 

optimum solution to the LP problem stored at 02000. The solution is a two-step 
process: first, a feasible solution is found, e. g. , one in which all the constraints 
and bounds are satisfied; second, the optimal solution is determined, e. g. , one 
which is feasible and for which the costs are minimized. Costs may be con- 
sidered in dollars, process time, or any other suitable measure of value. 

OUTPUTbbbbbbbbbbbb23 The third agendum card calls in the output reporting routine 
from disk. This will cause appropriate output cards to be punched so that a 
report on the solution to the LP problem can be printed. These output cards state 
specifically how much of each raw material is used in producing the desired end 
product. The number 2 indicates that these answers will be given to two decimal 
places, e. g. , ±xxxx.xx . In addition to the raw materials usage (those materials 
that are used are called the basic structural variables), the output routine also 
produces a card for each constraint row indicating the amount by which the end 
product differs from the constraint limitation on that row; e. g. , if the amount of 
iron must be less than 60 pounds out of 2000 pounds of end product, and only 40 
pounds were used, then there would be a card which would note that the end pro- 
duct contained 20 pounds less than the allowed limit. The number 3 on the 
agendum card means that these values will be shown to three decimal places, e. g. , 
ixxxx.xxx . A card is also punched for the objective function row giving the cost 
of the optimum solution. 

ENDJOB This agendum card turns control back to the 1620 Monitor. 

For this simple eight-row-by-eight-column example the 1620 LP System produces 17 cards 
stating the solution to the stored problem. The cards are so arranged as to provide appro- 
priate headers for the columns of data. The answers can be listed with a regular 80-column 
print board on a 407 or by a listing program on a 1440, 1401 or other appropriate machine. 
Through use of an ASSIGN agendum card (see "1620-1311 LP Agenda") the answers can be 
printed on the console typewriter (practical only for very small problems) or listed by the 
1620 on an online printer. 

While the problem is being solved, certain information is printed at the console typewriter. 
The 1620 Monitor notes operating data, while the LP supervisor lists each agendum card 
and key internal LP operations, like assignment of mantissa length and tolerances, use of 
the optimization and inversion routines, and, of course, minor or major errors encountered. 



26 



The agenda which have been covered in Sample Problem 2 are the fundamental agenda needed 
to prepare, solve, and present the results of an LP problem. While many special options 
and additional agenda are still to be covered, the routines called in solving problem 2 are 
the workhorses. These routines are always required to solve an LP problem. 



SAMPLE PROBLEM 3 



In Sample Problem 3 the matrix element values to be used are stored on disk from a pre- 
vious problem. The row names are used as an index to get at the matrix information effi- 
ciently. A new alloy is to be made which has different constraints from the previous alloy. 
It is desired to find the minimum cost (solution) for making this new alloy. It is also desir- 
able to be in a position to rapidly reanalyze the solution with changes in costs or constraints. 
The following instructions will solve this problem. 



* *XEQbLP1620 
INPUT. C03000ALLOYA 



data 
input 



ROW. IDD02000ALLOYA 
MATRIXD02000 
FIRST. BALOY2 

b b VALUE 

b b YIE LDb-D2000bbbbbb 

b b FEbbbb-Dbb60 



cu 


140 


MN 


40 


MG 


40 


AL 


1500 


SI 


300 



bbb250 



ENDATA 



MIN ALOY2 VALUE 

OUTPUTbbbbbbbbbbbb23 
SAVE. B. 10000 ALLOYX 
ENDJOB 



Quite a few cards are required in this problem since some of the input data is being specified 
explicitly rather than just through reference to previously stored information on disk. 
Notice that agendum cards are combined in the same deck with data cards and that the order 
of the cards is significant. Let us review each card to note its format and usage. 



**XEQbLP1620 



INPUT. C03000ALLOYA 



The same 1620 Monitor program execute card which initiates 
the LP system. 

The INPUT . agendum card calls in the input routine. In this 
case the C indicates that some input information is on cards 
and will directly follow the INPUT, agendum card. The user 
wants the problem stored beginning at location 03000 and the 
problem name to be assigned is ALLOYA. 



ROW. IDD02000ALLOY 



Following an INPUT. C agendum card there will be data input 
cards. Since different kinds of information can be inserted, it 
is necessary to have an indicator card to tell what kind of data 
will follow. ROW. ID is the name of an indicator card which 
precedes the row identification cards (i. e. , the cards containing 
the names of the objective function and constraint rows). In this 
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case the D02000 indicates that the row names are already stored 
on disk at location 02000; the problem name where they are 
stored is ALLOY. The name will be used as a double check 
before the new problem picks up the information from the refer- 
enced problem. The row names are used to index the matrix 
element values and right-hand-side element values. 

MATRIXD02000 ALLOY This is another input indicator card. MATRIX is the card name. 

The D02000ALLOY indicates that the matrix element values are 
already stored on disk at location 02000 and the referenced 
problem name is ALLOY. The matrix element values refer- 
enced (e. g. , the cost coefficients and the raw material compo- 
sitions) are those that will be used in solving the new problem. 
The information is not recopied; only appropriate linkages are 
made to the data. 

FIRST. BALOY2 The third indicator card is called FIRST. B . This is the card 

which introduces the right-hand-side element values. The B 
has been used historically to refer to these right-hand-side 
values; they are frequently called B values. ALOY2 is the 
name given to this special column. The right-hand-side name 
must be only five alphameric characters. 

b b VALUE, etc. These are the right-hand-side element cards. They give the 

actual constraint limits for each row in the matrix. The row 
name in the first card is VALUE: a six-position alphameric 
name is permitted. In the second card the name is YIELD. 
The 2000 indicates that the row called YIELD has a constraint 
limit of 2000. The row YIELD takes the following form: 

BIN1 + BIN2 + BIN3 + BIN4 + 
BIN5 + ALUM + SILCON = 2000 

So in this case the limit is an equality; i. e. , the sum of the 
weights used of each raw material must equal the required 
weight of the end product. 

The third through sixth rows (FE, CU, MN, and MG) also have 
constraint limits (60, 140, 40, 40 respectively). The constraint 
formula for FE takes this form: 

. 15(BIN1) + . 04(BIN2) + . 02(BIN3) + . 04(BIN4) 
+ . 02(BIN5) + . 01(ALUM) + . 03(SILCON) ^ 60 

The 60 in this case is an upper constraint limit. The same is 
true for the next three rows covered (CU, MN and MG). 

The next row, AL, states a lower constraint limit, 1500; i. e. , 
the amount of aluminum must equal or exceed 1500 pounds. 

Notice that regardless of whether the constraint limit is an 
equality, an upper constraint limit, or a lower constraint limit, 
the value is always placed in the same place on the right-hand- 
side element card. 
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In the last row there is a different situation. For SI , two 
values are given: 300 and 250. The 300 represents the upper 
constraint limit for SI, while 250 represents the lower con- 
straint limit. This is called a range constraint. In the rows 
with upper limits a lower limit of zero is implied. If there is 
an explicit upper and lower limit, the upper limit is placed first 
and the lower limit second. 

The field used for stating the limits is twelve positions wide. 
It is recommended that the following format be consistently 
followed: 

ixxxx.xxxxxx 

If the value is positive, no sign is required although a plus sign 
is often useful. For integer values the decimal point is not 
needed, though it too is often used. This arrangement tends to 
keep the values in a reasonable range for most efficient LP 
solution. 

ENDATA The indicator card which states the input data is now concluded. 

MIN ALOY2 VALUE MIN. . . is the agendum name which calls in the routine to 

determine the optimum mix of raw materials. In this case 
ALOY2 refers to the right-hand-side name of the values being 
used. If there were different alloys being smelted, it would be 
necessary to state which particular set of constraint limits 
should apply to the current problem solution. VALUE is the 
name of the objective function row used in evaluating the various 
possible solutions to pick the optimum. Since alternate objec- 
tive function rows may be used to state different cost criteria, 
the particular objective function row to be used in this problem 
solution must be named. 

OUTPUTb b23 The OUTPUT agendum card produces the same types of cards 

as indicated in Sample Problem 1: a card for each structural 
variable and one for each row indicating the difference between 
total usage and the constraint limits (to two decimal places). 

SAVE. B. 10000ALLOYX This is another agendum card. It is called SAVE. B to note that 

the user wishes to retain the current solution information in 
disk storage for possible later use. This solution is known as a 
basis; if relatively minor revisions are made to the matrix ele- 
ments or to the right-hand-side values or to the objective func- 
tion, a previous basis can provide a head start in solving the 
new problem. 10000 is the location where the problem informa- 
tion and the basis is to be stored, and ALLOYX is the name 
given to this stored version of the problem. 

Essentially, the output data is stored to form a possible starting 
basis. 

ENDJOB This agendum card concludes the LP solution process and re- 

turns control to the 1620 Monitor. 
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This sample problem has shown how some data from a previous 
problem may be used in conjunction with new constraint limits 
to set up a new problem. It has also indicated how data is 
entered into disk storage for later use. 

SAMPLE PROBLEM 4 

Usually the LP user is working with a previous problem, but has to make revisions to 
reflect changing inventories, costs, material composition, or constraints. 1620-1311 LP 
has been designed to allow changes to be made readily. The fourth sample problem illus- 
trates how some of these changes can be introduced. 

After analyzing the ALLOYA solution from Sample Problem 3, the aluminum smelter wants 
to determine whether it will pay to purchase 2000 pounds at $. 19 per pound or 2500 pounds 
at $. 18 per pound of a new scrap raw material; he is given its composition. He wishes to 
consider producing the original alloy (from Sample Problem 2) using the new material as a 
substitute for BIN5. 

He will examine the cost and usage of the new material. The problem is solved with the 
following set of input cards: 

*+XEQbLP1620 
INPUT. D02000ALLOY 
REVISE 
MATRIX 

b bBIN5bbVALUEbbbbbb. 19bbbbbb2000 

BIN5 FE . 05 

BIN5 CU . 08 

BIN5 MN .01 
BIN5 MG . 01 

BIN5 AL .84 
BIN5 SI .01 

ENDATA 

MIN ALOY1 VALUE 

OUTPUTb b23 

SAVE. B 

REVISE 

MATRIX 

b bBIN5bbVALUEbbbbb. 18bbbb2500 
ENDATA 

MIN ALOY1 VALUE 

OUTPUTb b23 
ERASEBD02000ALLOY 
END JOB 

**XEQbLP1620 The 1620 Monitor program execute card for 1620-1311 LP. 

INPUT. D02000ALLOY The INPUT, agendum card sets up the problem which is stored 

on disk at 02000 for re-solution. 

REVISE The agendum card, REVISE, serves to bring in a routine which 

permits the introduction of changes to the information stored 
for a particular problem. Only information in storage can be 
changed. No new information (new rows or columns, or addi- 
tional elements) can be added. REVISE must come only after 
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an INPUT, 
solution. 



It always revises the problem currently set up for 



MATRIX 



The MATRIX indicator card introduces the changed matrix 
element values. Only elements already in disk storage can be 
modified. 



b BIN5bbVALUEbbbbbb. 19bbbbbbbbb2000, etc. 



These are the matrix element cards. BIN5 is the structural 
variable. VALUE is the name of the objective function row. 
The element value is .19. The upper bound on this structural 
variable is 2000, i. e. , the quantity available to be purchased. 
If there were a lower bound it would follow the upper bound. 
The element value and the upper and lower bounds are all shown 
in the following form: 

±xxxx.xxxxxx 

The bounds for a single structural variable are shown on one 
matrix element card. 

The other matrix element cards state the composition of the 
proposed raw material; one card is required for each nonzero 
entry in a constraint row. 



END AT A 



This indicator card ends the data input for the REVISE agendum. 



MIN ALOY1 VALUE 



The MIN. . . agendum card optimizes the revised problem using 
the right-hand side called ALOY1 and the objective function 
called VALUE. The optimization starts from a previously 
saved basis to speed up the solution. This is not indicated 
explicitly, but MIN. , . will always start from any basis stored 
with the problem. 



OUTPUTb 



b23 



The OUTPUT agendum prepares the reports on raw material 
usage and total costs. 



SAVE. B 



The SAVE . B agendum stores the new solution basis with the 
current problem for later revise. 



REVISE 



The REVISE agendum is used again to make further changes in 
the actual material values. This is done in order to compare 
the usage of the BIN5 material at $.18 per pound and an avail- 
ability of 2500 pounds with the 2000 pounds at $. 19 per pound. 



MATRIX 



The MATRIX indicator card is used to introduce the value 
change for the matrix elements. 



BIN5bbVALUEbbbbb. 18bbbbbb2500 



The matrix element card changes the data value for the BIN5 
column in the row named VALUE (the objective function row); 
the data value is set at $. 18 pound. The upper bound is also 
changed for this structural variable to show an availability of 
2500 pounds. 
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ENDATA 



The data input termination indicator card ends the revised data. 



MIN ALOY. 1 VALUE 



OUTPUT. b_ 



Jb23 



ERASEB .02000 ALLOY 



END JOB 



The revised matrix is now optimized. The same right-hand 
side, ALOY1, and the same objective function, VALUE, are 
used. However, since the cost in the objective function row has 
been changed, it is possible that a different answer will result. 

The OUTPUT agendum prepares another report on raw material 
usage. The cost of this solution will be compared with the pre- 
vious solution to indicate the advisability of purchasing either 
2000 or 2500 pounds of the new raw material that is available. 
These results will also be compared with the previous solution 
(with the old BIN5 material) to decide whether the raw material 
should be purchased at all. 

The ERASEB agendum is the means by which a previous solu- 
tion basis can be eliminated from disk storage. Any inverse 
that has been stored will of course also be deleted. This is 
done only when the user expects that there will be such a signif- 
icant change in the right-hand-side or matrix elements that it 
would be better to start over again rather than use the previous 
results. If the basis is erased, the user cannot insert a 
revised basis. He would have to go through a normal input pro- 
cedure in order to get an initial basis. 

This agendum card ends the solution of this program and 
returns control to the 1620 Monitor. 



SAMPLE PROBLEM 5 

Often the most valuable aspect of using a good LP system is the ability to gain substantial 
insight into alternate solutions or sensitivity of the solution to changes in cost, composition 
or constraints. For efficiency, this type of analysis (sometimes called parametric pro- 
gramming) can be handled through special agenda cards. These cards can be used only 
after a MIN. . . or MAX. . . has been carried out. 

Let us assume that the alloy smelter decided to buy 2, 500 pounds of the new raw material 
at $ . 18 per pound. He now wants to examine in depth the use of this new raw material. 
The following instructions are one approach to this problem: 

**XEQbLP1620 

INPUT. D 02000 ALLOY 

MIN... 

COST. R 

DO. D/J 

CHECK 

ENDJOB 

Each of these cards will be reviewed: 

*+XEQbLP1620 This is the 1620 Monitor program execute card. 



INPUT. D02000ALLOY 



The input agendum card sets up the problem stored at 02000 
(with the problem name of ALLOY) for solution. 
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MIN The minimum value solution is determined. 

COST. R This agendum card calls in the program which determines the 

cost ranges for the structural variables which are in the optimal 
solution (called the optimal basis). For each such variable, 
COST.R determines the highest and lowest values of the objec- 
tive function coefficient which do not require any change in the 
solution activity level. In other words, this is the range over 
which the cost of the structural variable may change without 
changing the solution. The report also indicates what would 
happen if the cost of each variable were raised or lowered so 
as to be outside the stated cost range; it indicates which other 
variable would be affected by the change. This report, in 
effect, tells the user that as long as the costs remain within 
the stated boundaries he need not resolve the problem since he 
knows that the present answer will still be satisfactory. 

DO.D/J This agendum card calls in a program which will provide infor- 

mation on the economic effect of changes in the bounds and 
limits on certain structural variables and the rows. For each 
structural variable in the solution at its upper bound, it indi- 
cates how much the objective function value could be improved 
if the upper bound were increased one unit. It also indicates 
for each structural variable in the solution at its lower bound 
how much the objective function value could be improved if the 
lower bound were reduced one unit. For each row, the report 
indicates the additional profit to be obtained by increasing the 
upper constraint limit; similarly, the report states the value of 
decreasing the lower constraint limit. The DO.D/J report 
gives the user substantial additional insight into the sensitivity 
of the solution to changes in the bounds and limits. 

CHECK. This agendum card calls in a program which checks the accu- 

racy of the output which has been obtained. It does this by 
comparing the calculated sum for each row (coefficients times 
activity levels) and determining the difference between this total 
and the constraint limit. The largest row error is identified. 
This routine is used automatically in computing to make sure 
that error limits are not being exceeded. It may also be called 
for by the user to assist him in judging whether or not restart- 
ing or rescaling some rows may not provide an accurate solu- 
tion in less time. 

ENDJOB This agendum card indicates the end of the LP work and turns 

control back to the 1620 Monitor. 
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DATA PREPARATION 



The 1620-1311 Linear Programming System expects two input streams — the control stream 
and the data stream. The control stream contains agendum cards which direct the LP sys- 
tem in processing a particular problem. The data stream contains the data values to be 
operated on. Both types of input are introduced in the card reader. 

There are three main types of data input cards: 

• Indicator cards that announce the type of data following 

• Identification cards that either state the names of rows and columns used in the particu- 
lar problem or indicate the structural variable to be incorporated in the initial solution 
basis 

• Actual data value cards that introduce the values of matrix elements and right-hand-side 
elements 

The method of presentation is to describe the format and usage of each indicator card and 
then the identification cards and data value cards which follow it. The indicator cards are 
described in the recommended input sequence. 

To solve a linear programming problem, it is required that at least ROW. ID, MATRIX, 
FIRST. B, and ENDATA (or EOF) cards are all present in that order. The other indicator 
cards are optional. If the problem is concerned only with preparing output from a previously 
solved problem, then only a ROW. ID or a COL. ID and an ENDATA card are needed. 

The following states the sequence in which the data input cards are discussed: 

ROW. ID indicator 

Row identification 

COL. ID indicator 

Column identification 

MATRIX indicator 

Matrix elements 

FIRST. B and NEXT.B indicators 
Right-hand-side elements 
BASIS, indicator 
VARBLS indicator 

Basic structural variables identification 
SLACKS indicator 

Nonbasic logical variables identification 



ENDATA or EOF indicator 
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The 1620-1311 LP data input forms have been designed to be compatible with many currently 
used linear programming systems: 

1. SHARE standard matrix element cards will be accepted. 

2. 7040/7044 Linear Programming System (7040-CO-08X) data input decks will be accepted 
with four exceptions: 

a. On row identification cards 7040/44 LP permits a blank or zero for type of con- 
straint for both equality rows and objective function rows. In 1620-1311 LP each 
objective function row must have a blank in col. 12, and each structural row with an 
equality constraint must have a zero in col. 12. 

b. 7040/44 LP permits up to five row names per row identification card; only one row 
name per card is permitted in 1620-1311 LP. 

c. The 7040/44 LP basis cards have a completely different format and cannot be 
accepted. 

d. Certain data input cards in 7040/44 LP are not acceptable to 1620-1311 LP: SUM, 
CNT, CURTAIN, PARTITION, and SLACK. 

3. 1410 Linear Programming System (1410-CO-01X, 06X, 07X, 09X) matrix element cards 
will be accepted. 

1620/1311 LP data input will be accepted by the 7040/7044 LP system if: 

a. The cost rows are first in the ROW. ID deck. 

b. No bounds are used. 

c. No COL. ID's are used. 

d. NoperiodappearsinROW.ID, FIRST. B, NEXT.B. 

e. The first characters of the row and column names are blank or numeric. 

f . The right- hand-side names are numeric, and in ascending, but not necessarily consecu- 
tive, sequence. 

g. No basis indicator or identification cards are used. 

Since the data values are stored in a compressed form in disk storage, it is necessary to 
identify all matrix elements and right-hand- side elements that are expected to be used in 
later solutions of the problem. For instance, if the problem may need to be solved using 
different (but not yet defined) objective functions, extra objective function rows should be 
inserted so that appropriate values can be established later — similarly with alternate right- 
hand sides or initially zero-valued matrix elements. 

Two of the agendum cards produce output cards which can be recycled through the system. 
These agenda are OUTPUT and GETOFF . Some of the information on these recycled cards 
is ignored by the INPUT, and REVISE agenda. 

For a more complete discussion on how the data input cards are used to set up a problem for 
linear programming solution, refer to INPUT, and REVISE under "1620-1311 LP Agenda". 

A comments card which appears within a data input deck will not be typed. 
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ROW. ID INDICATOR 



Format 



Indicator 
Name 


Input Data 
Source 


Starting Sector 
Location 


Problem 
Name 


Comments 


1 6 


7 


8 12 


13 18 


19 80 


ROW. ID 
ROW. ID 


D 


sssss 


nnnnnn 





Indicator Name: ROW. ID; the period in column 4 is optional. 

Input Data Source: b (blank) indicates that row identification cards will follow. 

D indicates that the row identification file should be referenced from 
a previous problem. 

Starting Sector Address: A five-position numeric field which identifies the established 
location for the problem containing the row identification information desired. 

Problem Name: A six-position alphameric field containing the problem name which is 
checked against the problem name stored at the referenced sector location. 

Comments: Any desired information 

Usage 

The ROW. ID indicator card introduces the row names which are on the row identification 
cards or in disk storage. 

Only one ROW. ID card may be used per INPUT, card. It is always required to follow an 
INPUT. C card if MIN. . . or MAX. . . are to be executed for the problem. The ROW.IDb 
card must always be followed by the appropriate row identification cards. If there is a 
MATRIX card there must be a ROW. ID card preceding it. If ROW.IDD is used, the 
current problem refers to information stored with a previous problem. If that previous 
problem is destroyed (through writing a new problem over it), the second problem will not 
be usable. 

ROW IDENTIFICATION 
Format 



Blank 


Row Type 


Row Name 


Comments 


1 11 


12 


13 18 


19 80 


b b 


t 


rrrrrr 
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Blank: Must be all blanks 



Row Type: The type of logical variable to be generated for this row: 
b (blank) for objective function 
(zero) for equation 

+ (12 punch) for a less-than-or-equal-to (<) relation or a row with both 
an upper and lower constraint. 

(11 punch) for a greater-than-or-equal-to (>) relation or a row with 
both an upper and lower constraint. 
* (asterisk) for an omitted row (it can be revised later) . 



Row Name: A six-position alphameric field with special characters (other than record 
mark and groupmark) permitted. Blanks are significant characters; therefore, it is 
recommended that a period be used for internal spacing, for example, OB. ROW, not 
OBbROW. The names should all be left- justified. 

Comments: Any desired information 



Usage 

The row identification cards are used to define the type of row relationship (equality, less 
than or equal to, etc.), thereby specifying the type of logical variable (slack or artificial) to 
be generated. They are also used to name the rows (both objective function and constraint) 
for user reference. The row identification cards permit the selection of rows from a matrix 
by asterisking the unwanted row identification cards; this will omit them from the matrix. 
The order of row entry is immaterial, and objective function rows may be interspersed with 
constraint rows. 



A row identification card may be omitted rather than specifying an asterisk for it. This 
avoids using unnecessary storage. However, a row which is not defined, rather than being 
defined with an asterisk as an omitted row, cannot be later entered under the REVISE 
agendum. 

One row identification card is used for each row in the full matrix; every row must have a 
row identification card containing its row name. Omission of a row identification card will 
cause a minor error message when a matrix element card is read in for that row. 

The row type is a method for having 1620-1311 LP generate the artificial and slack columns 
by creating logical variables. These make it an equality; they are called artificial variables 
when placed in the objective function or in a constraint equality. They are required to pro- 
vide the basis (one column for each row) needed to start the linear programming solution 
process. An artificial variable can never re-enter the basis once it has left. 



Row type 

objective function 

equality constraint 

less-than constraint 
[A 1 X 1 +... <B] 



Col. 12 

code 

b (blank) 
(zero) 

+ (12 punch) 



Logical variable generated 

a free artificial unrestricted as to sign 

a zero-valued artificial 



a positive slack 
[A 1 X 1 +...+S=B] 

where S^O 
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Row type 



Col. 12 

code 



Logical variable generated 



greater-than constraint 
[A^H-... >B] 



- (11 punch) 



a negative slack 
[A^X^. . . -S=B] 
where S ^ 



range constraint 



+ (12 punch) 
or 

- (11 punch) 



treated as a positive slack with upper bound 
specified 

treated as a negative slack with upper bound 
specified 



omitted row 



* (asterisk) no logical variable 



COL.ID INDICATOR 
Format 



Indicator 
Name 


Input Data 
Source 


Starting Sector 
Location 


Problem 
Name 


Comments 


1 6 


7 


8 12 


13 18 


19 80 


COL.ID 
COL.ID 


D 


sssss 


nnnnnn 





Indicator Name: COL.ID; the period in column 4 is optional. 

Input Data Source: b (blank) indicates that column identification cards will follow. 

D indicates that the column identification file should be referenced from 
a previous problem. 

Starting Sector Address: A five-position numeric field which identifies the established 
location for the problem containing the column identification information desired. 

Problem Name: A six-position alphameric field containing the problem name which is 
checked against the problem name stored at the referenced storage location. 

Comments: Any desired information 

Usage 

The COL.ID indicator card introduces the column names which are on the column identifica- 
tion cards or in disk storage. Only one COL.ID card may be used per INPUT, card. The 
column identification file (COL.ID indicator card and column identification cards) is strictly 
optional. It is required only if the user wishes to eliminate certain columns from the matrix 
for a particular problem solution. If there is a COL.ID indicator card, every column must 
have a column identification card. If no column file is introduced, the 1620-1311 LP system 
will make one using the names in the matrix element files. 
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If COL.IDD is used, the current problem refers to information stored with a previous prob- 
lem. If that previous problem is destroyed (through writing a new problem over it), the 
second problem will not be usable. 

COLUMN IDENTIFICATION 

Format 



Blank 


Column 


Column 


Comments 




Usage 


Name 




1 11 


12 


13 18 


19 80 


b b 


t 


cccccc 





Blank: Must be all blanks 

Column Usage: b (blank) indicates that the column is to be used; it is an active structural 
variable . 

* (asterisk) indicates that the column is not to be used; it is omitted from 
the solution; it can be revised later to be considered in the solution. 

Column Name: A six-position alphameric field with special characters (other than record 
mark and groupmark) permitted. Since blanks are significant characters, it is recom- 
mended that a period be used for internal spacing, for example, ST. COL, not STbCOL. 
For consistency and accuracy the names should be all left-justified. 

Comments: Any desired information 

Usage 

The column identification cards are used to select certain columns (or structural variables) 
which will be in the current matrix. This permits the solution of various problems from one 
large matrix. The order of column entry is immaterial. Column identification cards are 
also used to set for selective analysis reports. 

One column identification card is used for each column in the matrix; if the column identifica- 
tion file (COL. ID indicator card and column identification cards) is used, every column in the 
matrix must be represented. Omission of a column identification card will cause a minor 
error message at the time that a matrix element card is read in for that column. 

The use of the asterisk makes it convenient to omit certain columns. By replacing the 
asterisked card with a blank card, the column can be included in the matrix. 

Note that column identification cards are not needed for the automatically generated logical 
variables, nor are they used for the right-hand sides. 
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MATRIX INDICATOR 
Format 



Indicator 
Name 


Input Data 
Source 


Starting Sector 
Location 


Problem 
Name 


Comments 


1 6 


7 


8 12 


13 18 


19 80 


MATRIX 
MATRIX 


D 


sssss 


nnnnnn 





Indicator Name: MATRIX 

Input Data Source: b (blank) indicates that matrix element cards will follow 

D indicates that the matrix element data values are to be referenced 
from a previous problem 

Starting Sector Location: A five-position numeric field which identifies the established 
location for the problem containing the desired matrix element data values. 

Problem Name: A six-position alphameric field containing the problem name which is 
checked against the problem name stored at the referenced storage location. 

Comments: Any desired information 

Usage 

The MATRIX indicator card introduces the matrix data values which are on the matrix 
element cards or in disk storage. Only one MATRIX card may be used per INPUT, card. 
The MATRIX file (MATRIX indicator card and matrix element cards, or MATRIXD indicator 
card) is required if a problem solution is to be obtained. If MATRIXD is used, any change 
to the referenced matrix may make the current problem unsolvable. 

MATRIX ELEMENT 

Format 



Blank 


Column 


Row 


Data 


Upper 


Lower 


Comments 




Name 


Name 


Value 


Bound 


Bound 




1 6 


7 12 


13 18 


19 30 


31 42 


43 54 


55 80 


b b 


cccccc 


rrrrrr 


±xxxx .xxxxxx 


iyyyy-yyyyyy 


±zzzz .zzzzzz 




b b 


cccccc 


rrrrrr 


±xxxx .xxxxxx 


±yyyy-yyyyyy 






b b 


cccccc 


rrrrrr 


±XXXX .XXXXXX 




±zzzz.zzzzzz 




b b 


cccccc 


rrrrrr 


±XXXX .XXXXXX 









40 



Blank: Must be all blanks 



Column Name: A six- position alphameric field containing the column name. Blanks are 
significant characters. Record mark and groupmark cannot be used. It is recommended 
that the name be left-justified and that periods be used for internal spacing. 

Row Name: A six-position alphameric field containing the row name as included in the 
row identification file 

Data Value: A twelve-position field permitting a ten-position numeric value with space 

for a sign (optional for positive numbers) and a decimal point. The format shown is 
recommended since it is used for SHARE standard input and 7040/7044 LP system 
(7040-CO-08X). The recommended format aids in obtaining reasonable scaling of values 
for more accurate solutions. All blanks are interpreted as zeros. There is no checking 
for alphabetic information in any data field. The user may vary the construction of this 
field provided the following rules are observed: 

• One to ten contiguous digits with or without a decimal point. The decimal point may 
appear prior to the first digit, between two digits, or after the last digit. The 
decimal point is not required for an integer. 

• A minus sign (11 punch) is required for all negative numbers. The plus sign (12 
punch) is optional for positive numbers. The sign may appear anywhere before the 
first digit or after the last digit and need not be contiguous with the number. In an 
objective function row, cost data values are usually entered as positive numbers and 
profit data values entered as negative numbers, but they may be reversed if the user 
is careful to use MAX. . . instead of MIN. . . 

Upper Bound: A twelve-position field permitting a ten-position numeric value with space 
for a sign and a decimal point. Same rules for construction as in data value. Usage is 
optional. An asterisk in the first position will delete a previously stored value. 

Lower Bound: A twelve-position field permitting a ten-position numeric value with space 
for a sign and a decimal point. Same construction rules as for data value. A lower 
bound is required if there is a negative upper bound; otherwise usage is optional. An 
asterisk in the first position will delete a previously stored value. 

Comments: Any desired information 

Usage 

There are four variations of matrix element cards: data value only; data value and upper 
bound; data value and lower bound; data value and both upper and lower bound. Every matrix 
element with a nonzero data value (coefficient) must have a matrix element card. Every 
matrix element which the user may wish to have inserted in the matrix at a later time must 
have a matrix element card when the problem is set up. The principal purpose of the matrix 
element cards is to establish the coefficients for the matrix to be solved. To do this it is 
required that the matrix element cards be arranged in column name sequence; for example, 
all matrix elements cards for a given column must be consecutive. The order of the rows 
within the column is not significant. All constraint rows and objective function rows are 
entered. The right-hand sides are not entered on the matrix element cards. 

A structural variable (i.e., the column name refers to the structural variable) may have its 
activity level bounded by the user: 

z i £X i £Y i 
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where: Xi represents a particular structural variable. 

represents an upper bound on that structural variable. 
Z-l represents a lower bound on that structural variable. 

The 1620-1311 LP System takes any explicit lower bound and adds it to the value of the struc- 
tural variable and to its upper bound (if any): 

if Z < X < Y 
111 

then < X 1 < (Y + Z ) where X f = X + Z 



and the amount Z\ times the matrix element data value is added to the right-hand side for 
each constraint in which that structural variable occurs . 



Both lower and upper bounds may be negative or positive, but the lower bound must, of 
course, be algebraically less than the upper bound (e.g. , Z^ < Y^); however, the system 
does not check for this. 



Only one matrix element card for a given column name should contain bounds. However, the 
bounds may appear on any one card for a particular structural variable. The first matrix 
element card for a column which contains a bound is assumed to correctly define both bounds 
for that structural variable. If no bounds are used, the variable is treated as having a lower 
bound of zero and no defined upper bound. 

If just the upper bound is used, it must be a positive value since the lower bound is assumed 
to be zero and the upper bound must be greater than the lower bound. If a negative upper 
bound is used, a negative lower bound is also required. 

If just the lower bound is used, the structural variable is treated as having no defined upper 
bound. The lower bound may be either positive or negative. 

If both bounds are used, both may be negative, both positive, or the lower negative and the 
upper positive. Extreme care must be taken in establishing upper and lower bounds to ensure 
that the signs are properly set. Normally, negative bounds are used only when there are both 
sending and receiving types of activities (receipts and disbursements, loans and payments, 
deposits and withdrawals). In these cases the user should set up the convention early as to 
which is the positive activity. There should be a careful check each time for consistency. 

For the objective function data values the user must remember that the practice has been to 
give costs a plus sign and profits a minus sign (just the opposite of what might seem natural) . 
The reasons for this practice are historical: the original simplex algorithms were aimed at 
minimizing the value of the objective function, that is, to make the objective function less 
positive or more negative. So if costs are plus values and profits minus values, minimiza- 
tion (the agendum, MIN. ..) would operate in the proper direction.. The user may, however, 
reverse this usage if he is careful to use the MAX. . . agendum instead of the MIN. . . agen- 
dum. In this case the user may make profits positive values and costs negative values. 

There may be multiple objective functions so that the same matrix may be evaluated against 
different criteria. While the discussion has represented the objective function as containing 
costs, it may use any consistent measure of value, e.g. , time used or quantity produced. In 
the use of MIN. . . and MAX. . . agenda a particular objective function is selected as the 
criteria for optimization. 

In revising an existing matrix the values of the previous upper and lower bounds may be 
changed as well as the previous data value. The 1620-1311 LP System assumes that the only 



42 



values to be changed are those that are explicitly stated on the matrix element card following 
the REVISE agendum card. For instance, if there were originally both an upper and lower 
bound and the revised matrix element card only contained an upper bound value, only the 
upper bound would be changed. If the user wishes to delete a previous bound, he should place 
an asterisk in the first position of the corresponding field on the revised matrix element 
card; e.g., an asterisk in column 31 would delete a previously stored upper bound. 

FIRST. B and NEXT.B INDICATORS 

Format 



cc 



Indicator 


Right- Hand- 


Comments 


Name 


Side Name 




1 7 


8 12 


13 80 


FIRST. B 






FIRST. B 


hhhhh 




NEXT.B 


hhhhh 





Indicator Name: FIRST. B represents the first (or only) right-hand side. The period 

in column 6 is optional. 

NEXT.B represents any additional right-hand side. 

Right-Hand-Side Name: A five-position alphameric field; blanks are significant. Record 
mark and groupmark cannot be used. The first (or only) right-hand side does not need 
to have a name; it is given a name of bbbbb (all blanks) if not otherwise specified. 



Comments: 



Any desired information 



Indicator 
Name 


Input Data 
Source 


Starting Sector 
Location 


Problem 
Name 


Comments 


1 6 


7 


8 12 


13 18 


19 80 


FIRST. 


D 


sssss 


nnnnnn 





cc 



Indicator Name: 



FIRST, represents the first (or only) right-hand side. 



Input Data Source: D indicates that the right-hand- side data values are to be referenced 

from a previous problem. 

Starting Sector Location: A five-position numeric field indicating the assigned location for the 
referenced right-hand- side data values. 

Problem Name: A six-position alphameric field stating the name of the referenced 

problem. This is checked against the problem name stored at the referenced location 
for accuracy. 



Comments: 



Any desired information 
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Usage 



The FIRST. B indicator card is used to introduce the first (or only) right-hand side values. 
The first right-hand side may be named bbbbb (all blanks). The NEXT.B indicator card 
introduces alternate right-hand- side values. Only one right-hand side will be used in a 
particular solution. Each additional right-hand side must have a name. 

There must be a FIRST. B or FIRST. D indicator card if the problem is to be solved. If 
there is a FIRST. B card, all of the data value cards for that right-hand side must follow 
behind it. If there are NEXT.B cards, each one must be followed by its appropriate data 
value cards. 

The FIRST. D indicator card references the data values stored with a previous problem. If 
the previous problem is overwritten, any problem which references this data will not be 
usable. With FIRST. D, all right-hand sides for the previous problem are referenced so that 
any of the previous right-hand sides may be used in the MIN. . . or MAX. . . agenda. 

The MIN. . . and MAX. . . agenda specify the particular right-hand side to be used in deter- 
mining feasibility of the solution. 

RIGHT-HAND-SIDE ELEMENT 

Format 



Blank 


Row 
Name 


Constraint 
Value 


Lower Range 
Limit 


Comments 


1 12 


13 18 


19 30 


31 42 


43 80 


b b 
b b 


rrrrrr 
rrrrrr 


±UUUU.UUUUUU 

ztaiuuu.uuuuuu 


iWW.VVWW 





Blank: Must be all blanks 

Row Name: A six-position alphameric field denoting the exact row name as 

specified in the row identification cards. 

Constraint Value: A twelve -position field permitting a ten-position numeric value with 
space for a sign (optional for positive number) and a decimal point. The format shown 
is recommended. Rules on construction are the same as those under "Matrix Element", 
data value. This field should be used to enter the constraint value if there is only one 
(whether it be an equality, a lower limit, or an upper limit) . It must be used for an 
upper range limit; An asterisk in the first position will delete any previous entry for 
that field. 

Lower Range Limit: A twelve-position field permitting a ten-position numeric value with 
space for a sign (optional for positive numbers) and a decimal point. Construction rules 
are the same as for Constraint Value. This field must be used for a lower range limit. 
It may be used (though not recommended) for any individual constraint value. An asterisk 
in the first position will delete any previous entry for that field. 

Comments: Any desired information 
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Usage 



A right-hand-side element card is required for every row where the value of the constraint 
is not equal to zero. There are five types of constraints which correspond to the row types 
specified on the row identification card: 

Type of Constraint Row Type 

objective function b (blank) 

equality (zero) 

less than + (12 punch) 

greater than - (11 punch) 

range + (12 punch) or - (11 punch) 

For objective functions it is not necessary to have a right-hand-side element card unless it is 
desired to insert a constant term in the value determination. The objective function row has 
the following construction: 

C.X., + C X ... C X + K = Z 
11 2 2 n n 

Assume that the costs (Ci, C£, etc.) are positive values for costs and negative values for 
profits. K represents a constant term. If it is a cost, it would be positive, and negative if it 
is a profit. But only coefficients for structural variables may be stated on matrix element 
cards, so the constant term cannot be covered. The constant term is therefore subtracted 
from both sides of the objective function equation, yielding: 

C^, + C X +... C X =Z-K 
11 2 2 n n 

The K value would therefore be shown on the right-hand-side element card for the objective 
function row, but the sign would be the reverse of a coefficient on the matrix element card. 
For example, if costs are + and profits -, the constant cost would be - and a constant 
profit + and MIN. . . would be used to solve the problem; but if costs are - and profits are +, 
a constant cost would be + and a constant profit - and MAX. . . would be used to solve the 
problem. 

For an equality, less-than, or greater-than constraint, the value should be inserted in the 
constraint value field. 

For a range constraint the upper limit should be placed in the Constraint Value field and the 
lower limit in the Lower Range Limit field. If the type of row (on the row identification 
card) indicates either a less-than or a greater-than constraint, it may be treated as a range 
constraint simply by inserting both an upper and a lower limit. If the lower and upper limits 
are equal, the constraint operates as an equality but the solution process is less efficient than 
if the row type were indicated as an equality. Care must be taken to ensure that the lower 
limit is algebraically smaller than the upper limit; however, the system does not check for 
this. By proper choice of values, a less-than constraint can be changed effectively into a 
greater-than constraint without changing the row identification card. For a row classified as 
+ (less-than constraint) with a value of 3. in the constraint value field, by inserting a large 
negative number in the lower range limit field, the 3. will be interpreted as an upper bound 
with effectively an infinite lower bound. 

The right-hand- side element cards are placed directly behind the appropriate FIRST. B or 
NEXT.B indicator card. Row order is not significant. Right-hand sides can be added during 
INPUT, even when the row identification and matrix files are referenced from disk storage. 
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It is also possible to insert additional right-hand sides during the original input and then 
REVISE these values as desired for additional problem solutions. To modify a right-hand- 
side element value, a new right- hand-side element card is needed. It will change the value 
only of those fields explicitly stated on the new card. If there were previously both an upper 
and lower constraint limit and the new card contained only an upper limit, just the upper limit 
would be modified. If the user wishes to delete a value he inserts an asterisk in the first 
position of the corresponding field on the new card. 

BASIS. INDICATOR 

Format 



Indicator 
Name 


Input Data 
Source 


Starting Sector 
Location 


Problem 
Name 


Comments 


1 6 


7 


8 12 


13 18 


19 80 


BASIS. 
BASIS. 


D 


sssss 


nnnnnn 





Indicator Name: BASIS. The period in column 6 is optional. 

Input Data Source: b (blank) indicates that basis cards will follow. 

D indicates that the basis will be referenced from a previous problem. 

Starting Sector Location: A five-position numeric field indicating the disk storage location 
of the problem containing the desired basis. 

Problem Name: A six-position alphameric field containing the problem name which is 
checked against the problem name stored at the referenced storage location. 

Blank: Must be all blanks 

Comments: Any desired information 

Usage 

Basis cards are used to provide an advanced start in solving a linear programming problem. 
The basis is a statement of the particular structural variables which are in the initial solution, 
plus information as to which logical variables are to be omitted from the initial solution. 
Normally, basis cards are available from previous LP solutions either of the entire problem 
or of portions of the total problem. However, with care a user can insert an initial basis 
directly. 

Basis cards for recycling come from OUTPUT and GETOFF, two of the agenda in 1620-1311 
LP. The BASIS, indicator card must precede both the VARBLS indicator card, which 
introduces the structural variables to be placed in the initial solution, and the SLACKS 
indicator card, which introduces the logical variables to be omitted from the basis. 

The BASIS, indicator card provided by OUTPUT contains information on mantissa length 
and tolerances which is used by the program described in "1620-1311 LP Agenda" under 
"Output". 
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If the user is preparing the BASIS, indicator card manually, he may provide this information 
by means of an ASSIGN agendum card which is covered later under "ASSIGN". Experience 
with certain problems will indicate which values may lead to efficient solution. 

The BASIS. D card references a basis established in a previous problem. No other basis 
cards may follow a BASIS. D card. 

VARBLS INDICATOR 

Format 



Indicator Name 


Comments 


1 6 


7 80 


VARBLS 





Indicator Name: VARBLS 

Comments: Any desired information 

Usage 

The VARBLS indicator card precedes the basic structural variable cards. Only one VARBLS 
indicator is permitted per BASIS, indicator card. The VARBLS indicator card provided by 
OUTPUT has column headings in the comments field (e.g. , TYPE, NAME, ACTIVITY 
LEVEL). 

BASIC STRUCTURAL VARIABLES IDENTIFICATION 
Format 



Blank 


Basis Entry 
Type 


Column 
Name 


Comments 


1 11 


12 


13 18 


19 80 


b b 


e 


cccccc 





Blank: Must be all blanks 

Basis Entry Type: There are three levels at which a structural variable may be entered 
into a basis: 

W - The structural variable is entered at its lower bound; these are not 
required for basic structural variables. Remember that there is 
always an implied lower bound of zero unless otherwise specified. 

G - The structural variable is entered at its upper bound. 

F - The structural variable is entered at an intermediate level; an un- 
bounded variable could be noted only as type F. 
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Column Name: A six-position alphameric field stating the name of a structural variable 
which is to be in the initial basis. 

Comments: Any desired information 

Usage 

If the basis cards come from OUTPUT, they will also contain the solution activity level in 
columns 20-33. W-type cards (structural variables in at their lower bound) will also be 
included unless there is an implicit lower bound of zero. These W cards are ignored in 
establishing the initial basis. 

If the user is trying to manually provide an advanced starting basis, he should select those 
structural variables which he feels quite certain will be present in any reasonable solution 
to the problem. A selected structural variable might be a low-cost ingredient, or an item 
which is the only source for a certain component of the final product, or a limiting resource 
(a particular machine or process) . If an improperly formed basis is entered, a poor starting 
solution may result. It is therefore suggested that the new LP user avoid manually preparing 
a basis until he has gained considerable experience in the characteristics of the problem 
with which he will be dealing. 

The activity level, other than upper (G) or intermediate (F), is not needed since it is 
calculated by 1620-1311 LP System automatically after the problem is entered. If basis 
cards have been produced by more than one previous problem, they can be combined by 
inserting all the basic structural variables cards behind one VARBLS indicator card, and 
nonbasic logical variables cards behind the SLACKS indicator card. 

SLACKS INDICATOR 

Format 



Indicator Name 


Comments 


1 6 


7 80 


SLACKS 





Indicator Name: SLACKS 

Comments: Any desired information 

Usage 

The SLACKS indicator card precedes the nonbasic logical variable cards. Only one 
SLACKS indicator is permitted per BASIS, indicator card; it must come after the basic 
structural variables cards. The SLACKS indicator card provided by OUTPUT has column 
headings in the Comments field (for example, TYPE, NAME, ACTIVITY LEVEL, 
SIMPLEX MULT.). 
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NONBASIC LOGICAL VARIABLES IDENTIFICATION 
Format 



Blank 


Basis Entry 
Type 


Row 
Name 


Comments 


1 11 


12 


13 18 


19 80 


b b 


e 


rrrrrr 





Blank: Must be all blanks 

Basis Entry Type: W - A logical variable at its lower bound. 

G - A logical variable at its upper bound. 

F - A logical variable at an intermediate level. 

Row Name: A six-position alphameric field containing the row name corresponding to the 
logical variable which is to be removed from the basis. 

Comment: Any desired information 

Usage 

For every structural variable inserted in the basis, a corresponding logical variable must 
be removed from the basis. Removing the logical variable is accomplished by naming the 
constraint row in which that logical variable occurs. The row named should be one in which 
the basic structural variable has a nonzero coefficient. 

The OUTPUT and GETOFF agenda produce nonbasic logical variable cards which can be 
recycled. These cards also contain certain values (activity level and simplex multiplier); 
these values will be ignored when the cards are recycled. All F-type cards will be ignored. 

ENDATA OR EOF INDICATOR 

Format 



Indicator Name 


Comments 


1 6 


7 80 


ENDATA 
EOF 





Indicator Name: ENDATA or EOF. 
Comments: Any desired information 

Usage 

ENDATA or EOF is used as a terminal indicator. It must be placed at the end of a set of 
data input cards. Before any agendum card is used, following the data input cards, either 
an ENDATA or EOF indicator card is required. 
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1620-1311 LP AGENDA 



This section of the manual gives the details of each of the agendum statements. Thus, it is 
more of a reference than a tutorial section. It contains four major categories, each with an 
introduction, followed by the agendum statements: 

1. Data Preparation 

INPUT. 

REVISE 

2. Optimization 

MAX... 
MIN. . . 

INVERT 

3. Report Preparation 

OUTPUT 
DO.D/J 
COST.R 
CHECK. 

4. Control and Data Maintenance 

ASSIGN 

ENDJOB 

ERASEp 

GETOFF 

MAP... 

PAUSE . 

SAVE.p 

DATA PREPARATION 

It is necessary to provide data in order to initiate problem solution. There are two agendum 
statements to accomplish this: INPUT, and REVISE. INPUT, brings new data into disk 
storage or references data already in storage. REVISE changes information in disk storage. 
These agenda have the data cards following them. The data input cards have been explained 
in detail in "Data Preparation". 

INPUT . 

Format 



Agendum 


Data 


Starting Sector 


Problem 


Comments 


Name 


Source 


Address 


Name 




1 6 


7 


8 12 


13 18 


19 80 


INPUT. 


C 


sssss 


nnnnnn 




INPUT. 


D 


sssss 


nnnnnn 
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Agendum Name: INPUT. 



Data Source: C indicates that input data cards will follow. 

D indicates that the input should be taken from a previous problem 

stored on disk. 

Starting Sector Address: A five-position numeric field identifying the initial sector location 
for storing (or retrieving) the problem data. The address must be 01805 or greater. 

Problem Name: A six-position alphameric field; any symbols may be used except record 
mark and groupmark. Blanks are permitted but not recommended as internal charac- 
ters of a name. For disk input the problem name must correspond exactly with the 
original problem name at the referenced location. 

Comments: Any desired information 

Examples 

INPUT. C03000BLENDG 
INPUT. D12000ALLOC 

Description 

Portions of disk module are used for table, program, and monitor residence, and are 
not available for storage of problem data. Sectors 01805 through 04799 are available. If 
additional room is necessary, the user must specify, on an ASSIGN card, which sectors 
beyond the programs and below the monitor should be made available. A description of the 
ASSIGN card is given later in this section, and a list of the unavailable areas on module 
appears immediately under "Disk Storage Utilization". 

The INPUT, agendum is used to set up the linear programming model for solution or to set 
up data for special LP analysis or output. This agendum will establish the input files in 
disk storage from card input, or will prepare to work with the referenced input files in 
disk storage. INPUT, will reference a previously stored problem, or will input one or more 
card files. These files are: 

• Row identification file, preceded by ROW. ID indicator card 

• Column identification file, preceded by COL. ID indicator card. This file is optional 
and is used primarily to select a subset of variables from the problem to be used during 
optimization or for selective output. 

• Matrix data file, preceded by a MATRIX indicator card. This file defines the data 
value (coefficients) and bounds of the structural variables . If no column identification 
file is given, a column identification file is generated from the matrix file. 

• Right-hand-side file is preceded by a FIRST. B indicator card. This establishes the 
data value (coefficients) of the principal or only right-hand side. Additional right-hand 
sides can be defined; each is preceded by a NEXT. B indicator card. 

• Basis identification files. These files are preceded by a BASIS, indicator card. These 
files are optional. They are used to establish an advanced starting basis. The struc- 
tural variable file is preceded by a VARBLS indicator card; the logical variable file is 
preceded by a SLACKS indicator card. This file was probably generated by a previous 
OUTPUT agendum. 
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• Input data is terminated by an EOF or ENDATA indicator card. 

An input card which references disk storage (INPUT. D) completely defines the input; no 
indicator card or card file can follow this card. 

If any input files are on cards, an INPUT. C card must be used. Data may also be inserted 
through special disk cards, e.g., ROW.IDD, COL.IDD, MATRIXD, FIRST. D, and BASIS. D. 
All the problem data may be read in from cards, but it is usually more efficient to define 
a problem wholly or in part by reference to data stored on disk. The data file for the var- 
iables need be read only once. Selective column and row identification files then serve to 
select the data as required for each particular problem. Right-hand sides and objective 
functions can also be varied. By this approach, data maintenance is simplified since one 
master data value file is maintained. 

The following illustration shows a complete set of card input identification and data files. 
Unless noted otherwise, all cards are required to set up a new problem. If the indicator 
card is present, the file cards following the indicator card are required unless the indicator 
card refers to information in disk storage. Each of the cards is described in detail in 
"Data Preparation". 

REVISE 

Format 



Agendum Name 


Comments 


1 6 


7 80 


REVISE 





Agendum Name: REVISE 

Comments: Any desired information 

Example 

REVISE 
Description 

The REVISE agendum permits the modification of existing problem data. REVISE works 
only on the current problem, so it must be preceded by an INPUT, card, though other cards 
may intervene. For example, INPUT. C, MIN. . ., OUTPUT and then REVISE, is an 
acceptable sequence. 

REVISE may be followed by any of the input indicator cards, in turn followed by the desired 
revision data and identification cards: ROW. ID, row identification, COL. ID, column 
identification, MATRIX, matrix elements, FIRST. B, right-hand-side elements , BASIS., 
VARBLS, basic structural variables, SLACKS, nonbasic logical variables, EOF or 
ENDATA. 

The order and use of the data input cards is the same as with INPUT. , except that no disk 
references may be made (for example, no ROW.IDD or MATRIXD entries). 
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next agendum card 



ENDATA or EOF 



■> (optional) 



nonbasic logical 
variables 



SLACKS 



-► (optional) 



(optional) 



starting basis structural 
variables 



VARBLS 



BASIS. 



Input 

Indicator 

Cards 



-> (optional) 



additional RHS 
element file 



NEXT. B 



J 



1st RHS element file 



FIRST. B 



A 



matrix element file 
(coefficients and bounds) 



MATRIX 



A 



(optional) 



column identification file 



COL. ID 



Agendum 
Card 



row identification file 



ROW. ID 



INPUT. C sssss nnnnnn 



repeated for 
each RHS 
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Since REVISE works on the current problem, it must be preceded at some previous point by 
INPUT. REVISE actually changes the identification information and data stored for a particu- 
lar problem, so if any other problem uses that data, the second problem will in turn be 
affected by the data. 

REVISE can only change information already in disk storage. It will alter any entry in any 
of the files of the current problem, but it cannot add any new row or column or elements 
unless provision has been previously made for the entry. 

• Row identification file. — Any row type can be revised. No rows can be added. Rows 
can be omitted by using an asterisk (*) on a row identification card. The type of gener- 
ated logical variable can be changed by varying the row type . Rows can be inserted by 
changing an asterisk row type to any other row type. 

• Column identification file . — Any structural variable can be revised and inserted in the 
matrix by putting in a column identification card with a blank for column usage . A 
column may be omitted by putting an asterisk in for column usage. No column names 
can be added. 

• Matrix element file. — Any entered data value or bound can be revised. It is not pos- 
sible to revise an implicit zero data value or to specify a bound for a previously 
unbounded variable. All values for a given matrix element card must be included in 
the revised data input. 

• Right-hand-side element file. — Any entered constraint value can be revised. It is also 
possible to change a simple inequality to a range constraint. The first indicator card 
must be a FIRST. B indicator card. The right-hand- side name can be the name of any 
right-hand side. 

• Basis identification file. — (1) If basis was saved by SAVE. B or SAVE. I, basis entry 
type (F and G only) of structural variables can be revised. Similarly, basis entry type 
(W and G only) of logical variables can be revised. (2) If basis was inserted as part of 
an input deck, any of the entries may be revised. 

• Data revision must be terminated by an EOF or ENDATA indicator card. 

Based on information from the various output-producing agenda, the user may want to revise 
various problem elements. Or, as the result of changing conditions, various factors may 
have to be changed. 

The user may omit rows or columns that are not highly significant to the solution. This is 
done by inserting a card with the particular row or column name with an asterisk for row 
type or column usage . The bounds on a structural variable may need to be changed to reflect 
a different inventory position or technological requirement. The matrix element values may 
have to be modified to reflect changing use of resources to accomplish different activities. 
The right-hand- side element values may change to recognize different capacities or demands. 
The objective function values may need to be changed to consider different prices or resource 
costs. 

For solution efficiency it may be desirable to modify a previous basis to compensate for 
changes in the matrix elements, bounds, or constraints. 

REVISE provides a method for changing right-hand-side values and objective functions. This 
multiple examination can also be made through the use of alternate right-hand-sides and 
objective functions in the original problem input. 
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Revised problems can usually be solved rapidly by using the basis from a previous solution 
of the problem. The basis can be retained for later use by the SAVE. B agendum. Then this 
basis will be automatically used as the starting basis for re-solution of the problem. If the 
user does not want to start from a saved basis, he must use the ERASEB agendum before 
beginning the solution. 

OPTIMIZATION 

The optimization phase of a linear programming problem takes the data provided by the input 
programs or the revise programs and determines the solution which has the maximum profit 
or the minimum cost, whichever is specified by the user. This phase of processing takes 
considerably longer than the input or output phases except for very small problems. The 
number of rows and the number of nonzero elements will have a significant effect upon the 
time required to solve a problem, as does the variation in the size of the numbers. It is 
extremely difficult to estimate the running time for a problem, even when the number of 
rows and variables and the density are known. 

Unless a starting basis is provided by the user, the 1620-1311 Linear Programming System 
will start from an all-slack basis which it generates. Providing a good basis from a 
previously solved variation of the problem will usually dramatically reduce the solution 
time. 

The 1620-1311 LP System uses the Dual algorithm to produce a feasible solution. The 
Primal algorithm then automatically takes over to obtain an optimal feasible solution. A 
sophisticated matrix inversion routine is also used. Inversion is called for automatically 
every 15 iterations unless the user chooses to change the inversion frequency with an 
ASSIGN card. In addition, INVERT will be called when the Dual or Primal algorithms 
consider it necessary to maintain accuracy. 

The optimization procedure, called by a MIN. . . or MAX. . . agendum card, will automatically 
calculate and assign a mantissa length according to the problem it is to solve. During the 
course of the calculation, when the CHECK, routine is executed before each inversion, the 
mantissa length will automatically be increased if necessary to maintain the accuracy 
required. The user may override the mantissa length the system selects, by using an 
ASSIGN card. The user may also stipulate the accuracy he requires by assigning a maxi- 
mum error tolerance (see "ASSIGN" page 2. 15.39). All these features are automatic within 
the LP system, unless the user chooses to override them because of special knowledge 
about the problem or special requirements. 

In addition to reducing the size of the numbers used in the calculation by using an appropriate 
mantissa length to fit the problem, a further speed increase is made by calculating, when 
possible, with the normally smaller original input numbers, maybe two or three digits, 
even when the mantissa length may be eight or ten. This gain is made possible by the nature 
of the machine and use of the revised simplex method, product form of the inverse. Addi- 
tionally, the inversion method is designed to maintain sparsity, thus reducing the number 
of calculations. See "Mathematical Description of 1620-1311 LP" for a detailed mathematical 
description. An additional feature, the bounded variable algorithm, significantly improves 
performance for problems which have bounded variables and/or range constraints. 



55 



MAX..., MIN... 



Format 



Agendum 
Name 


Spacer 


Right- Hand- 
Side Name 


Objective Function 
Name 


Comments 


1 6 


7 


8 12 


13 18 


19 80 


MAX... 
MAX... 
MIN. . . 
MIN. . . 


• 
• 


hhhhh 
hhhhh 


ffffff 
ffffff 





Agendum Name: MAX. . . or MIN. . . 

Spacer: An unused position. A period is suggested to ensure column alignment. 

Right-Hand-Side Name: The name of the right-hand side to be used in the optimization. 
The right-hand-side name has a maximum of five characters. If the field is left blank, 
the first (or only) right-hand side is used. 

Objective Function Name: The name of the objective function row to be used in the optimi- 
zation. Whenever the field is left blank, the first (or only) objective function is used. 

Examples 

MAX... 

MAX ALOY2 

MIN BEEFbCOST2 

MAX bbbbbPROFIT 

Description 

MAX. . . or MIN. . . (maximize or minimize) calls the main solution procedures to compute 
the optimal value of the selected objective function. MAX. . . is generally used with profits 
expressed as positive numbers and costs expressed as negative numbers in the objective 
function; the user wishes to maximize profit. MIN. . . is used to minimize cost and the signs 
of the objective function would be reversed. 

The 1620-1311 Linear Programming System permits the user to include multiple 
right-hand sides and objective functions. For any optimization, the user selects the 
right-hand side and objective function he desires by name. Note that right-hand 
sides may have only five-character names. If the user does not specify a right-hand 
side, the first one will be used; the same is true for the objective function row name. 
Neither, either, or both may be specified. 
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INVERT 



Format 



Agendum 
Name 


Spacer 


Right-Hand- 
Side Name 


Objective Function 
Name 


Comments 


1 6 


7 


8 12 


13 18 


19 80 


INVERT 




hhhhh 


ffffff 





Agendum Name: INVERT 

Spacer: An unused position. A period is suggested to ensure column alignment. 

Right-Hand-Side Name: The name of the right-hand side to be used in the inversion. 

The right-hand- side name has a maximum of five characters. If the field is left blank, 
the first (or only) right-hand side is used. 

Objective Function Name: The name of the objective function row to be used in the inversion. 
If the field is left blank, the first (or only) objective function row is used. 

Examples 

INVERT 

INVERT. ALOY2 
INVERT. ALOY2COST2 
INVERT, bbbbb PROFIT 

Description 

The INVERT agendum is used to solve a set of simultaneous equations. It is never neces- 
sary for the 1620-1311 LP user to place an INVERT card in a linear programming agenda. 
If it is used in a linear program agenda, that is, preceding a MIN. . . or MAX. . . agendum 
card, it will force an initial pass through the invert routine. This pass is then bypassed 
when the MIN. . . or MAX. . . card is read. If no INVERT agendum precedes the MIN. . . or 
MAX. . . card, the LP system will take an initial pass through the invert routine. 

The purpose of the INVERT agendum is: 

1. To solve simultaneous equations 

2. To be compatible with other linear programming systems which require the INVERT 
agendum to precede their optimization agendum when a basis is given 

The solution to a set of simultaneous equations requires: 

1. An INPUT, agendum and input data 

a. Row identification file, including an objective function row 

b. Matrix 
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c. Right-hand side(s) 

d. Basis* 

2. An INVERT agendum card 

3. An OUTPUT agendum card 

4. (Optional) A CHECK, agendum to test inversion accuracy 

DO. D/J and COST. R should not be used (they are meaningless in this case). 
REPORT PREPARATION 

Various reports are available from the IBM 1620-1311 Linear Programming System. These 
reports contain, among other things, the objective function value, activity levels, simplex 
multipliers, reduced costs (Dj's), marginal values, row errors, submatrix summaries, 
and cost ranges. The particular reports are called by agendum cards (OUTPUT, DO. D/J, 
COST. R, CHECK.). Any or all of these reports can be obtained after an optimal solution 
is reached. They may be obtained in any order desired by the user. In addition, an iteration 
log is maintained to show how the solution is proceeding. 

Before describing the reports and how they are obtained, a few basic nonmathematical 
definitions are presented. For fuller definitions, see the glossary of the IBM manual "An 
Introduction to Linear Programming" (E20-8171). 

Activity level — the number of units of a structural variable or a slack in the current 
solution. 

Objective function value — the total cost (or profit) of the current solution. The alge- 
braic sum of the unit costs and/or profits times their activity level yields the 
objective function value. 

Simplex multiplier — the objective function row of the basis inverse. This row, except 
for signs, is the set of marginal values. The simplex multipliers are used 
internally to calculate the reduced costs. 

Marginal value — the amount the objective function value would increase if one more 
unit of the desirable right- hand-side element were allowable; or, for undesirable 
right-hand-side restrictions, the amount the objective function value would increase 
if one less unit were required. 

Reduced cost (Dj) — the objective function coefficient change required to alter the 
solution activity level of a variable. (1) For a variable without bounds, it is the 
reduction in cost (or increase in profit) required to tie for a position in the solution. 
A change in objective function coefficient exceeding the reduced cost will force the 
variable into the optimum basis. (2) For variables in the solution at a specified 
lower bound, the same definition applies relative to an increase in solution activity 
level. The reduced cost of a variable at lower bound can also be interpreted as 
the unit cost of forcing an additional unit into the solution. (3) For a variable at its 
upper bound, it is the unit profit of allowing an additional unit into the solution. 



* If the basis does not include all of the variables, see the description of the INVERT agen- 
dum in "Mathematical Description of 1620-1311 LP" for an explanation of the results. 
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Row error — the computational error associated with a row, caused by accumulated 
roundoff. The row error is calculated by summing the products of each row coeffi- 
cient by its activity level and subtracting from the right-hand side element. 
Excessive row errors, when not handled automatically by the system, can generally 
be handled by increasing the floating-point precision. 

Submatrix summary — the product of the solution level of selected variables with row 

coefficients summed for the row. One or more rows may be selected. For example, 
in a feed blending problem this could be used to determine the protein and fat due 
to vegetable ingredients. 

Cost range — the upper and lower limits which an objective function element may take 
on and remain in the optimal basis. If the cost (or profit) coefficient changes to a 
point outside the range, a basis change will occur resulting in an increase (or 
decrease) in solution value of the variable being analyzed. 

There are two report modes: (1) where full reports are desired, and (2) where reports of 
only selected portions are desired. Full or partial results can be obtained from each of the 
reports. While the details of each of the reports follow, a quick summary is given to show 
what values are available and how to get them. 

The full reports are obtained by simply placing the agendum card(s) in the agenda deck 
following the MIN. . . or MAX. . . cards (other cards may intervene). The four primary 
reports are called for as follows and contain the significant values indicated. 



1 6 
OUTPUT 



The OUTPUT report contains as the first card (or a line of print on the typewriter Or 
printer) a BASIS, indicator card. This deck can be used as an advanced basis in solving 
subsequent related problems when placed in the input deck, since the format is compatible. 
For the basic structural variables, the activity levels are given. For the logical variables 
(slacks and artificials), the activity levels are given for basic variables, and the simplex 
multipliers for the nonbasic variables. For details see "Output", later in this section. 



16 
DO.D/J 

The DO.D/J report contains the reduced costs (Dj's) for nonbasic variables and basic 
variables which are forced into the basis because of upper or lower bounds. This report 
also gives the current cost and the basis value (difference between current cost and reduced 
cost). Each of the rows in the problem is given the type of row (+, -, =, RANGE), the row 
name, and the marginal value. For details see "DO. D/J" later in this section. 



1 6 
COST.R 
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The COST.R report contains the cost (or profit) interval over which the objective function 
coefficient may range and remain in the optimal basis. In addition to the upper and lower 
limit of the range, the report will give the variable that will enter if the upper bound were 
exceeded (by a cost change), and the change in basis that would occur if the cost were to go 
below the lower bound. For details see "COST.R" later in this section. 



1 6 
CHECK. 



The CHECK, report will give the row errors and the maximum error (the row error of the 
row having the largest error); in addition, the CHECK, report includes the specified row 
limit or limits and the solution row activity level. For details see "CHECK." later in this 
section. 

The illustration below shows how a user might set up his deck when he desires the full 
reports. The order of the report agendum cards is not significant. 



ENDJOB 
DO. D/J 



COST.R 



CHECK. 



OUTPUT 



MIN. . . 




Input deck 



INPUT. C 



+ +XEQSLP1620 

The format (and method of obtaining selective output from the four output reports) is similar 
and is briefly summarized. A more detailed explanation is given with each agendum 
description. 



Agendum 
Name 


Data 
Source 


Starting Sector 
Location 


Problem 
Name 


Decimal Digits 


Comments 


Field 1 


Field 2 


1 6 


7 


8 12 


13 18 


19 


20 


21 80 


OUTPUT 


D 


sssss 


nnnnnn 


dl 


d 2 




CHECK. 


D 


sssss 


nnnnnn 


dl 


d 2 




DO. D/J 


D 


sssss 


nnnnnn 


dl 






COST.R 


D 


sssss 


nnnnnn 


dl 
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Agendum Name: 



The agendum call card OUTPUT, CHECK., DO.D/J, or COST.R. 



Data Source: D indicates that the referenced identification files have been previously 

placed on the disk by INPUT. . 

Starting Sector Location: This five-digit field is used to give the sector address of the row 
and/or row column identification files containing the names of the rows and/or columns 
(a subset of the entire row and/or column identification files for the problem) for which 
output is to be prepared for the report. The address must correspond to the address 
given on the INPUT. C card for the desired file. 

Problem Name: Used to specify the name of the subproblem to be reported on. The 
name must correspond to the name given on the INPUT. C card for the desired file. If 
the name does not take the full six characters, the names should line up on a character- 
by-character basis; for example, PROD lb is not the same as bPRODl. 

Decimal Digits: The number of desired decimal digits to follow the decimal point for 
the fixed-point output. 

Field 1 blank: Three decimal digits are given as standard for activity level, reduced 

cost, row solution level, and all costs for OUTPUT, DO.D/J, CHECK., and COST.R 
respectively. A digit in this position specifies the number of decimal digits desired. 

Field 2 specifies the number of decimal digits for row error and simplex multipliers in 
CHECK, and OUTPUT respectively. Standard for row error and simplex multipliers 
is 8 and 3, which is given if the field is left blank. Does not apply to DO.D/J and 
CHECK. 

The various output reports available to the 1620-1311 linear program user may be obtained 
on the typewriter, on cards, or on an online printer. The standard output device is the 
card punch. If another device is desired, the user may assign it using the ASSIGN agendum 
card (see "ASSIGN"). For example, to cause the output results to be printed on an online 
printer, the user would enter the following agendum card: 



cc 



OUTPUT 



ASSIGN 9 
16 16 



To put the results on the typewriter the user would precede any output agendum cards with 
an assign card, as follows: 



DO.D/J 

ASSIGN 



cc 1 6 16 
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To reassign to the card punch would call for an ASSIGN card with a 4 punch in column 16. 



COST. R 



ASSIGN 



1 (typewriter) 



DO.D/J 



ASSIGN 



4 (cards) 



OUTPUT 



ASSIGN 



9 (printer) 



MAX... 



INPUT. 



The agenda above would put the normal OUTPUT on the printer, the Dj T s (reduced costs) on 
cards, and the cost ranges on the typewriter. Note that any INPUT, card will reinitialize the 
output unit to card. For assigning output devices for the iteration log and the agenda log 
see "ASSIGN". 

The Iteration Log, unlike the other output reports, is not controlled by agendum cards. At 
the completion of each simplex iteration a one-line item is prepared to indicate what hap- 
pened on that iteration. Through the use of the ASSIGN card, the user can obtain this 
information directly on the typewriter, on the printer, or punch a card deck. If the user 
chooses not to log every iteration, he can, through the ASSIGN card, assign the frequency 
of logging that he chooses. Additional control of the iteration log is provided through console 
switch 1. By putting console switch 1 ON, every iteration will be logged on the typewriter. 
When it is turned OFF, the normal logging frequency is maintained. The standard frequency 
for logging is every iteration. The standard output device for the iteration log is the 
typewriter. 

The headings and the meaning of the values to which they apply follow: 
ITER: The iteration number. 

VARBL. EXIT: Name of the variable leaving the basis. "U.B." indicates that the 
variable is leaving at its upper bound. 

VARBL. ENTR: Name of the variable entering the basis. "U.B." or "L.B." indi- 
cates that the variable is entering at either its upper or lower bound. 

CNT: While in the dual algorithm, the number indicates the number of infeasibilities 
at the beginning of the iteration. After the solution has become feasible and the 
primal algorithm takes over, the field indicates the number of desirable variables 
on the last iteration. 

OBJ. FUNCTION: The current value of the objective function. 
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OUTPUT 



Format 



Agendum 
Name 


Data 
Source 


Starting Sector 
Location 


Problem 
Name 


Decimal Digits 
Field 1 Field 2 


Comments 


1 6 


7 


8 12 


13 18 


19 20 


21 80 


OUTPUT 
OUTPUT 


D 


sssss 


nnnnnn 


d l d 2 
d l d 2 





Agendum Name: OUTPUT 

Data Source: blank - complete output of all basic variables and all slacks plus 

simplex multipliers for all nonbasic slacks. If Data Source is blank, the Sector 
Location and Problem Name fields will also be blank. 

D-Output will be selected from a file specified on the disk. 

Starting Sector Location: The address specified on the INPUT. C card for the row and/or 
column identification files which provide the names of the desired basic structural 
variables and/or the names of the desired rows for the output report. 

Problem Name: Provides the name of the file specified on the INPUT. C card which 

gives the names of the desired variables and/ or rows for the output report. The problem 
name is used to check that the user has specified the correct Starting Sector Location. 
An error condition will be indicated if there is not an exact match. 

Decimal Digits: If either Field 1 or Field 2 is blank, three decimal digits will follow 
the decimal point. 

Field 1 — The digit specified will indicate the number of desired digits following 
the decimal point for the activity level field in the output report. 

Field 2 — The digit specified will indicate the desired number of decimal digits to 
follow the decimal point in the simplex multiplier field on the output 
report. 



Examples 

OUTPUT 

OUTPUT 25 

OUTPUTD03000SHIPNG 

OUTPUTD02500PRODlbl6 
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Description 



The OUTPUT report contains three primary pieces of information: (1) the activity level for 
basic variables and slacks, (2) the values of the simplex multiplier for the rows containing 
noribasic slacks, and (3) the current value of the objective function. The activity level indi- 
cates the number of units of the variable in the optimal solution. The activity level of a slack 
indicates the difference between the sum of activity levels for a row and the right-hand-side 
value. For example, in a resource availability constraint row, the slack would indicate the 
amount of the resource not used. The simplex multipliers (the row evaluators) have the 
same absolute value as the marginal values and indicate how uneconomic it is to bring the 
corresponding slack variable into the basis. 

In addition to displaying the primary solution to a linear programming problem, the OUTPUT 
agendum also performs another valuable function. If the output is placed on cards, the same 
cards may later be reintroduced as input to provide an advanced starting basis for re-solution, 
Starting with a good advanced basis will substantially reduce the time to resolve the problem 
even when some significant changes have been made to the objective function, right-hand side, 
or matrix. To accomplish this, three input indicator cards are produced by the OUTPUT 
routine. These cards also provide the headings for the OUTPUT report. The indicator cards 
are BASIS. , VARBLS, and SLACKS. 

The output deck can then be taken intact and placed behind the matrix deck when inputting a 
subsequent problem. 

If the output is placed on cards, sequence numbers are punched in columns 76-80. 
Sample Report 



OUTPUT 

BASIS._ MANTISSA 10 

VARBLS TYPE NAME ACTIVITY LEVEL 

FBIN2 190.678_ 

WBIN3 500.000 

_ FBIN4 ^233 .051„ 

FALUM 961 .864 

FSILC0N. _ 1 14.407 
SLACKS TYPE NAME ACTIVITY LEVEL 

FVALUE 373..686 _ 

FVALUE2 373.686- 

. ___0G YIELD .000_ 

+WFE 

+WCU _ 

+FMN 22.712 

+FMG 24.280^ 

-FAL 54.619 

_.. ... -WS I 



TOLERANCES 05 05 05 05 05 



SIMPLEX MU LTV 



.264- 
"2.843" 
2.754 



.208- 



Report Description 



The report resulting from the OUTPUT agendum appears in two sections. Each section of the 
report is preceded by its indicator/heading card. In addition there is a BASIS, indicator 
card which precedes the entire OUTPUT report. It is this indicator card that is recognized 
by the INPUT agendum card when the output deck is being recycled to provide an advanced 
starting basis. In addition, this card indicates the mantissa length and tolerances which were 
set during input processing or during the calculation. (See " FORMAT " for definitions of 
these lengths and tolerances. ) 
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Since the OUTPUT report, when punched on cards, can be recycled as input, significant card 
columns will be indicated. 



"RASTS 


MANTTSSA tcx 


TOLERANCES t t tt tt tt tt xxxxx 
11 22 33 44 55 


rol 1 -fi 


BASIS. 


TnHipatoT'/Tifiadinp' 


rnl 90-97 


MANTISSA 


Title 


col 29-30 


XX 


Mantissa length 


col. 35-44 


TOLERANCES 


Title 


col. 47-48 


h 


The element tolerance 


col. 50-51 


h 


The pivot tolerance 


col. 53-54 




The feasibility tolerance 


col. 56-57 




The objective function tolerance 


col. 59-60 




The row error tolerance 


col. 76-80 


xxxxx 


The output card sequence number 



After the BASIS, indicator card, the variables section of the output report appears preceded 
by its indicator card. Except for the indicator VARBLS in columns 1-6, the card contains 
just column headings. 



VARBLS: Indicates that the structural variables follow. 

TYPE: Indicates the type of structural variable. Type W indicates a variable at its lower 
bound (implicit lower bound of zero excepted). Type G indicates a variable in the basis 
at its upper bound. Type F indicates a variable at an intermediate level. 

NAME: Contains the name of the structural variable. 

ACTIVITY LEVEL: Indicates the number of units of the named variable in the current solu- 
tion. The activity level will be to three decimal digits unless otherwise specified in 
column 19 of the OUTPUT agendum card. 

Card Count: The output cards are numbered sequentially. 

The second section of the report is preceded by a SLACKS indicator in positions 1-6. This 
section contains information about the logical or slack variables. 

SLACKS: Indicates that the slack variables follow. 

TYPE: Indicates the type of slack variable — type F when at an intermediate level, type W 
when at the lower bound, and type G when at the upper bound. Range constraints and 
equations can be type G. 

NAME: Contains the name of the row for the slack (logical) variable. 

ACTIVITY LEVEL: Indicates the number of units of the slack variable in the basis. First 
entry in this section of the report generally contains the value of the objective function. 
There will not be a value in this field for type W slacks. 

SIMPLEX MULT: This field contains the value of the row simplex multiplier. SLACKS at 
an intermediate value, type F, since they are in the basis, will not have a simplex 
multiplier value. The simplex multiplier for a logical (slack) variable indicates how 
uneconomical it would be to bring the slack variable, identified by its row name, into 
the basis. The simplex multipliers for the implicit constraints expressed by the bounds 
are available from DO. D/J. 
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Card Count: The output cards are numbered sequentially. 

The OUTPUT agendum is sometimes used when an optimal solution has not been obtained. 
This might happen when the problem has not been formulated properly or some unexpected 
difficulty has occurred. The output cards can be recycled to restart the solution from an 
advanced basis. 

DO. D/J 

Format 



Agendum 


Data 


Starting Sector 


Problem 


Decimal 


Comments 


Name 


Source 


Location 


Name 


Digits 




1 6 


7 


8 12 


13 18 


19 


20 80 


DO. D/J 








d 




DO. D/J 


D 


sssss 


nnnnnn 


d 





Agendum Name: DO. D/J 

Data Source: Blank. — All Dj's (reduced costs) will be produced. If Data Source is blank 
then Starting Sector Location and Problem Name will also be blank, 
D - Dj's will be selected from a file specified on disk. 

Starting Sector Location: The address specified on the INPUT. C card which contains the 
identifications of the desired variables and/or rows. 

Problem Name: The name of the file specified on the INPUT. C card. The name of the file 
on the disk must match the name on the DO. D/J card or an error condition will exist. 

Decimal Digits: The number of digits the user wishes\following the decimal point for the 
reduced costs. If blank, three digits will be given. 

Examples 

DO. D/J 

DO. D/J 6 
DO. D/JD03000SHTPNG 
DO. D/JD02700DISTRB2 

Description 

Dj's (the reduced costs) indicate the uneconomic change in the value of the objective function 
if a unit of the undesirable activity (variable) is brought into the solution. Or to look at it 
another way, the reduced cost indicates the difference between the current cost of a variable 
and the lower cost at which it would, if in the basis, yield the same optimal objective function 
value. This lower cost is called the basis value in the report. If the cost were reduced still 
more, the variable would enter the optimal basis and the objective function value would be 
decreased (when minimizing). When applied to profit rather than cost, the Dj's indicate the 
increase in profitability required to make the variable tie for a position in the basis. If the 
profitability is increased more it will enter the optimal solution and the maximum profit will 
be increased. 
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Sample Report 



DO.D/J 

NO £ SE D IG I T ._ .__ '_ . 

MANTISSA LENGTH 09 

SUBR- SET- 02 _ ..... 

VBLS TYPE NAME CURRENT COST REDUCED COST BASIS VALUE 

WBIN1 .030 .254 _ _ _.224-_ 

WBIN5 J50 ,015 .135 

ROWS TYPE NAME__ INCR_B_VALUE DECR B VALUE 

Y I ELD • Ol4- 
+FE 2..5.68_ 



+CU 

. . +MN : .S^k 

+MG 

-AL . .' .252 

RANGE+S I i .^ 8 5 



Report Description 

The report resulting from the DO. D/J agendum appears in two sections, each with appropri- 
ate report headings. 

The first section contains the structural variables which are not in the basis and the struc- 
tural variables which are uneconomic ally in the basis at either their lower or upper bound. 
The output sequence is the same as the input sequence. 

VARBLS: Indicates that the variable group follows. 

TYPE: Type W indicates the type of a structural variable which is either in the basis at its 
lower bound or is not in the basis. Type G indicates a structural variable in the basis 
at its upper bound. 

NAME: Contains the name of the structural variable. 

REDUCED COST: This field contains the reduced cost. For nonbasic variables (Type W) 
the meaning given above applies. For a basic variable at its lower bound (also type W), 
the economic interpretation is the cost of forcing one unit of the variable into the solu- 
tion by increasing the bound. For a variable in the basis at its upper bound (type G) , 
the economic interpretation is the reduction in the objective function (when minimizing) 
that would be obtained if one more unit of the variable were permitted in the solution. 
When applied to profit items, the opposite economic interpretation holds. The reduced 
cost will be shown to three decimal digits unless otherwise specified in column 19 of the 
DO. D/J agendum card. 

BASIS VALUE: This field contains the value that a variable must attain in order to econom- 
ically tie for entry in the optimal solution. It is computed as BASIS VALUE=CURRENT 
COST-REDUCED COST. 

Card Count: The output cards are numbered sequentially. 

The second section of the report provides the incremental or decremental marginal values. 
The incremental marginal value tells how much the objective function will increase cor- 
responding to an increase of one unit of the right-hand side (within a certain range). The 
range is not provided. The decremental marginal value indicates the increase in the objec- 
tive function value for each unit decrease in the right-hand-side constraint. 
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ROWS: Indicates that the row group follows. 



TYPE: Identifies the type of row: + for a ^ constraint row, - for a ^ constraint row, and 
for an equality row. If the row also has a range specification on the input cards, the 
word RANGE will precede the indicator. 

NAME: The row name corresponding to the logical variable is given. 

INCR B VALUE: The marginal value prescribing the amount by which the objective function 
would be increased if the right-hand-side restriction were increased by one unit. 

DECR B VALUE: The marginal value prescribing the amount by which the objective function 
would be increased if the right-hand-side restriction were decreased by one unit. 

Card Count: The output cards are numbered sequentially. 

COST. R 

Format 



Agendum 


Data 


Starting Sector 


Problem 


Decimal 


Comments 


Name 


Source 


Location 


Name 


Digits 




1 6 


7 


8 12 


13 18 


19 


20 80 


COST. R 








d 




COST.R 


D 


sssss 


nnnnnn 


d 





Agendum Name: COST. R 

Data Source: blank. - Cost ranges for all objective function coefficients will be produced. 
When Data Source is blank, so also will be Starting Sector Location and Problem Name. 

D - Cost ranges will be given only for those variables specified in the 
referenced file on the disk. 

Starting Sector Location: The address specified on the INPUT. C card which provides the 
names of the variables for which the user desires cost ranges. Only those variables 
in the file which are in the basic solution will have cost ranges calculated. 

Problem Name: Specifies the name of the file which was on the INPUT, card. If the name 
specified on the COST. R agendum card does not match the name on the sector specified, 
an error condition will result. 

Decimal Digits: The number of digits the user wishes following the decimal point. If blan T 
three digits will be given. 

Examples 

COST. R 

COST. R 6 
COST. RD03000SHIPNG2 
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Description 



Cost ranges have meaning only for structural variables that are in the basis. The ranges 
apply to the coefficients of the objective functions and indicate the upper and lower limit for 
prices which would yield the same activity level as the optimal solution. If the price were 
to go above the upper limit or below the lower limit, the activity level of that variable will 
change. It provides a good measure of the sensitivity of the optimal solution to price 
changes. 

Sample Report 



cos", a 

NO iSE DIGIT 0. ._. 

MANTISSA LENGTH 09 

SUBR SET 02 

COST.R NAME CURRENT COST HIGHEST COST HI-VAR LO-VAR LOWEST COST 

BIN2 .080 .089 MN BIN1 .017 

~ BIN3 " .170 .179 3!N5 MN ' " .160 

BIN4 .120 oW BIN1 MN _ .109 

ALUM " ~ .210 ~ ' .226 MN ". AL .189 

SILC0N .380 .467 MN GSI .148 



Report Description 

The report resulting from the COST. R agendum contains a heading line and a series of detail 
lines corresponding to all, or a subset, of the basic structural variables. The output se- 
quence will be the same as the input sequence. This page describes the report if COST. R 
follows a MIN 



COST. R: Indicates that the following report contains the cost ranges 

NAME: The name of the structural variable to which the cost ranges apply. The cost 
ranges apply only to variables in the basis at an intermediate level (type F variables). 
Cost range information about other variables is obtained from DO. D/J. 

CURRENT COST: Current value of the objective function coefficient. 

HIGHEST COST: The highest value the cost of the variable may take on and still yield the 
current activity level. 

HI-VAR: The basis type if at its upper bound, and the name of the variable which limits the 
highest cost of this variable; that is, if the cost of variable whose cost range we have 
obtained increases above the upper cost limit, its activity level will decrease, and the 
activity level for the variable named in this field will increase. * 

LO-VAR: The basis type if at its upper bound, and the name of the variable which limits 
the lowest cost; that is, the variable named in this field will increase* if the variable 
whose range we are obtaining goes below the lowest cost. 

LOWEST COST: The lowest value the cost of the variable may attain and still yield the 
current activity level. 

Card Count: The output cards are numbered sequentially. 



*If the variable is already at its upper bound, it will decrease rather than increase. 
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If COST.R follows a MAX. . . , then highest cost is equivalent to lowest profit, and lowest 
cost is equivalent to highest profit. 

CHECK. 

Format 



Agendum 
Name 


Data 
Source 


Starting Sector 
Location 


Problem 
Name 


Decimal Digit 
Field 1 Field 2 


Comments 


1 6 


7 


8 12 


13 18 


19 20 


21 80 


CHECK. 
CHECK. 


D 


sssss 


nnnnnn 


d l d 2 
d l d 2 





Agendum Name: CHECK. 

Data Source: blank. - Complete CHECK, of all rows and columns. If Data Source is blank, 
the Starting Sector Location and Problem Name field will also be blank. 

D - The row and column identification files will be obtained from disk. 

Starting Sector Location: The address specified on the INPUT. C agendum card which 
contains the names of the desired rows and columns. 

Problem Name: Provides the name to correspond with that on the INPUT. C agendum card 
which contains the names of the desired rows and variables. The name of the file on 
the disk must match the name on the CHECK, card or an error condition will exist. 

Decimal Digits: 

Field 1 - the number of digits desired following the decimal point for the upper limit, 
solution value, and lower limit field. If blank, three digits will be given. 

Field 2 - the number of digits desired for the row error. If blank, eight digits will 
be given. 

Examples 
CHECK. 

CHECK . D02500SHIPNGb6 
CHECK . D03000PROTEN15 

Description 

The CHECK, agendum is used in three different ways: (1) to provide a row check of the 
solution accuracy for the optimal solution, (2) to provide row summaries for selected vari- 
ables, and (3) internally by the system to ensure that the solution is proceeding within the 
defined accuracy. 

The report contains a line item for each row of the problem or a subset of the rows if so 
directed by a row identification file selected through the agendum card. The report contains 
the row name, the upper limit and/or lower limit for the row as indicated by the right-hand 
side or range constraint, the computed solution value, and the row error. The row error is 
computed by taking the sum of the products of the activity level times the coefficients and 
subtracting from the row limit. 
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A submatrix summary can be obtained by indicating the desired row and column identifica- 
tion (from row and column identification files) for the desired results. In this case the 
summation will occur only for the columns selected, that is, all points of intersection in 
the matrix. 

Sample Report 



CHECK. . . 
NOISE DIGIT 
MANTI SSA_LENGTH 09. 



-CHECK. — BOW NAME 
VALUE 
VALUE2 


_ UPPER LIMIT 


SOL. VALUE 
296.217""" 
296.217- 


LOWER LIMIT 
.000 
.000 


ROW ERROR 
"".00000000"" - 
.00000000 


YIELD 

FE __ .. 


2000 .000 
60 .,000 


2000 .000 
60 .000 




.00002000 
.00000010- 


CU 

MN 

MG 
AL 


100.000 
JfO.,000 
30 .000 


83.968— 
4*0 .000 
— " T9.960— 
1 500 .000 


1500 .000 


.00000050"" 
.00000010- 
"".00000010-"* 
.00001000 


SI 

* MAX ERROR ? __ .. 


300 .000 
.,00002000 


250 .000 


250 .000 


•000001U0 


Report Description 










CHECK: Identifies the type of report 









ROW NAME: Contains the name of the row. 

UPPER LIMIT: The upper limit of the row, either the right-hand-side value in an equal- 
to-or-less-than row, or the range constraint. 



SOL. VALUE: The sum of the products of the row coefficients times the activity levels for 
the selected columns. If all columns are used — that is, if a submatrix summary is 
not being obtained — the value should differ from the upper or lower limit by the amount 
of the row error (if any). 

LOWER LIMIT: The lower limit for the row — either the right-hand side for an equal-to- 
or-greater-than row or the lower limit of a range constraint. R-H-S objective function 
entries are reported in this field. 

ROW ERROR: The total accumulated roundoff error for the row, calculated as the differ- 
ence between the right-hand side and the sum of the products of row coefficients times 
the activity levels. 



CONTROL AND DATA MAINTENANCE 



In addition to the main functions of solving a linear programming problem — input, optimi- 
zation, output — a number of other functions are needed in the control language to make a 
complete system. These include various control statements and data maintenance facilities. 
Included in this section, in alphabetic order, are: 



ASSIGN 
ENDJOB 
ERASEp 
GETOFF 
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MAP. . . 

PAUSE. 

SAVE.p 

ASSIGN 
Format 
Card Columns: 



1-6 ASSIGN - agendum name 

7-8 Mantissa length: standard= as calculated; may be assigned values 05-18. 

9 Master tolerance: if used, all tolerances will be set to this value; as for all 
tolerances, it is the negative power of 10 desired. 

10 Element tolerance: standard = as calculated 

11 Pivot tolerance: standard = as calculated 

12 Feasibility tolerance: standard = as calculated 

13 Objective function tolerance: standard = as calculated 

14 Maximum error tolerance: standard = as calculated 

15 Iteration log output unit: standard = typewriter: 1 for typewriter, 4 for punch, 
9 for printer 

16 Reports output unit: standard = punch: 1 for typewriter, 4 for punch, 9 for 
printer 

17 Agenda log unit: standard = typewriter: 1 for typewriter, 4 for punch, 9 for 
printer. 

18-22 Sector address-DIM: where DIM table, programs and data terminate 

23-27 Sector address upper limit: usable for linear programming data up to this 
location 

28-29 Inversion frequency: standard = 15 
30-32 Iteration log frequency: standard = 1 
33-37 Common computation area sector address 
38-80 Comments 



Agendum Name: ASSIGN 

Mantissa Length: The size of the fractional part of a floating-point number. The system 
will initially calculate a mantissa length based upon the problem to be solved and the range 
of numbers encountered by the INPUT, agendum. As its solution proceeds, a recalculation 
is made after the CHECK, routine is executed prior to each inversion. If the error buildup 
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warrants it, the mantissa length is increased automatically. The dual and primal optimiza- 
tion routines continually check for error buildup; if it becomes excessive, control will also 
revert to a routine to recompute a proper mantissa length. If the user chooses, he may 
override this automatic mechanism and specify a mantissa length before the optimization 
is initiated by the MIN. . . or MAX. . . agendum card. To do this, the ASSIGN card should 
be placed after the INPUT, card. 

Master Tolerance: A single control is provided to allow the setting of all tolerances to a 
single value. It would give the same results as writing the number in this field in the next 
five fields. All tolerances are calculated by the system automatically, and generally are 
changed by the user only when the user has special requirements or special knowledge about 
the problem, or when he encounters difficulty in solving a particular problem. 

Element Tolerance: When any computed matrix value becomes smaller than the element 
tolerance, it is treated as "essentially zero" and given a value of zero. 

Pivot Tolerance: Whenever a candidate pivot element is found to be lower than tolerance 
level, it is treated as zero and ignored. 

Feasibility Tolerance: The feasibility tolerance is the amount a right-hand-side element 
may be negative and still be considered to yield a feasible result. 

Objective Function Tolerance: If during the "pricing" operation, the pricing criterion 
values becomes less than this tolerance level, it is considered "essentially zero" and 
rejected. 

Maximum Error Tolerance: If a row error exceeds this amount when automatically cal- 
culated by the CHECK, routine, the mantissa length will be increased prior to further 
iterations. 

Iteration Log Output Unit: The output unit to receive the iteration log; the standard unit is 
the typewriter. A digit 1, 4, or 9 will put the log on the typewriter, punch, or printer re- 
spectively. In addition, the user may get every iteration logged on the typewriter during 
MIN. . . or MAX. . . by placing console switch 1 in the ON position. 

Reports Output Unit: The standard output unit for reports is the card punch. They may be 
placed on the typewriter, card punch, or printer by assigning a 1, 4, or 9 respectively in 
this position. 

Agenda Log Unit: The agenda and comments cards in the agendum stream are always 
typed. To place this information on the printer or card punch, assign a 4 or 9 respectively 
in this position. 

Sector Address - DIM: This field is used to specify the lower limit where LP data may be 
placed following the DIM table and user programs. This field should be used to provide a 
lower limit for LP data following the work cylinders and the DIM table. The user must 
explicitly indicate available space above the work cylinders, DIM tables, and other user 
information before the LP system will use the area. He must also make an entry in the 
next field. 

Sector Address-Upper Limit: This field is used to specify the highest sector address 
usable by the LP system for data storage. If this and the prior field are not specified, the 
LP system will terminate if the problem data exceeds the work cylinders. Use of this field 
and the preceding field allows the user to file-protect part of the disk. 



73 



Inversion Frequency: The user may specify how often he wants to take inversions. If not 
specified, inversion will take place every 15 iterations. Only those iterations in which a 
variable leaves the basis are counted. 

Iteration Log Frequency: The user may specify the desired frequency of iteration logging. 
Standard logging frequency is 1 and, unless otherwise specified, it will be on the typewriter. 
If console switch 1 is ON, every iteration will be logged on the typewriter, in addition to the 
unit specified on the ASSIGN card. 

Common Computation Area Sector Address: A considerable saving of disk space can be 
achieved when many problems are being stored and solved by specifying a single common 
computation area. All computational results will be started at this location if an entry is 
made. 

Examples 

ASSIGN097bbbbbl91bbbbbbbbbb20005bbbbb 

ASSIGNbbbbbbbb999080001600006500 

ASSIGNbbbbbbb5 

Description 

The 1620-1311 Linear Programming System has been designed to allow simple hands-off 
operation. Many automatic features have been implemented to accomplish this goal. How- 
ever, there are times when it is desirable for the user to be able to control the system more 
directly. The ASSIGN agendum allows him to do this by giving him access to a number of the 
internal controls. These controls include the mantissa length, inversion frequency, toler- 
ances, output devices, and the reservation of disk storage. 

A number of things should be kept in mind: 

• The system will operate with standard settings unless otherwise directed by ASSIGN. 

• INPUT. C will reintialize to the standard settings except for agenda log unit, lower and 
upper available disk storage, and the start of the common computational area. 

• The contents of an ASSIGN card will stay in effect until another ASSIGN card changes 
these, or another INPUT. C agendum is executed which will initialize all values except 
those noted above. 

• When the SAVE, p agendum is used to save a problem, the contents of the ASSIGN card 
in effect will be saved along with the problem. If an INPUT. D agendum later refer- 
ences the saved problem, the ASSIGN controls which were saved will be restored. Thus, 
the problem can be reoptimized without another ASSIGN card being used. 

• Any or all of the fields on an ASSIGN card may be filled in. 

• Multiple ASSIGN cards are permitted and they may be interspersed among the agendum 
cards. 
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END J OB 



Format 



Agendum Name 


Comments 


1 6 


7 80 


END JOB 





Agendum Name: ENDJOB 
Comments: Any desired information 
Example: ENDJOB 
Description 

This card is used to terminate the linear programming run(s) . Control will be returned to 
the 1620 Monitor. 

ERASEp 

Format 



Agendum 


Storage 


Starting Sector 


Problem 


Comments 


Name 


Code 


Location 


Name 




1 6 


7 


8 12 


13 18 


19 80 


ERASEB 










ERASEB 


D 


sssss 


nnnnnn 




ERASEI 










ERASEI 


D 


sssss 


nnnnnn 





Agendum Name: ERASEB is used to delete a previously saved basis or a previously saved 
basis and inverse. 

ERASEI is used to delete a previously saved inverse (the basis is still 

retained) . 

Storage Code: blank - erases current problem basis and/or inverse 

D - erases basis and/or inverse which is at specified disk location. 

Starting Sector Location: A five- position numeric field which identifies the disk storage 
location where the referenced basis (and inverse) is stored. 

Problem Name: A six-position alphameric field containing the problem name which is 
checked against the problem name stored at the referenced sector location. 

Comments: Any desired information 
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Examples 



ERASEB 04300PRODMX 
ERASEI 08750DISTBR 

Description 

ERASEp is the means by which a previously retained basis (and inverse) may be deleted. If 
the basis will no longer be useful, it must be erased before re-solution of the problem. If 
just the basis has been saved (through a SAVE. B), an ERASEB is all that is needed. If a 
SAVE. I has been used to save both basis and inverse, the user has to decide whether the 
basis might still be of value. If not, then an ERASEB is used which eliminates both the basis 
and the inverse. The ERASEp agendum is not used very often, since a new basis and in- 
verse can be saved right over an old one. 

If a previously retained inverse is destroyed by another input or usage of the disk area, this 
destruction must be indicated to the system by an ERASEI before its attempted use in 
optimization. 

GETOFF 

Format 



Agendum Name 


Comments 


1 6 


7 80 


GETOFF 





Agendum Name: GETOFF 
Comments: Any desired information 
Example 

GETOFF 
Description 

GETOFF provides a means of interrupting a linear programming solution for a higher- 
priority job or to provide a restart basis if problem solution information (from iteration log) 
indicates that solution is not progressing satisfactorily. GETOFF is meaningful only if a 
basis has been established. It will generate basis cards (similar to OUTPUT except that 
values are not inserted): 

BASIS. 
VARBLS 

Basic structural variables 
SLACKS 

Nonbasic logical variables 

If sense switch 2 is ON, all cards will be bypassed until either a GETOFF or ENDJOB card 
is read. 
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MAP. 



Format 



Agendum 
Name 


Data 
Source 


Starting Sector 
Location 


Problem 
Name 


Comments 


1 6 


7 


8, 12 


13 18 


19 80 


MAP. . . 










MAP. . * 


D 


sssss 


nnnnnn 





Agendum Name: MAP. . . 

Data Source: blank. - The MAP. . . report is desired for the current problem. 

D - The MAP. . . report is desired for a specified problem which is 
stored on the disk. 

Starting Sector Location: The address specified on the INPUT. C card for the row and/or 
column identification files. 

Problem Name: The name of the problem specified on the INPUT. C card. This name must 
match the name of the problem at the Starting Sector Location, or an error condition 
will exist. 

Examples 

MAP. 

MAP. . . D02000ALLOYA 
Description 

The MAP report specifies the disk storage that contains the current problem or a previous 
problem. The report always contains ( 1) the starting and final sector locations occupied by 
the input, (2) the problem name, and (3) the number of entries in the row and column identi- 
fication files. In addition, if the input data includes both a matrix and a right-hand side, the 
report contains the current mantissa length and the number of sectors required for the in- 
verse and beta file computation area (computation file area). This area is in addition to the 
input area. 

If the report is requested before optimization for a problem with which no basis was input, 
the input area listed does not include that which would be necessary for saving a basis. If 
the report is requested after optimization, or if a basis was input with the problem, the input 
area does include the basis. 

Format of Report 

Columns 

1-6 MAP. . . - agendum name 

7 S 
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8-12 
13 



Starting sector location occupied by the input 

/ 



14-18 Final sector location occupied by the input 

19 b 

20-25 Problem name 

26-27 . b 

28-30 Number of entries in the ROW. ID file 

31-34 ROW, 

35-38 Number of entries in the COLJD file 

39-41 COL 

42-48 . MANSA- 

49-50 Current mantissa length 

51-59 -COMPUTE- 

60-65 Number of sectors occupied by the computation file area 

66-74 -SECTORS. 

Columns 42-74 appear only if the input includes a matrix and right-hand side. 
Sample Report 

MAP...D03600ALLOYC 

MAP...S03600/03605 ALL0YC».009RQW,0007CQl, t l!lANSA-10-C0MPUTE-000086-SECT0RS. 
MAP...DCM00ALL0YD 

MAP . . • S0A40Q/0V»0& -ALLQYft - JQ.Q0 RQW »Q.Q0 2 CO I „ 



The compute area varies with machine size. The illustration given above is for a 40K 
machine. For a 2 OK machine the area would be less; for a 6 OK machine it would be greater. 

The following chart was developed from the MAP. . . reports for several problems. It gives 
an indication of the number of sectors required for the files of problems of various size 
problems. 

Several of the MAP. . . reports were produced before optimization; several after. This is 
indicated by the columns titled Number of Sectors — No basis, and Incl. basis, respectively. 
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Problem 


Rows 


Cols. 


Number of Sectors 


Mantissa 
Length 


No 
basis 


Incl. 

basis 


Compute 
Area 


Share 21 


245 


352 


971 




8958 


10 


Share 5 


399 


344 


1175 




17156 


9 




63 


22 


188 




426 


11 




39 


13 


96 


104 


232 


11 




28 


13 


104 


110 


201 


11 




42 


107 


264 




361 


9 




6 


3 




25 


142 


9 



PAUSE . 
Format 



Agendum Name 


Comments 


1 6 


7 80 


PAUSE. 





Agendum Name: PAUSE. : period in column 6 is not required. 

Comments: Any desired information 

Example 

PAUSE. 
Description 

This agendum serves to stop the computer, and leaves it ready to continue LP processing 
when the start key is depressed. At this point, the next agendum is read and processing 
continues. It might be used between two LP jobs when it is desirable to look at the results 
of the first, before the second is started. If it were decided to conclude LP processing, an 
ENDJOB card would be placed in the reader, and upon depressing the start key, control 
would be returned to Monitor I. 
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SAVE, p 



Format 



Agendum 


Storage 


Starting Sector 


Problem 


Comments 


Name 


Code 


Location 


Name 




1 6 


7 


8 12 


13 18 


19 80 


SAVE. B 










SAVE. B 


D 


sssss 


nnnnnn 




SAVE. I 










SAVE. I 


D 


sssss 


nnnnnn 





Agendum Name: SAVE. B indicates that the current basis should be retained in disk storage. 

SAVE. I indicates that the current basis and its inverse should be retained 
in disk storage. 

Storage Code: blank. - Saves basis (and inverse) with current problem. 

D - Saves basis (and inverse) at specified disk location. 

Starting Sector Location: A five-position numeric field which identifies the disk storage 
location where the current basis should be stored. 

Problem Name: A six-position alphameric field containing the problem name which is used 
to identify the problem when referenced at a later time. 

Comments: Any desired information 

Examples 

SAVE. B04600PROB1 
SAVE. I04200BLENDG 

Description 

SAVE.p is a method for retaining an advanced basis in disk storage and, if desired, the 
latest inverse from which to begin a subsequent solution. If the matrix data values are 
revised between runs, it is not worth storing the inverse; but unless the changes are very 
significant, it is often valuable to store the basis for subsequent use. If the basis is not 
saved on disk or recorded on cards through OUTPUT or GETOFF, any subsequent re- 
solution must start from an initial basis of all logical variables. 

To be able to use BASIS. D it is necessary to have used SAVE. B or SAVE. I when that prob- 
lem was previously solved, or to have previously input a basis. 

In executing SAVE.p, the 1620-1311 LP system stores the entire communication region at 
the referenced location. Therefore the entire problem is effectively stored at the new loca- 
tion. This problem can then be prepared for re-solution through an INPUT. D agendum card. 
The ASSIGN controls that were in effect when the problem was saved will be restored by the 
INPUT. D agendum. 
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FORMAT SUMMARY 



The next three pages contain an agendum card summary (see "1620-1311 LP Agenda" for 
detailed explanation) . Following this is a two -page summary of data input cards (see "Data 
Preparation" for detailed explanation) . Below is an index of all the symbols used in the 
summaries. 

SYMBOL EXPLANATION 

Column name 



cccccc 
dl 
d2 



Number of decimal places in primary output value field 

Number of decimal places in secondary output value field 

Basis classifications 

W Variable basic at lower bound 
G Variable basic at upper bound 
F Variable basic at intermediate level 

Name of objective function row. If this field is blank the primary (first) 
objective function will be used during optimization. 

Right-hand-side name 

Problem name 

Row name 

Disk sector of data to be stored, referenced, or deleted 

Type or classification. 

Row classifications 

b (blank) identifies an objective function 
(zero) identifies an equation 

+ (plus) identifies a less-than inequality or a range 
constraint 

- (minus) identifies a greater-than inequality 

* (asterisk) identifies an omitted row 

Column classifications 

b (blank) identifies a structural variable to be included in the 
matrix 

* (asterisk) identifies a structural variable to be omitted from 

the matrix 

±uuuu. uuuuuu Only entry of a single limit right-hand side, or upper row limit of a range 
constraint 



ffffff 

hhhhh 

nnnnnn 

rrrrrr 

sssss 

t 
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AGENDUM CARDS SUMMARY 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 




COMMENTS 


1 


7 I ■ . ■ . 


13 l , , ■ i 


191,,, , 2 . 4 


31 |36 |4I |46 |5I 


|3o |o 1 |oo |73 ( ( 80 


t.t.X.E.Q, 


LiPy.6,2,0 


1 . , , i 


1 , , i , • 




1 1 1 l_i L..1 1 . 1 1 






1 . , . , 





1 1 1 


—J 1 1 1 , . , , , , , 


* 


1 


1 . , , , 


1 i ,,, , 




rAiMMCMT^ fiADh TflMMCklTC Tkl /*AI. 9 -Oft 

vwiwintN'» vhk-u • wiffnirtcN » o in Cum w ( ( 




1 


1 , , . i 


1 . i , . • i . . . . , 




III t 

— 1 1 1 1 


1.W.P.U.T,. 


ClS.S.S.S.S 


MFtliTliTli'Ylilt 


1 i i i i • i i , i i i 




IfllrV 1 — ttbtNIDUM. 1NWUT INDLCAIftEl CftBftX roLLfXtl 


I.N.P.U.T,. 


Dis,s.s.s& 


hhtiliiliiTtiT'i 


1 i i . i • i i i i i , 




ililr" 1 — MUtwinuiw . man I . , . , , , , 








1 , 




lit i 

— I 1 1 1 . . . . ... 


R.E.V.I.S.E 


i 


1 , , , . 


1 ,,,,,, 


1 1 1 1 


RErt/l^E. AGENDUM INPUT IhinTfiiTdta rjvt?n^ ctsiiaw 







1 ■ i i i 


i 




III i 

— i i — i i i i i . . . . 


M.l.N 


.ih.h.Vi.kk 




1 1 1 1 . 




MIlM A£EKirHllM*MTMTMll?( CiRTPPTTlUt PllNf TMSW 
pixiim w»^rauuim«i»iinixi'iij.iic ppjcuuvt r\m\^fiVM, 


M.l,N,.,.,. 






1 i , , , 


1 1 1 1 


Mlikl ACENDUiM. MlNIMiIZE (&B3EC.TilVE.njUCTIfllN 


i i i i i 


1 


i i i i i 





1 1 1 1 


1 1 t 1 i i i i • • , 


M.A.X 


-iK,K,K,K,K 


fif.f.*,?,? 


1 . ... , 


1 1 1 1 


MAX A&ENhUiM.MAVIMiIZF ABJFCTIil/t FUNr.T.lOKl 


M.A.X 




i , . . , 


1 , , , , • 




NAM AGFNDUiM.MA¥lMiI2F /2)MFr.TlA/E FUN/TOON. 


■ i i i i 


i 


i , , , , 




1 1 1 1 


i i i I ...... , 


0.U.T.P.U.T 




i , . . , 


iikil i ■ 


1 1 1 1 


0UITPUT AGiENinUM. I i 


fi.U.T.P.U.T 


Dis.-s.s.&.s 


hnY.TCh.n.-h 


<l«l<k ■ • ■ • 




pUlTDUT AGE iN HUM . (ftU ITPUT SELECTED .FTMJM .MSK , 




i i i i i 


I i , i i 


1 i i i i*i i i i . i 




1 1 1 1 


0.A..D./.J 


i i i i i 


i . , . , 


J«i 


1 1 1 1 


WlTyjAGEMttUM.0UTPiUT ALL REhUCiEO .C.OST.S. . 


d.<&...i>./.j 


D|S,S.&,S.3 




i«l ... 


1 1 1 1 


Mhty.T AGENDUM. 0UTPDT SELECTED «0M VP* . . 




1 1 1 1 1 


1 i i i i 


I ■ 




I i i i 


C.0.S.T...R 


I I 1 1 1 


1 , , ■ i 


J 




C&ST.R AGErtDUM.0UTtPOT ALLf.0STi RANGES . . 


C0.S.T...R 


Dta.s.s.s.s 




, , 




ttiST.B AdEWWJM.pUTiPUT SELFCTED FROM. DISK . 


i i i i i 


I i i i • 


1 i . i . 


i i i . i , 




i i i I 


C.H.E.CX.. 


1 i i i i 


i i i i i 


J.Kit. 




CHiECK- AGF.NDUM. CHECK ALL R0WS 




DlS,6.S.&.3 




IiiJl 




lECK . AGENiDUM. HD\tA. I D FTI F FBAM .DIS1C. . . 


■ i i i i 


1 i i i i 


i , . , , 


« ' 




i i i 






1 . . i » 


1 • • 




1 t 1 1 i i ii i i . 



AGENDUM CARDS SUMMARY (Continued) 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
N A ME 


VALUE 




COMMENTS 


', I , I I 


?l , . . . 


l3 l , , , , 


l9 l , . . , 2 ' 4 , . . , . , 


31 |36 |4I |46 (51 


|56 |6I ,66 ,73, 80 


E.N.DJ ,tf,B 


1 . . , , 


l.l.l 


1 1 1 1 ■• 


1 1 1 1 


TERMINATE. Si LP RUNS,. END JOB ACENDM .... 


■ i i i i 


1 . I 1 i 


1 . 1 1 . 


1 ■ ■ ■ I' 


1 1 1 1 


1 1 ' 1 


S,A,V,t,.,B 


DlA,a.».s.s 




Ill 11*1111 1 I 


1 1 1 1 


SAi/E AGENDUM. CURRENT BASIS ■ 


S.A,V,E..,B 


1 . . . . 


1 1 1 1 1 


1 1 1 1 l-l 1 1 1 ■ 1 


I I I I 


SAiVE AGINDiUM.CuRKiNT BASIS , 


S.A.Y.E...I 


Dlfi.S.&.S.S 




III 11*1111 1 1 


1 1 , 1 1 


SA.VE ACENDiUM. CURRENT RASJS # ilNVERSE. . . . 


5.A,v,e.i.ii 


1 , , , , 


1 , , , , 


1 1 1 ■ l-l 1 1 . 1 ■ 


1 1 1 1 


SAiVE ACEMDiUM. CURRENT BASIS 4. INVERSE. . , 





1 . . . . 


I,,,, 


1 ■ 1 1 I* ■ 


1 1 1 1 


1 1 1 1 . . . ... . 


E.R.A.S.E.B 


Dl 3,S,S,S,5 


h|n,y»,»,x,yi 




1 1 1 1 


■RlASE A&ENlDUM. RASlIS CANO INV«^; .... 


E.R.A.Sf ,B 


1 I I . I 


1 .11. 


1 1 ■ 1 I'l ■ 1 ■ 1 ■ | 




1 1 1 L 


ERiASE AC EH|DOK|. BASILS (AND INV^R*^) 


E.R.A.S.E.t 






-I I . ■ 1 . 


1 1 1 L 


ERpASE AuEtyDOM . lNy EKSt , 


E.e.A.S.E.I 


I.... 


1 1 . . . 


1 1 ' 1 1 


1 1 1 1 


ERASE A^EIJIDOM -IN/ERSE , 


i i i i i 


1 i i i i 


1 ■ , , , 


'■■■''■'■Ill 


1 L J 1 


i ' i ' ■ ■ 


T fcl V/ C P T 


• iK ,K, K,K ,K 


p r p p 1* 


1 1 1 1 '•' ' 1 ' 1 ' 


1 1 1 1 


INiVERT AGfNDUM.SpLVES SIWOLJAN^OUS _ _ 


— i i i i i_ 


1 i . . . 


I.... 




1 1 1 1 


i i i . . . 


w,c , i .(p.r.r 


1 . .. i . 


1 _ i_i , i.„.i 






CEJTOFF ACENDUW. Lp DJTERRJJPTpM 


— i — i — i i i_ 


i , . , . 


1 1 1 1 ■ 






1 l' 1 1 


M.A.P, 


I.... 


1 ... 1 






MAP AGENDUM. CURRENT PR0BLEW 


M.A.P 


Dis.sa.s a 








MAP AtENOqM. PROBLEM SELECTED FRCM DISK , 


■ i i i i 


1 . . . . 


1 ... 1 






1 I 1 I . . i . . i . 


P.A.U.S.E.. 


I11.. 


1 ... 1 






PAUSE , , 


i i i i i 


1.... 


III.. 






1 1 1 1 


■ i > i i 


1 . . . . 


1 . . . . 






■ ' 1 i ...... . 


■ i i i i 


1 1 . . 1 


1 , , , . 






1 i 1 


i > iii 


1 , . . . 


1 1 1 1 1 






1 I l 1 ■ ■ 


— i — i — i i i 


1 . . . , 


1 , , I I 






1 1 i Iii 


i i — i i i 


1 , . . . 


1 , , . I 






1 I 1 I . . . . . . , 


— i — i i i i 


I . ■ • 


1 1 1 1 1 






' i I 




1 . . 1 . 
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AGENDUM CARDS SUMMARY (Continued) 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 




COMMENTS 




11111 


l3 l , , , 1 


191 2 . 4 


31 |36 1 41 |46 |5I 


I 56 I 61 |66 ,73, _ 80 


AHTfu 

H,0 , ,1 ,l» ,N 




11111 


' ' ' ' '*' 


1 1 1 1 


A^ICN ACEVIDUM i i , 


' 1 ' ' 1 






' ' ' ' '*' ' ' ' ' ' 


1 1 1 1 


MANTISSA i LENGTH . i 


11111 


1 i'i i i 


11111 


' " "'' 


1 1 1 1 


MASTER TlOLERAMCE. i . . . , . 


1 1 1 ' ' 


Y 


11111 


' ' ' ' '*' ' ' ' ' ' 


1 1 1 1 


ELfcMtNT T|OLERANC£| , , , 


11111 


. v 
1 i i i "i 


11111 


' ' ' ' '* 


1 1 1 1 


PIlVOT TOLlERAKICE. i i 


11111 


Y 

1 i i i i* 




' " " ' ' 


1 1 1 1 


FEASIBILITY, TOLERAiNCE i , . 


1 ' ' 1 ' 


1 ' ' ' ' 


v. 

"1 i i i i 




| 1 | | 


O&tfECTIVE 1 FUNCTION TOLtfcANCG 


11 11 1 


' ' ' ' ' 


,Y 

I'ii ii 


' ' " ' ' 


1 1 1 1 


MAYIMUM lERROR (TOLERANCE 1 


J 1111 


11111 


• V 

1 iX i ii 


' ' ' ' '*' ' ' ' ' ' 


1 1 1 1 


OOlTPUT JIMltT- ITERATION LOG 1 , . , , , , , 


' ' ' ' ' 


' ' ' ' ' 


1 i 1X1 , 


''''''' 


1 1 1 1 


OUTPUT UNIT- REPQRTS 1 , , , , , 


- 


l,i,i 


1 i i i*i 






aqenda lqg unit i i 


_J_I L_X_1_ 


11111 


1 , i i iA 


y,v y V 

AIAiaiAi i a | i i i i i 


1 1 1 1 


jfcJLTOR ADDKh5S-fc.NO OF V\\A TAnlfclj PROGRAMS. PA t A 


11111 


11111 


11111 


■ V V Y Y V 

1 i i iAiAiAiAiXi i i 


1 1 1 1 


iinnrn HTCW. 1 It'/lTT — il P h/VTA 


' 1 ' ' ' 


1 ' 1 ' ' 


' ' ' ' ' 


1 , i . .*i . , .XiX i 




Xiif | V - ^ i i^^Jf^^'v^ t | 








i , i , i • , , . , , ,y 


XY i | i | 


ITEEAT\e>Ni .LOG FRBQUE.N«=V . 

r- 1 U | i ■ i i iii 










YXXlYX i i i 


COMMON coy\PuTATtqM AREA 5ECtfj>R AppRESS 












1 1 1 I 












1 i I i 












1 1 1 1 , 












1 1 1 1 , , 












1 1 1 1 












1 1 1 1 












1 1 1 1 












1 1 1 1 












i 1 1 1 ,i 












1 1 1 1 i , , , , ■ , 












I 1 1 
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DATA INPUT CARDS SUMMARY 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






win m u ii i v 


-I-J — I — i I 1 


Zl — I I 1 I 


l3 l i . , , 


I9| f j f 2* 


31 |36 |4I 


|46 |5I 


i 56 i 61 I 66 I 73 ?o 


R.0.W...I.D 


_l I i ■ ■ 


I.I.I 


' ' ' ■ ■'■ ■ ■ ■ ■ ■ 


1 1 


1 1 


W»W ID INDIcja.TOC.MW imckiTTCTrAnoii cabds tf^ww 


l?0.w.. .In 


Dis.s.s.s.s 


m i"n,T>i>i,vi 


i ■ t 


1 1 


1 1 


ID indicator . QESK 1 . . , 


__ i i i ii 


i , . . .t 




' 1 1 1 I'l 1 1 1 ■ ■ 


1 1 


1 1 


R0|W IDENTIFICATION , 


■ i i i i 


1 ... . 


1 1 1 1 1 


1 1 


1 1 


1 1 


—1- 1 1 " 


C,0,L,.,I,1> 


1 . . . . 


1 1 1 1 1 


1 ' 1 ■ I'l ■ ■ ■ 1 1 


1 1 


1 1 


Cd|L ID I NDICjATDg . Q^i. jtDSNTIFlCATt^U £Aftlfi f^fUL^W 


c.a,L...i.D 


Dls .a .s iS .s 


h hi mm m>m 




1 1 


1 1 


C0ILID IND10AToe.T>l6* 





i , . . .i 


Ciccccc 


1 ■ ■ ' I,'! ■ 1 1 ■ 1 


1 1 


1 1 


C<*IL I-DtNTUFlcATION | 


1. J_l 1 1 


111.. 


1 . . . . 


» ' ■ ' '•' ■ ' ' ' ' 1 


1 1 


1 1 


1 1 1 1 


M.A.T.C.l.X 


III.. 


1 . i . . 





1 1 — 


1 1 


MAITRTX IkiDlOATOR.MATRTX ELEMENT f»B» . 


MA.T.R.I.X 


Dis.-s.s.s.s 


Hr>,"v*n,-n,T-\ 






1 L 


MATRIX TNniCJCT^B. IDBK. 1 


■ i i i i 


ClCc.c.^.c 


UK ,*.Y. t,v 


-ix.x.x.x, • v. x.x,* 


- yyy v 1. yyyy \ yy 




MATPTX EUEMeNT 1 I 


..... 


Cic.c.c.c ,c 


M*.i*.h .i-.i- 


tlX.X.X.Xi •iX.X.X.X ,y,x 


t yyyy i. vyvy tyy 


— 1 1 


MA-TRTX ELiEMENT I i 


— i — i — i i — i .,, 




W 


t|X,X,X.X, • iXrXiXiXiX iX 


1 1 


tza 12.2 •» £ |*<*« 


MAlTRIX E.Ll£N\£NT I i ....... , 


..... 


JCCCC.t 


tih.f.Ki'.f 


iix.x.x.x, . .x.v.x.x.x 


1 1 


L 1 


MAfrf?IX EU6.MENT I , 


— i — i — i i i 


1 1 1 1 1 


■ tiii 








1 1 1 1 . . . . . . . 


r.T.R.srr. . 


Si i i i , 


1 . . , . 








FIIRST CKklHT HANOI SIDE. TNDICJAT^K .... 


F.I.K.S.T.. 


BtK.k.h.K.h 


1 , . . , 








FIIPST RIGlHT HAUDI SIDE lNDlClATtf*. . . . . 


N.E.X.T...B 


•ih.K.h.h.h 


1 . . . . 








AftDmoMAL lEfciHT HAND SIDE lNDffCAXOfc . . . 


F.I.R/S.T.. 




K|-*.-nv..-*.-r> 








FTic^T rtcwt Wand i srne: iNmciThR . P15K- . 


L— — 1 1 1 1 


I.... 


Mr-.lr.lr.lr.lr 


-IU.U.U.U. • .U.U.U.U.U.U 






R£«MT HAMD| SIDE. EL|EMCK1T | , 




1 . . . . 




'lU.u.u.u, • ,au, u,u,u,u 


-VVVVI-VV-VVIV-V 




R1IG.HT HANlb SIDE E1L.EN1E.NX 1 


—t 1 — ii 1 1 


1 . . , . 


1 . , . . 








1 1 1 1 


B.A.S.I ,S. - 


1 ... . 


Ii... 








BASIS INDIfCAT^R . ^A-SIS CARDS | , , 


B.A.S.I.S.. 


Dis.s.s.s.s 










BAS"*S DODffCAT^It. (DTLSte 1 . 


— i i i i i 


1 . . . . 


I , . , . 








1 1 1 l.ii.,,, 


V.A.I?.RL.S 


1 ... . 


i i i i i 








VAlPTAR\_ESl IkiDlCAIT^K 1 , 


— i i i i i 


1 . . . ■* 


qc.cccc 








RAS.TC STRlUCTOtAL VAetARL.ESIDEMTlPJCAT.l^M. 


i i i i i 


1 


l.i.i 








1 1 1 1 1 I ■ ■ ■ l_L_ 



DATA INPUT CARDS SUMMARY (Continued) 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 






COMMENTS 


'•■•ii 


_l J 1 1 l_ 


'3| . . . . 


.9, , , , , 2 . 4 


31 |36 |4I |46 |SI 


|56 |6I (66 I 73 ! i . . , , , 80 


si-.A.e.ie.s 




I.,., 





1 1 1 1 


&ltACkSXm>TCA.T£>R-| 


.... i 


1 . . . .e 


►IKK h.l-.h 





1 1 1 1 


NC>N- BASlCf UOGXCAUi VARiA3U£ S jtptHXttrlJpAT^o 











I , , . ,. 


1 1 1 1 


" 1 1 1 


E.W.D.AT.A 


i . . . . 


1 ... . 


I . , . 


1 1 1 1 


TqTCr4lNATiS& InPuTi PATA | 


E i < A f =i . i 


I.... 


1 . . , , 


1 . . 


1 1 1 1 


TqCMINATRS 3W"PUT| PATA 1 , . , ■ . , . 


• ■■■I 


I . . ■ . 


1 , ■ , ■ 


J _L i I 1^ 


1 1 1 1 


1 1 1 1 


■ ■ ■ i . 


i . . . . 


1 . . . . 


1 . . . , 


1 1 1 1 


• • 1 1 


■ i i i i 


1 . . . . 


1 . . . . 


1 


1 II 1 


I 1 1 1 


■ ... i 


i . . . . 


1 . . . . 


1 


1 1 1 1 


1 1 • 1 • 


■ . . i i 


1 . ■ . ■ 


1 . . . . 


1 


1 1 1 1 


1 1 1 1 . . ■ . . . . 


..... 


I . ■ . ■ 


1 . . . . 


1 . . . 


1 1 1 1 


1 1 1 1 i i i i iii 


ii iii 





1 . , . . 


1 


1 1 1 1 


1 ' l • 


i i i i i 


I . . . , 


1 . . . . 





1 1 1 1 


1 " 1 


. . . . i 


i . . . . 


1 . . . . 


1 . . . 


1 1 1 1 


1 1 1 1 








1 . . . . 


1 


1 1 1 1 


1 1 1 1 ...... . 












" l 1 1 












1 1 1 1 . . i . i . . 












1 1 1 1 i . i i . • i 
























1 1 1 1 . i • • • i • 












1 1 1 1 i I i i i i . 












1 l 1 1 












1 1 1 1 












1 1 1 1 i . 












1 1 1 1 












1 1 1 1 ■ 












1 1 1 1 ii 


* 1 ■ ■ ■ 










» « » « 



±vvw. www Lower row limit of a range constraint 

±xxxx. xxxxxx Data value of a structural variable 

± yyyy- yyyyyy Upper bound of a structural variable 

±zzzz. zzzzzz Lower bound of a structural variable 
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SAMPLE PROBLEM 



The sample problem is a blending problem. It is constructed to illustrate the use of the 
various agenda and data input formats in the 1620-1311 LP system. 

STATEMENT OF PROBLEM 

An aluminum alloy smelter wishes to produce a particular alloy at minimum cost. The alloy 
must, however, meet certain chemical constraints. The smelter has available various scrap 
materials and some pure metal stocks. There are five scrap materials which are stored in 
separate bins; each scrap material has been chemically analyzed. The inventory of each bin 
is known, as is the cost of each scrap material. The pure metal has a stated cost and can be 
purchased as required. Two of the bins contain powdered scrap; since these oxidize rapidly 
it is required that at least a certain amount from each of these bins be used in the solution. 

The original matrix is displayed on the next page. To assist in understanding the problem, 
consider these illustrations: 

BIN3, for example, has an inventory of 800 lbs. The end product must include 400 lbs. 
of this scrap. The composition of the scrap is 2% Fe (iron), 8% Cu (copper), 1% Mn 
(manganese), no measurable Mg (magnesium), 80% Al (aluminum), and 8% Si (silicon). 
The market value is estimated at $. 17 per pound. 

A ton of the end product is required. The composition of the end product (ALOY1) must 
be less than or equal to: 3% Fe, 5% Cu, 2% Mn, 1. 5% Mg. ALOY1 must contain at least 
75% aluminum. It must have between 12. 5% and 15% silicon. 

DESCRIPTION OF SAMPLE PROBLEM INPUT 

The sample problem consists of five segments which together use virtually every program in 
1620-1311 LP. Each segment is briefly explained below. (These explanations should be 
examined in conjunction with a review of the input forms which directly follow them. ) 

Segment 1 ALLOYA 

Input (a row identification, matrix element, and one right -hand-side element file) is read 
from cards and stored on disk. The solution minimizes the objective function, VALUE, for 
right-hand side, ALOY1. The solution basis is stored on disk. The solution values are 
reported. Check is performed on all the rows to determine the row computational error. 
DO. D/J is performed outputting structural variables at their lower and upper bounds, and 
marginal values for each row. 

Segment 2 ALLOYA 

The current problem files are revised: the row type of SI is changed from a range con- 
straint (positive slack) to a negative slack; BIN1 is eliminated from the matrix. The alum- 
inum (AL) content of BIN2 is changed to .77, and the inventory set to 3000. The BIN2 copper 
content (CU) is changed to . 03. Change the cost of BIN5 material in objective functions 
VALUE and VALUE 2 to $. 13 and -. 13 respectively, and lower its inventory to 1000. Raise 
the required amount of BIN3 to 500. Change the right-hand-side value of FE to 50, and indi- 
cate a lower bound of 275 for SI. Take BIN4 out of the basis and put AL into the basis. 



88 



STRUCTURAL VARIABLES 


RIGHT-HAND SIDES 




BIN1 


BIN2 


BIN3 


BIN4 


BIN5 


ALUM 


SILCON 


ALOY1 


ALOY2 


ALOY3 


( VALUE 
Objective J 


. 03 


. 08 


. 17 


. 12 


. 15 


. 21 


. 38 








FunCti ° nS ) VALUE2 
V 


-. 03 


-. 08 


-.17 


-. 12 


-. 15 


-.21 


-.38 








/ YIELD 


1 


1 


1 


1 


1 


1 


1 


2000 


2000 


2000 


1 FE 


. 15 


. 04 


. 02 


04 


. 02 


. 01 


. 03 


<gn 


<40 


£50 


\ CU 


• uo 


. uo 


. UO 


• uz 








<i fin 


— DU 


-DU 


Constraints / MN 


.02 


.04 


. 01 


. 02 


.02 






<40 


<40 


<35 


] MG 


.02 


.03 






. 01 






<30 


<30 


<35 


/ AL 


.70 


. 75 


. 80 


. 75 


. 80 


. 97 




>1500 


>1600 


>1500 


\ SI 


. 02 


. 06 


.08 


. 12 


. 02 


.01 


.97 


250 to 
300 


200 to 
300 


300 to 
400 


Inventory 
Bounds j (Upper) 


200 


2500 


800 


700 


1500 












\ Minimum 
r Required 
^ (Lower) 






400 


100 















Segment 3 ALLOYD 

Input from cards a reference column file for later use as selective output. 
Segment 4 ALLOYB 

Input a new column file. Input the revised matrix files from the disk, and input three right- 
hand sides. Maximize objective function VALUE2 for right-hand-side ALOY2. Save the 
solution inverse and basis. Output the solution values. Minimize objective function VALUE 
for right-hand-side ALOY1. Output the solution value. Input from disk the saved files, and 
minimize the objective function VALUE for right-hand- side ALOY3 using the saved basis and 
inverse. Output the solution values and perform a check on the problem. 

Segment 5 ALLOYC 

Input a set of files from previous problems saved on the disk. Minimize objective function 
VALUE using right-hand-side ALOY1. Perform a selected output and do Dj using the refer- 
ence file for segment 3, ALLOYD. Punch out a set of solution basis cards and return to 
monitor. The basis is erased because major subsequent changes to the matrix are expected 
to be made. 

The actual input sheets for the sample problem are on the next seven pages. Following these 
pages is the typewriter output which was obtained during the running of the sample problem. 
The sheets are annotated to describe what the information means. If different ASSIGN cards 
had been used, the output reports would have appeared on cards or on an online printer. 
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IBM 

JOB nd 727 



1620/131 1 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 



KEY PUNCH- PUNCH DEC. PT. IN COL. 24 



JOB tit, f SAMPLE PROBLEM • SEGMENT 1 anai y*t 



DATE_lZi^_ 



PAGE _i 0F_L 



INDICATOR 
CARDS 


COL.UM 
NAME 


N 


ROW 
NAME 


VALUE 






COMMENTS 


1 , i i i i 






l3 1 , , , , 


191,,, , 2 . 4 


31 1 I 36 , , i , I 41 , 


. . I 46 . . . . > 51 , . . 


|56 t |6I i ,66 _ |73 i 80 


3,4-,0,0,0.3 


2lO,0,T 


O, I 


3lfc ,0,0,0,3 


210,0,7,0,2^,9,0,2,^,0 


Ei5,l ,1 ,9,6,3,6,1 ,1 .3,6 


Oil, 0,2, ,. 


M0NIT0IE COLD SITART 1 1 1 


*.*,J,0,B, 


Si , , 


1 


1,,, 


i , , , 


• 


• 


J0B C0 lNT(?Di_ I.I 1 I 


*,*,X.E,Q, 


LlP,l .fo 


2,0 


1 




i , , , , 


1. ' 


EXECUTE, LP I (oZO |PR06»RA|M , | 


l.N.P.UJ,. 


CiO.2.0 


0,0 


AiLLd.Y.A 


i , , , ,• 


i , , , 


I , 


1 1 1 l l 


R,0,W,.,IJ> 


1 i i 


i 


I , , , , 




i , . , ,• 


l , , , 


1 1 i 1 1 





1 . i i 




ViA.L.O.E, 


i , , ■ 


i 


I , , , ■• 


WW IDlENTIFIClATIfcM IFTL.E I I 





1 i . 




V|A,L,U,£,Z 


I . . . 


i . . . 


I . . . 


1 1 1 1 1 


i i i i i 


1 , , , 


,0 


Y.IEL.D, 


• 








1 1 1 1 1 


i i i i i 


1 i i 


,* 


F I E , , , , 






1 i i i 1 * I i i i i i 


1 1 1 1 1 


i i i i i . 


i , , , 


,+ 


C|0, , , , 


• i 


1 . ■ ■ ■ 


1 r , , , , , 


1 1 1 1 — u 


i i i i i 




,+ 


W|N ( , , , 


1 i i i i • i i i i i i 


>■■■■' 


■■It. 


1 1 1 1 U 





i , , , 


+ 


Ht&, , i . 


1 i i > i ' 


1 i i i i'l 


1 1 1 1 I'l 1 I 1 , . 


1 1 1 1 U 


I 1 ,1 1 - 


1 i i 




A|L, , , , 


1 i i i i • i i i i i i 


1 , , , ,-, I , , , , 


1 1 I 1 I'l 1 1 1 I I 


1 1 1 1 1— 


1 1 l 1 — 1 — 


1 . . 




S|l, , , , 


• 


• 


' 


1 1 1 1 U 


M.A.T.R.I.X 


I < i 














■ , , i i 


B|I,H,» 




ViA.l.U.E. 


.,0,3, , , , 






MATRIX i ELEMENT, RLE. i i r 




B|I ,N,i 




YiX.E.l-.D, 


, , , ,1,. 


i ,2,0,0,., , , , , , 








BiI.N.I 




F|E, . ■ , 


i 1.5. . . . 










B,I,N,I 




C|U, , , , 


i 0,5. . , , 










BiW 




M|N, , , , 


i , , , ,.,0,2, , , , 








i i i i i — 


B,I,N,I 




M.&. , . , 


OA , , , 








i i i i i 


BiI.N.I 




A,L, , , , 


7,0, , , , 











6,1 N, \ 




s l x , , , , 


.,0,Z, , , , 








i i i t i 








-. 0,3, , , , 






i l l l i 


* * ' i — i — 






ViA.L.U.E", 


i , , , ...G.S, , , , 








— 






Y|I .E.L.D, 


i , , 


,2,5,0,0 






— i — i i i i 






, , , 


i , , , ,.,0,H, , , , 















C I U , , , , 


i . . . ...o, 5 . , , , 









gj^g 1620/ 13 1 1 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 

JOB no 727 ,,q R TIT , f SAMPLE. PROBLEM • SEGHEMT I anai yst hatp PAGE OF J7_ 



KEYPUNCH-- PUNCH DEC. PT. IN COL. 24 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






COMMENTS 


1 , i i i i 


7 


13 1 , . , , 


191,,, , 2 . 4 


*. . , . 136. , , ,. ,41, 


|46 _ . -51. 


|56 |6I . 166 . t 73 . 80 


..... 


■ 


M.N. . . . 


i . , .-...O.-f. , . , 


1 i i ' i * 


1 i i i i * i > i . 


1 1 I I 1 


..... 




W.6. . . . 


1 . , . ,-A3, , , , 


1 , , . .*.. i i , , 


1 1 i 1 i'i 1 1 II i 


1 1 11 1 


i i i l — 1 


B.I.N.2. . 


A,L, . , , 


I , , . ,.,7,5. , , , 


1 , , , i * , , i I , i 


Iiiii'iiiii. 


1 1 1 1 1 


, , ■ ■ ■ 


B|X,N,2, , 


S|I. , , , 


..afc . ■ ■ 


1 i i i i * l i i i i i 


• 


1 I I 1 1 




Bil.N.Z, .. 


ViA.L.0,6,2 


-i , , ■ ,.,0,8, , , , 


1 1 1 1 1*1 1 , " , I , 


1 i i i , * i i . . , . 


1 1 1 1 1 


i l_J I — 1 — 


BiI.N.3, , 


tfA.LD.E, 


. , , , ...1.7. . , . 


1 I , , l*i 1 1 , 1 i 


1 , , , i * i i i . i , 


1 1 1 1 1 


, , , . . 


M.N.3. , 




> ■ • .». 


1 .s.o.q. 


1 .f.o.o,* 


1 I I 1 1 


■ i i i i 


B|W , 


^E, , , , 


..o.e, , , , 


1 1 I 1 1 * 1 1 1 1 1 1 


1 1 I 1 1 * 1 


1 I I 1 1 


■ i ■ i i 




C|U, , , , 


1 0,*, , ■ , 


1 1 1 1 1 . 1 1 1 1 1 1 


1 1 1 1 1 * 1 1 1 1 1 1 


1 1 1 1 1 






W|M, , , , 


i 0,1. . , . 


1 1 1 I 1 * 1 1 1 1 — 1 — 1 — 


1 1 1 1 1 * ■ 


J 1 1 1 L_ 





BiI.N.fc. , 


, , , 


«,,,,. ,?,o. , , , 


1 1 , 1 I * 1 1 I 1 1 1 


■'II' II 


1 J 1 1 U 


i i ■ i i 


BiI.N.S, , 


S| x . , , , 


i . . , ,.,o,s, , , , 


1 1 1 1 1 . 1 1 1 1 1 1 


1 1 1 I 1 * 1 1 , I , I 


1 — 1 1 1 1_ 




P|I.N,3, , 


ViA.L.U.E.Z 


- . . . 


1 1 1 1 i-*J — 1 1 1 1 1 


1 1 1 , 1 * 1 1 1 1 I 1 


1 1 1 1 U 




BiE.N.4-. . 


V.A.L.U.E, 


.... ...),2. , , , 






1 1 1 1 U 




D T M li 


v.t p i r» 
rii-.t.j-.D, 


1 


7 fin 


1 1 1 1 U l u l . 1 I 1 , 1 , 






B.1.N.M-. . 


, , , 


I . . . , , , 










*lW. . 


C.U. . . , 


0,2, , , , 


1 I, 1 ,. I 1 , 1 1 1 


1 1 1 1 1 • 1 J 1 1 1 1 








MiN , . , 


o.z. . . . 




1 1 1 1 1 • 1 1 1 1 1 1 


1 1 1 1 L_ 




¥M , 


A,L, , , , 


7.5, . . . 






1 1 1 1 L 




•VM.H. , 


5,1, , , , 


, , J.Z, . , , 


I. .1 1 1 1 *.l 1 1 1 . J L 




1 1 1 1 1 




B.W, ■ 


ViA.L.O.E.Z 


-1 . . . 


1 , 1 , , * , 1 1 . 1 1 


1 I 1 1 1 . 1 1 1 1 1 1 


III I I 




B.1.H5. . 


ViA.LU.E. 


. . . 


1 1 . , , « . ii, i , 


1 1 , I , • 1 1 , 1 , 1 


1 1 1 1 1 




B|I,W,5, , 




. ■ • .1 


1 1, 5, 0,0, 


1 1 1 I 1 . 1 , 1 1 1 1 


1 1 1 1 1 




BiI.N.S. . 


^, , , , 


0.2. . , . 


Ill 11*1111 1 — L 


1 1 1 I I * 1 1 1 1 1 1 


1 1 1 1 1 




&I 1 N.5. , 




0.<o. . . . 


1 I i I lJLi — i i I i — i 


1 1 I I 1 • I I I I 1 I 


1 1 I 1 1 — 




B|W, , 


K«N. . . , 


o.z. . . . 




1 , 1 I I • 1 1 1 1 1 , 


■ I I 1 I 




Bil.N.S, , 


*C. , , , 


o,i. . , , 


1 I 1 1 I I I 1 — 




1 1 1 1 1 






A|L, , , , 


. . . 


1 1 1 1 1 * 1 I 1 1 1 1 


* ■ 11 1 * 1 11 1 1 ■ 


1 I I 1 t _ 



gjyi 1620/1311 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 

JOB NO. -I*! JOB TIT. F SAMPLE PROBLEM - SEGMEMT 1 AM a, vc T DATE 7/IO 3 np 7 



KEYPUNCH- PUNCH DEC. PT. IN COL. 24 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






COMMENTS 


1 

1 1 1 1 1 


7| 


l3 l , , , , 


191,,, , 2 . 4 


31 136 |4I 

1 1 1 1 1 L J [ 1_ 


|46 |5I 
. . 1 . . . . r . . . 


|56 (61 166 |73 80 




B,I,N5 , 


, . , 


1 , , , ,.,o,z. . . , 






1 1 — » 1 1— 

1 1 1 _ 1 1_ 




B,T.W,5, . 


ViA.L.U.E.Z 


-1 , . , ...1,5, , , . 




m 


1 1 1 1_ 


1 1 1 1 1 


AiL.u.m , 


tfi/U-.U.E. 


1 . , , , • £,l , , . . 


1 


m 


1 1 




Aii-.U.Ml , 


Y,I,£,L 1 D, 




1 




1 1 1 


1 1 t 1 1 


M.U.M. . 


F,E„ , , , 


.0.1. . , . 


1 




1 1 1 




A.L.U.M, , 


c,i> , , , 


.,0,1, , , , 


# 






1 1 1 1 1 


AiL.U.M, , 


A.L, , , , 


.... ,.,9,7, , , , 





1 




1 1 1 1 1 


A.L.U.M. , 


Sit, . , , 


.,0,1, , , , 








1 1 1 1 1 




A.L.U.M. , 


Y.A.L.U.E* 


-1 , , , ,.,?,), , , , 








1 1 l 1 1 





S.I.L.C.#JN 


YiA.L.U.E. 


3.S . , . 


1 


m 


1 1 l 1 1 




5.I.L,C.0,N 


Y.I.E.l.D, 


1 . . ,1 








1 1 1 1 1 

1 1 1 1 1 


1 1 1 1 1 


S.l.L.C.fc.M 


F,E, , , , 


1 0,3, , . . 









1 1 1 1 1 


SiI.L.C.M 


S,T, , , , 


.,9,-7, , , . 


1 





1 1 1 l 1 




SiI.L.C.^.N 


ViA,L,U,G,2 


-. , , . ...3,* , , , 





1 


1 l 1 1 1 


F.I.R.S.T.. 


BiA,L,0.Y,l 


1 , . . . 






, 


RIGHT jUAND 5 ( IDE ELEiMENT iFILE 1 




1 , , , . 


YiI.E.L.O. 


,2,0, 0,0,. 










1 1 1 , 1 


F.e, , . , 


■ . .6,0 










1 . . . 1 


c,u, . , , 


, ,1,0,0,. 








— 1 — 1 i 1 1 


1 , . . . 


*iN. . , , 


1 ■ .*o 










I , , , , 


M.G, , , , 


, , ,3,0 






1 1 1 1 i_ 




1 . . , . 


A.L, , . , 


,1,5,0,0,. 






1 1 I I U 


1 J 1 1 — 1 — 


1 . , . . 


Si I, , , , 


, ,3,0,0 


1 ,2,5,0,. 








1 1 1 1 1 











END 0F, INPUT iDATA 1 1 I 


A.S.S.I.G.N 


0,9,3,5, . 


i , , 


1 


Oil, , , ,. 






M.I.W.-.-.- 


•lA.LfiJ.Y.I 


V.A.L.U.t. 










S.A.V.£...B 


1 , 1 , , 


1 , , , . 








5 AVE .BASIS 1 1 1 1 


G.U.T.P.U.T 


1 , , 1 . 


1 . 1 > 1 












1 , . 1 . 


l 1 1 1 1 


1 , 1 1 1 * L 1 1 1 1 _l. 









IBM 



JOB N0._l5=l 



1620/1311 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 

JOB TITLE SAMPLE PROBLEM -SEooaeMT e,3 ANAI Y<;T 



r>ATF 7 ^'° PAGE J±_ of:JL 



KEYPUNCH- PUNCH DEC. PT. IN COL. 24 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






COMMENTS 


■.III! 


7 | I . 1 1 


l3 l I 1 1 1 


q. 2* 




|46 j5l 


|56 |6I |66 |73 80 
1 1 1 1 i 


C.0.S.T...R 


1 . . 1 I 


1 I 1 i.l. 




1 1 1 1 ■•' 1 ■ 1 1 1 


1 1 1 1 '•' 1 1 1 1 1 


1 1 — 1 1 L_ 


D.0...D./.J 


1 . ■ . . 


1 . 1 1 • 


C 

VI i . ■ .*...,. . 




1 1 1 1 '•' 1 1 11 1 


tfUTPUT 1 S ID£C1MALIS 1 I 


ft , . ,C,H 


EiC.K, 


TiP.U.T, .A 


N)|D, ,M,A,K,ti iKiC.,V,L 


oil .(JiN.s, • rr,{s, ,i,n.p 


UiT ,D,A,TA , , , , 


CARD 1 1 i i 


P.A.U.S.E.. 


_l i i i i 


I ■ ■ ■ i 




1 L. 1 1 1*1 1 1 1 


1 1 1 1 ■*' 1 ' 1 1 ' 


1 1 1 1 1_ 


R.E.V.L&E 


l i i i i 


1 i i i i 


1 " " 1 


' ' ' ' '•' ' ' ' ' ' 


1 1 1 1 '•' 1 1 ' 1 ' 


REVISE iCURREMTI PR<6BLE|N\ Pili_e.s I 


R.&W...I.D 


1 i i ■ i 


1 i i i i 


1 i i i i • i i i i i i 


1 1 L—l L*J 1 1 1 1 1 


1 1 1 1 ''■ 1 1 1 1 1 


L 1 1 1 U 


■ ■ i i i 


1 i . i i- 


C.T 


1 1 l l l • l l l l 1 L 


1 1 1 1 ■ * ' 1 1 1 1 1 


1 1 1 1 1 1 ' 1 ' 


1 1 1 1 _L 


t.O.L.. ,I,D 


1 . . . > 


i , . . . 


I ■ ■ ■ ■ * 1 


1'' ii*iii. . . 


1 1 ' 1 '*' ' 1 ' ' 1 


1 1 1 1_ 1_ 


■ ■ i i i 


1 . . , ,* 


Bawl, , 


1 1 1 1 1 1 1 


1 l 1 . 1 * 1 1 1 1 1 1 


1 1111 


1 1 1 1 U 


M /% "v O T* 

W.ft.T.R.i.X 


1 i ■ . ■ 


i i i i i 




1 1 1 1 1 * 1 1 1 1 1 1 


*'■!■'■ 


1 1 1 1 1 


i i i i i 




Ail 

HI*-, , . i 


1 i ■ i i • i 1 1 >■ ■ i ■ 


|3,000 


1 1 1 1 1 1 1 1 1 1 1 


1 1 1 I I 


i i i 


ft I N 2 


1*1 VJ, , , , 


i ..A3. 




' '' ' 




■ i i i i 


P|.»- 1 l\|,\> 1 , 


v/.A I ilC 
VIAil.|V>iC, 


lii i . • i iH i i i 


,100,0,., 




Ill U 


> i . i ■ 




T|X,C, 1— , w t 






1 .w.w ....... 




— i — i — i — i — i — 


&1.M.5. . 




-i . . , ...1.3. . . . 








F.I. B. ST.. 


BiA.L.Gi.Y.I 


U i i . . 




1 1 1 11*1 — 1 1 1 1 — 1 — 




i i i i r 


I I 1 .1. 1 


1 , , , , 


F,C. . . , 


i . .5.0. 


1 ' . . . * 




i i i i i 




i . . . . 


s,i, , . , 


i .2.7.5. 


^1 1.11* 




i i i i i. 


— I — i 1 — I — 1 — 

B.A&.T5.. 


1 . . . . 


J 1 1 1 — 1 — 


i ■ i i i ,* i i i i i — i — 


1 1 1 1 1 * 1 1 1 1 1 — 1 




i i i i i 




i . . . . 


1 .... 1 


1 1 1 L 1 * 1 1 1 1 1 1 — 


III 11*1 III . 1 




i i i i i 




1 , . . ,w 


Biim , 


1 l 1 i i • l i l • 1 1 




1 i l l i • I l l . 1 i 


i i i i i 


i i i i i 

SiL iAiC iKiS 


1 . . . . 


i i . . i 


1 1 1 1 1_JJ 1 — 1 1 1 1 


1 1 , 1 . . I 1 1 . 1 1 




i i i i i 


ill _i — i — 


1 , , ■ F 


AiL. . . ■ 


1 ll*i i*l il i 1 l 


1 1 ■ ■ ■ * 




i ii i i 


E.ftF, , , 


i i i i • 


I i i i i 






III l l • l III 1 i 


END (4FI INPUT 1 DATA 1 1 1 


I.M.P.U.T.. 


:io.4-.«f.o.o 




1 l ■ i i • i i i l l l 


1 1 ....... 1 . . 


1 1 i i i * i i i ■ i i 


STOCE 1 CGLOWU4 1 PILE | f&H \ i 


C.0.L...I.D 


i ii — i — i — 


1 . . . . 


Ii. i i » l 1 1 1 1 l 








■ i i i i 


i i i i i 


AiL.un , 








1 1 1 1 1 


i i i i i 


i . i i . 


S.I.L.C.0.M 


1 ... . •! i 1 1- 1 I 


1 1 111*1111 — 1 1 




1 1 I I I . 



gjyg 1620/ 131 1 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 

JOB NO 1^2 JOB TITI F -SAMPLE PROBLEM -SEGivibmt H- amai yct DATE_2^_ PAGE 5 



KEYPUNCH-- PUNCH DEC. PT. IN COL. 24 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






COMMENTS 


1 


7| , . . , 


'3. , , , , 


191 2 . 4 


31 1 . i 36 . . , , i 4 ', 


i 46 ... i 51 . 


|56 |6I m |66 _ |73 t 80 


E.C&F, , , 


I , , , . 


1 , , . , 


1 ' i i i * i 








ENfc API INPUT I DATA i i , 


I,N,P,U,V,. 


Cl0,2,«3.O.0 


A.L,u,a>,YB 


1 , , , 










R.fc.W,. ,l,t> 


D.0,2.0.0,0 


AiL.L.O.YA 


1 i i i i * i 1 i I , i 


i , , , ,. 


i , , , ,. 


f?0W ID. FILE. .DISK iREFCeeiNCC , 


C,0,L,. ,1,0 


1 i i . i 


. . , , , 


1 i i i . * i i i i , , 


i , , , 




l 1 l l I 





1 , , , , 


3iI.M,l. . 


1 i i i i * i I , i i i 


i . 





1 1 1 I 1 


, i , , i 


1 , i i i 


3.i,N,e, . 


— 1 — i i I i i i i i_ 


i , , , ,. 





1 1 1 1 l 


, , , , , 


..... 


Bil.N.3, , 


1,111*1,111, 


. , 


I,,,,.,,,,,, 


I 1 1 1 1 


,11,1 


1 , . . . 


*il,N,4-. , 


1 i ' i i * i 


i , , , ,. 


i 




iiii, 


1 , . . , 


W.N.S, , 


1 ' i i i * 


i 


i ,...*,.... . 


1 1 1 1 1 


, . , . , 


1 i i i i 


A|L,U,M, , 


1 i i i i • i i i i i i 


i , , , , 


i 


11.11 





1 . ■ . ■ 


S.t.L.C.tf.N 


1 i i i i * i 


i , , , , 


i . , , ,., 


1 1 1 I 1 


M.A",T,R,I,X 


1*0,2,0,0,0 


A|L,L,0,Y,A 




. , , , , 





MATRIX 1ELEMEN1T PILE (DISK RE|FERENC6| 


F.I.fc.S.T,. 


3A,L.,0,Y,l 


1 i i i i 


— 1 l i i i_*_l 1 l l i i 


i , , , , 


- 


1 II I 1 





1 , . , ■ 


Y|I.E,L.,D, 


i2,0,0,0. ., 


i , , . ,. 


i , , , ,. 


i 1 i i I 





l , , , , 


F,E, , , , 


I ■ .4.0 


i 





. 1 1 I l_ 




1 , , , . 


C|U, , , , 


i .'.o,o 










1 , , , , 


M|M. , , , 


i ■ ,V> 








• , i , , 


l . . , , 


, , , 


i . .3,0 











1 , , . . 


A,L. . , , 


.'.5,0,0 










I , . , , 


S.I. , . , 


I ,3,0,0,* 


i ,2,5,0,. 








.lA.Ld.Y.fc 


i i , , i 












i , . , . 




12,0.0,0,. 










l . . . . 


FiE. . . , 


i , ,V> 






1 1 1 I L_ 


i i i , i 


1 . . . . 


«|U, , , , 


I , ,6,0 










..... 


-1iN. . , . 


1 , ,4,0,. 








_i t i i i 





*I.G. . . . 


i , ,3,0 










. , , , , 


, , , 


.i,t,o,o, 






1 1 1 1 L_ 




. , , , . 


, , , 


, ,3,0,0 


, ,2,0,0,. 




1 1 1 I 1 _ 



IBM 

JOB NO. 



1620/1311 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 

727_ job TITLE SAMPLE PROBLEM -SE6MLNT H- amai y^t 



DATE _Z^i£__ PAGE OF_Z_ 



KEY PUNCH- PUNCH DEC. PT. IN COL. 24 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 










'.iiii 


71 ... . 


'31 , , , , 


»i . . . ? 


31 1 I 36 . . . . I 41 . 


i 46 > 51 . 


|56 ' t |6I _ |66 _ |73 § 80 


N,E,X,T,.,B 


.iA.L.U.Y.3 


1 , , , , 




i 





1 I 1 1 1 




i , , , , 


Yil.E.L.D, 


,2,0.0,0.. 








I I 1 1 1 


■ i i i i 


i , , . , 


FiE. . . . 


i . .5.0. • 








1 1 1 1 1 


1 i i i i 


i , . . . 


CiU, . , , 


i . .4..0,. 


, i , 






> I 1 t 1 


i . , , . 


4,N, , , , 


. , .3,5, 




i , , , ,. 


1 1 1 1 1 


i 1 i i I 


i ... i 


%C. , , , 


i , ,3,5,. 





, , , , 




1 i i i i 


1 i ■ ■ i 


AiL. , , , 


i 1.5,0,0,., 








1 1 1 1 1 


i i i i 1 


i , . , . 


8.1. . . . 


, .f.0,0, 


1 ,3,0,0,. 


i , , ,, * , i , . i , 


1 1 I I 1 


B.A.S.I.S.. 


DiO.2,0,0,0 


A.l.L,<*,Y.A 


, ,,,,,, 


i , , , , 


i i i i i • 


1 1 1 1 1 


E.4.F. . . 


1 , . . . 


. , , , . 


, , , , , 







END 0F .INPUT iDATA , , , 


A.S.S.lAN 


1.0.5. . , 


I , ,i, , 


. , , ,1.2,0,0.0, . .0 


Oil, , , , , , 





1 i I i i 


M.A.Y, 


..A.L.a.y.2 


VlAJL.U.E.Z. 





, , 


i i i , , • i i , i . i 


i i i ii 


S.A.V.&...1 


fill. 0,0.0 .0 


l.W,V,E,R,S 








i i . i i • i i i i . i 


SAVE I.NVERSE I AND BAi&IS *N D,1SK , 


fl.U.T.P.U.T 


1 i , , . 


i . , . , 


i . 








iiii i 


M.I.N 


.iA.L.fl.Y.1 


V.Ai,U.E, 


• i i i ' • 







i i i i i 


0.U.T.P.U.T 


1 i i , . 


i . . . 


i i i , , • . . . , , , 






i i i i r 


I.N.P.U.T,. 


Di 1,0,0,0,0 


IiN.V.E.R.S 


i i . . . • . . i . , . 






INPUT .PROBLEM. DATA iFRflM PJSK , 


M.I.N. 


..A.L.0.Y.3 


V.A.L.U.E. 


— i — i — i i i i i i_j 






1 1 1 1 L_ 


0.U.T.P.U.T 




i , . . , 


1 L—l 1 l_Lj 1 1 ■ ■ 1 






i i I i i 


C.H.E.C.K.. 


i ■ . . . 


i ... i 


1 1 . 1 1 • . . I 1 . . 






t 1 1 1 l 1— 


l.N.P.U.T,. 


010,3,6,0,0 


A.L ,L,0,Y,C 


1 1 I 1 1 • 1 1 1 1 1 1 






START \<&F SEGMENT \B i . 


R,<&,W.,1,D 


PiO.2 .0.0.0 


A,L,L,a.Y,A 


1 1 1 1 1_1J 1 1 L__l 1 






i i i i i 


C.0.L...I.B 


DiO.Z. 1.0.0 


A.L.L.0YB 


1 1 1 1 1 • 1 1 1 , 11 






l l l i i — 


M.A.T.R.l.Y 


DiO.Z ,0.0.0 


A.l.l.tf.Y.A 


. 1 J 1 . 1 J_i_J 1. 1 L. .1 1_. 






i ii I i_ 


F.I.R.S.T,. 


D|0.2 .0.0.0 


AiL.L,0,Y,4 








I 1 1 1 L_ 


B.A.S.I.S.. 


Di 1,0.0,0.0 


l.N,V,EJ?,S 








1 1 1 1 L. 


F.Q5.F, , , 


1 . . i > 


I , i , , 








END 0F, INPUT ,X>ATA , i , 


A.S.S.I.G.N 


0,9.3,5, , 


1 , ■'■ , 


, , . , . . .0 


0,1 




1 II 1 L_ 



gjyg 1620/1311 LINEAR PROGRAMMING SYSTEM 

INPUT DATA FORMAT 

JOB mh 727 ,ior t.ti p SAMPLE PROBLEM - SEGMENT 5 AMfll y * r nATF 7/l0 page J_ f_L 



KEYPUNCH^ PUNCH DEC. PT. IN COL. 24 



INDICATOR 
CARDS 


COLUMN 
NAME 


ROW 
NAME 


VALUE 






COMMENTS 


' 


7 I , , , , 


'31,,,, 


19, , , , , 2 . 4 


31 1 I 36 , , . . I 41 . 


|46 p5l 


|56 ,61 |66 _ ,73 _ 80 


M.I.N. 


..A.l.fil.Y.I 


V.A.L.U.E. 


1 . , , ,. 


1 , , , 


1 , , , 


1 1 1 1 1 


0,U,T,P,U,T 


PiO.M-,4-,0,0 




• 


1 , . , ,• 





1 1 1 1 1 


P.0,.,D,/,J 


D.O.^f.0,0 


AAL.tf.Y.D 


1 , 


1 . . . .• 


1 . . . ■• 


1 1 1 1 1 


C.C&S.T,.,R 


D.0.4-.4-.0.0 


A.L.L.fc.Y.D 


1 . . . 


1 ■ . ■ 


1 . . . 


1 1 1 1 1 


CET.CJ.F.F 


1 i i > i 


i . , , , 


I.i..- 


1 , 


1 . . . ■• 


PUNCH I0UT SdlLOTlOM l BASIS OiM CARDS i 


E.R,A,S,E,B 


1 , , i i 


i , , , , 


• 


1 . . . .• 


1 , , , ,• 


ERASE. iBASTS l l i i 


MAP,-,.- 


D.0.2.0.0.0 


AiL.L,(5,Y.A 


1 . . . 


1 . . , ■• 


1 . . . 


l 1 1 II 


M.A.P - 


DiO,2, "2.0,0 


AiLL&YB 


1 . . . 


1 , . . ,• 


1 . . . ,• 


1 1 III 


M.A.P 


DiO.3,6,0,0 


AL.L ,05 Y ,C 


1 . . . ,• 


1 . . . ,• 


1 . . . 


1 1 1 l l 


M.A.P. 


DiOAH-.O.O 


A.L,L,0,Y.D 


1 . . \ . • 


1 • ■ ■ .• 


1 . . . . •. 


1 1 1 l 1 


E.N-.D.J.0.B 


1 • . . i 


l i , , > 


1 ■ ■ ■ ■• 


1 , 


1 ■ ■ . 


END |0F PROBLEM l i i 


i i i i i . 


l , , i i 


i i i i i 


* 


1 , , , I'l 1 , . , , 


1 , ■ . 1*1 


I i I i i 


i i j — i — , — 


i . , , , 


i , , , , 


1 ■ ■ ■ 


1 ■ ■ . 


1 . . . 


I 1 l i i 


. . i i i 


i . . . . 


i . . . , 


■ . . . 


1 , , , 


* 


i i i i i 


i i i i i 


i . . . . 


i , i , , 


1 





* 


i I I i i 


i , — , — t — i — 


1 , , , , 


I i . . , 








l i l i r 





1 , , , , 


i , , . ■ 








iii.i 


i , — , — i — i — 


i . , ■ , 


i , . , , 








i i i i i 


_i i i — i — , — 


i . . . . 


i , , , , 








1 1 1 1 L 


i , — , — i — i — 


I , , , , 


i , , , , 








1 1 1 1 1— 


i i i i i . 




i , , , , 








l i I i I 


> i — i — , — i — 


i . , . . 


i . , . , 








i i i i i 


i i i i i 


1 , , i i 


i , , , . 








i I I i i 


■ i i i i 


i . , . . 


_l L_J 1 1 








ill I i 




i , . . . 


1 11.. J , 








I I 1 i i 


i i — i — > — i — 

. . . i i 


i . . . . 


1 . , . , 








I l l l i 


i i , — i — , — 


1,.,. 


1 , , , . 








i i I I i 


i , i i i — 


i , 1 L.._._ 


1 . , , , 








i 1 J 1 U 



TYPEWRITER OUTPUT - SAMPLE PROBLEM 



ttJOB 5 

ttXEQ LP1620 
EXECUTION 

NOISE DIGIT 
MANTISSA LENGTH 08 
SUBR SET 02 

UPPER LIMIT DIM = 07799 
LP1 620- 1 62 1 



JOB CONTROL 

EXECUTE LP1620 PROGRAM 



INPUT .C02000ALL0YA 

INPUT, 

ROW. ID 

MATRIX 

F IRSTB 

ASSIGN0935 1 
M IN . . . .AL0Y1 VALUE 
LP! 62 1 TO INVERT 

NOISE DIGIT 
MANTISSA LENGTH 09 
SUBR SET 02 
LP 1 62 1 TO DUAL 
NOISE DIGIT 
MANTISSA LENGTH 09 
SUBR SET 02 



001 



ITER 


VARBL. 


EXIT 


VARBL. 


FNTR 


CNT 


OBJ. FUNCTION 


0001 






BIN1 


U.B. 


003 


86 .000 


0002 


YIELD 


U.B. 


BIN2 




003 


1 90 .000 


0003 


SI 


U.B. 


SILCON 




005 


230.879 


0004 






B IN4 


U.B. 


004 


243 .010 


0005 






BIN1 


L.B. 


004 


250.373 


0006 


AL 




ALUM 




004 


288.007 


0007 


FE 




BIN4 




001 


295 .298 


0008 


MN 




BIN3 




001 


296.216 


FEASIBLE 












OPTIMUM 












SAVE.B 














OUTPUT 















SAVE BASIS 



NOISE DIGIT 
MANTISSA LENGTH 09 
SUBR SET 02 



97 



OUTPUT 
BASIS. 
VARBLS 



SLACKS 



MANTISSA 09 
TYPE NAME^ ACTIVITY LEVEL 
FBIN2 665.343 
FBIN3 490.253 
FBIN4 424.188 
FALUM 299.639 
FSILCON 120.578 
TYPE NAME ACTIVITY LEVEL 
296.217 
296.217- 
.000 



TOLERANCES 05 03 03 03 03 



output UfoKr 



FVALUE 
FVALUE2 
OGYIELD 
+WFE 
+FCU 
+WMN 
+FMG 
-WAL 
+GSI 

CHECK. 

NOISE DIGIT 

MANTISSA LENGTH 09 

SUBR SET 02 

CHECK. ROW NAME 
VALUE 
VALUE2 
YIELD 
FE 
CU 
MN 
MG 
AL 
SI 

* MAX ERROR = 
COST.R 

NOISE DIGIT 
MANTISSA LENGTH 09 
SUBR SET 02 



16.032 

1 .040 
50.000 



UPPER LIMIT 



2000 .000 
60 .000 
100.000 
40 .000 
30 .000 

300 .000 
.00002000 



SIMPLEX MULT 



.014 
2.568 



.544 

.252- 
.485- 



SOL. VALUE 
296.217 
296.217- 
2000 .000 
60 .000 
83.968 
40 .000 
19.960 
1500 .000 
250 .000 



LOWER LIMIT 
.000 
.000 



1 500 .000 
250 .000 



ctteck Rcfeftr 

ROW ERROR 
.00000000 
.00000000 
.00002000 
,00000010- 
.00000050 
.00000010- 
.00000010- 
.00001000 
.00000300 



COST.R NAME 


CURRENT COST 


HIGHEST COST 


HI-VAR 


LO-VAR 


LOWEST COST 


BIN2 


.080 


.089 


MN 


BIN1 


.017 


BIN3 


.170 


.179 


BIN5 


MN 


.160 


BIN4 


.120 


.147 


BIN1 


MN 


.109 


ALUM 


.210 


.226 


MN 


AL 


.189 


SILCON 


.380 


.467 


MN 


GSI 


.148 



98 



DO.D/ J 

NOISE DIGIT 
MANTISSA LENGTH 09 
SUBR SET 02 
VBLS 



ROWS 



vo* d/i retort 



1 Yrt NAnfc. 


CUKRhNI CUbl 


REDUCED CUST 


BASIS VALUE 


WD l IN 1 




OCT/. 




WBIN5 


.150 


.015 


.135 


TYPE NAME 


INCR B VALUE 


DECR B VALUE 




OYIELD 


.014 






+FE 


2.568 






+CU 








+MN 


.544 






+MG 








-AL 




.252 




RANGE+S 1 




.485 





* CHECK OUTPUT AND MAKE REVISIONS TO INPUT DATA 

PAUSE. 

REVISE 



COMMENTS CARD 



REVROW 
REVCOL 
REVMAT 
REVFST 
REVBAS 
REVSLK 



INPUT.C04400ALL0YD Seenwr 5 

INPUT. 

COL. ID 

INPUT.C02800ALL0YB sgtME^H 

INPUT. 

ROWDSK 

COL. ID 

MTXDSK 

FIRSTB 

BASDSK 

ASSIGN105 1 12000 001 



99 



MAX....AL0Y2VALUE2 
LP1621 TO INVERT 



NOISE DIGIT 
MANTISSA LENGTH 10 
SUBR SET 02 
LP1621 TO DUAL 
NOISE DIGIT 
MANTISSA LENGTH 10 
SUBR SET 02 

ITER VARBL. EXIT 

0001 

0002 

0003 BIN3 

0004 AL 

0005 CU 

0006 BIN5 
FEASIBLE 
OPTIMUM 

SAVE • ID10000 INVERS 
OUTPUT 

NOISE DIGIT 
MANTISSA LENGTH 10 
SUBR SET 02 
OUTPUT 
BASIS. 
VARBLS 



TYPE 
FB 
WB 
FB 



SLACKS 



NAME 
N2 
N3 
N4 
FALUM 
FSILCON 
TYPE NAME 
FVALUE 
FVALUE2 

OGYIELD 

+WFE 

+WCU 

+FMN 

+FMG 

-FAL 

- WSI 



VARBL. 


ENTR 


CNT 


OBJ. FUNCTION 


SI 


L.B. 


003 


499.560- 


BIN5 


U.B. 


003 


423 .560- 


MN 




003 


335 .755- 


BIN5 




003 


354.631- 


AL 




001 


371 .486- 


B IN4 




001 


373.686- 



MANTISSA 10 
ACTIVITY LEVEL 
1 90 .678 
500 .000 
233.051 
961 .864 
1 14.407 
ACTIVITY LEVEL 
373.686 
373.686- 
.000 



22.712 
24.280 
54.619 



TOLERANCES 



output Re Port 

05 05 05 05 05 



SIMPLEX MULT. 



.264- 
2.843 
2.754 



.208- 



100 



M IN , . . .AL0Y1 VALUE 
LP1621 NEW RHS 
NOISE DIGIT 
MANTISSA LENGTH 10 
SUBR SET 02 
LP1621 TO DUAL 
NOISE DIGIT 
MANTISSA LENGTH 10 
SUBR SET 02 

« m L - EX,T 

0002 ALUM 

0003 AL 

0004 MN 
FEASIBLE 
OPTIMUM 
OUTPUT 

NOISE DIGIT 
MANTISSA LENGTH 10 
SUBR SET 02 
OUTPUT 
BASIS. 

VARBLS TYPE NAME 
FBIN2 
WBIN3 
FBIN4 
FBIN5 
FALUM 
FSILCON 
SLACKS TYPE NAME 
FVALUE 
FVALUE2 
OGYIELD 
+WFE 
+FCU 
+WMN 
+FMG 
-WAL 
-WS I 

INPUT. D10000 INVERS 
INPUT. 

MIN....AL0Y3VALUE 
LP1621 NEW RHS 
NOISE DIGIT 
MANTISSA LENGTH 10 
SUBR SET 02 



cu 

ALUM 
BIN4 



ENTR 



003 
002 
001 



OBJ. 



273.623 
287 .986 
292 .607 



MANTISSA 10 
ACTIVITY LEVEL 
574 .696 
500 .000 
469.636 
130.972 
206.680 
1 18.016 
ACTIVITY LEVEL 
292.607 
292 .607- 
.000 

23.441 

1 1 M9 



TOLERANCES 05 05 05 05 05 



SIMPLEX MULT 



.036 
2.324 



.769 

.273- 
.501- 



101 



10 

EXIT 
U.B. 



LP 1 62 1 TO DUAL 
NOISE DIGIT 
MANTISSA LENGTH 
SUBR SET 02 
ITER VARBL. 

0001 BIN4 

0002 BIN2 

0003 BIN1 
FEASIBLE 
OPTIMUM 
OUTPUT 

NOISE DIGIT 
MANTISSA LENGTH 10 
SUBR SET 02 
OUTPUT 
BASIS. 

VARBLS TYPE NAME 
WBIN3 
FBIN4 
FALUM 
FSILCON 
SLACKS TYPE NAME 
FVALUE 
FVALUE2 
OGYIELD 
+FFE 
+WCU 
+FMN 
+FMG 
-FAL 
- WSI 

CHECK. 

NOISE DIGIT 
MANTISSA LENGTH 10 
SUBR SET 02 

CHECK. ROW NAME 
VALUE 
VALUE2 

YIELD 

FE 

CU 

MN 

MG 

AL 

SI 

* MAX ERROR = 



VARBL. 

BIN1 

FE 

BIN4 



FNTR 



CNT 
002 
001 
001 



OBJ 



FUNCTION 
368.091 
369.010 
369.205 



MANTISSA 10 
ACTIVITY LEVEL 
500 .000 
677.570 
644 .860 
177 .570 
ACTIVITY LEVEL 
369.206 
369.206 
.000 
1 .121 

16.449 
35 .000 
33.692 



UPPER LIMIT 



2000 .000 
50 ,000 
60 .000 
35 .000 
35 .000 

400 .000 
.00000600- 



TOLERANCES 05 05 05 05 05 



SIMPLEX MULT. 

.307- 
9.82? 

.075- 



SOL. VALUE 
369.206 
369.206- 
2000 .000 
48.879 
60 .000 
18.551 
.000 
1533.692 
300 .000 



LOWER LIMIT 
.000 
.000 



1500 .000 
300 .000 



ROW ERROR 
.00000140- 
.00000140 
.00000600- 

.00000015- 

.00000084- 

.00000005- 

.00000002 

.00000537- 

.00000090- 
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INPUT o C03600ALL0YC 

INPUT. 

ROWDSK 

COLDSK 

MTXDSK 

FSTDSK 

BASDSK 

ASSIGN0935 -1 
M I N . . . .ALOY1 VALUE 
LP1621 TO INVERT 



SEGMENT 5 



001 



NOISE DIGIT 
MANTISSA LENGTH 09 
SUBR SET 02 
LP1621 TO DUAL 
NOISE DIGIT 
MANTISSA LENGTH 09 
SUBR SET 02 
ITER VARBL. EXIT 

0001 BIN** 

0002 AL 
FEASIBLE 
OPTIMUM 

0UTPUTD04400ALL0YD 
NOISE DIGIT 
MANTISSA LENGTH 09 
SUBR SET 02 
OUTPUT 
BASIS. 
VARBLS 



VARBL. 

BIN5 

CU 



ENTR 



CNT 
005 
001 



OBJ. 



FUNCTION 

316.236 

333.611 



TYPE NAME 
FALUM 
FSILCON 
D0.D/JD04400ALL0YD 
NOISE DIGIT 
MANTISSA LENGTH 09 
SUBR SET 02 
VBLS TYPE NAME 
ROWS TYPE NAME 
COST .RDO 4400 A LLOYD 
NOISE DIGIT 
MANTISSA LENGTH 09 
SUBR SET 02 

COST.R NAME CURRENT 
ALUM 
SILCON 



MANTISSA 09 
ACTIVITY LEVEL 
416.612 
187.993 



TOLERANCES 05 03 03 03 



OVTPl/T 

03 



fePORT 



oo, D/r report 



CURRENT COST 
I NCR B VALUE 



REDUCED COST 
DECR B VALUE 



BASIS VALUE 



COST 
.210 
.380 



HIGHEST 



COST 
.230 
.480 



HI-VAR 

BIN3 

BIN3 



LO-VAR 

AL 

SI 



LOWEST 



COST 
.158 
.174- 
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GETOFF 
GETOFF 
ERASE? 

MAP...D02000ALLOYA n ^ ' 

MAP...S02000/02036 ALLOYA. 009ROW,0007COL.MANSA- 09-C0MPUTE-000086- SECTORS . 
MAP...D02800ALLOYB 

MAP...S02800/02813 ALLOYB. 009ROW/)007COL.MANSA-09- COMPUTE-000086- SECTORS . 
MAP...D03600ALLOYC 

MAP...S03600/03605 ALLOYC. 009ROW,0007COL,MANSA- 10- COMPUTE-000086-SECTORS • 
MAP...DOM*OOALLOYD 

MAP...S04400/0Mf06 ALLOYO. 000R0W,0002C0L 
ENDJOB 

END OF JOB £ NB of SAnfLE 
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1620-1311 LP SYSTEM ORGANIZATION AND CONTROL 



PROGRAM DECK 

The 1620-1311 Linear Programming System is supplied to the user in the form of a self- 
loading card deck, which is to be loaded onto disk under control of the 1620 Monitor System. 
To load the program, place the deck in the card reader, press RESET on the console, and 
then press LOAD on the card reader. 

The program deck consists of a series of object decks for the individual routines included in 
the 1620-1311 LP System. The first three cards of each object deck are: 

Monitor JOB 5 card 
Monitor **DUPcard 
Monitor *DLOAD card 

See Bibliography Reference 12, for Monitor card format. 

These cards call in Monitor, which will load the individual routine on disk and place the 
appropriate entries into the Monitor DIM Table and Equivalence Table. The entire 1620-1311 
LP System program deck is arranged as follows: 

Monitor cold start card 

+ + JOB 5 card v These four cards 

4=4= DUP card ( MUST appear here. 

*DLABL 11111 card ' 



Object 


deck, 


INPUT. 


Object 


deck, 


ROW. ID 


Object 


deck, 


COL. ID 


Object 


deck, 


MATRIX 


Object 


deck, 


FIRSTB 


Object 


deck, 


BASIS. 


Object 


deck, 


INPUTA 


Object 


deck, 


INPUTB 


Object 


deck, 


INPUTC 


Object 


deck, 


INPUTD 


Object 


deck, 


MTXDSK 


Object 


deck, 


REVISE 


Object 


deck, 


REVROW 


Object 


deck, 


REVCOL 


Object 


deck, 


REVMAT 


Object 


deck, 


REVSLK 


Object 


deck, 


REVFST 


Object 


deck, 


REVBAS 


Object 


deck, 


SAVE .B 


Object 


deck, 


OUTPUT 


Object 


deck, 


DO.D/J 


Object 


deck, 


GETOFF 


Object 


deck, 


COST.R 


Object 


deck, 


CHECK. 


Object 


deck, 


INVRT1 


Object 


deck, 


INV002 


Object 


deck, 


INV003 


Object 


deck, 


INV004 
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Object deck, INV005 
Object deck, IVLAST 
Object deck, NEWRHS 
Object deck, LP1620 
Object deck, LP1621 
Object deck, 1DUAL 
Monitor # * + * card 

SYSTEM ORGANIZATION 

All the routines included in the 1620-1311 LP System reside on the disk; they are called into 
memory individually as required. The sequence in which routines are called in and executed 
depends completely on the control cards which form the input deck to an LP run; these cards 
are called agendum cards. 

The executive program, called LP1621, reads an agendum card and calls in the required 
routine or series of routines. For example, when the DO.D/J agendum card is used, the 
DO.D/J routine is executed, and LP1621 reads in the next agendum cards. When MIN. . . 
(or MAX. . .) is used, however, LP1621 initiates an optimizing loop consisting of INVERT, 
DUAL- PRIMAL, CHECK. The loop continues until optimum solution is reached (unless the 
problem goes infeasible or unbounded, or a major error occurs); it is not until optimum 
solution is reached that another agendum card is read in. 

The agendum cards are always entered through the card reader. The data cards may be 
entered through the card reader, in which case the INPUT, agendum card will read INPUT. C 
in columns 1-7. However, the data used in a particular run may already be on the disk, in 
which case an INPUT. D agendum card is used. When input data is on disk, selected values 
may be changed through use of the REVISE agendum card without requiring the entire data file 
to be re-entered from the card reader. This means that the same basic data files stored on 
disk may be used in a number of LP runs. 

INPUT/OUTPUT ROUTINES 

All the programs of the 1620-1311 LP System use the GET and PUT statements of the 1620 
Monitor System for input/output. Thus the Monitor System I/O routines are in memory at 
all times, regardless of which LP System routine is being executed. All I/O messages are 
Monitor messages. 



OPERATING INSTRUCTIONS 



LOADING THE SYSTEM 

Before an LP problem can be run, the 1620-1311 Linear Programming System must be on the 
disk; once the system has been loaded, it need not be reloaded before each LP run. 

To load the system onto disk: 

• Load IBM 1620-1311 Monitor 1 on the disk 

• Place the LP actual deck in the card reader 

• Set all console sense switches OFF 
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• Press RESET on the console 



• Press LOAD on the card reader 



When the LP System has been loaded, the message END OF JOB will be typed out. 
RUNNING AN LP PROBLEM 

1. Place the LP deck in the card reader. The deck should be arranged as follows: 



When LP deck is first in the 
stack of Monitor jobs 

Cold start card 

+ 4= JOB 5 card 

+ +XEQ LP1620 card 



ENDJOB card 
PA US (optional) 



When LP deck is stacked 
behind another Monitor job 

**JOB 5 card 

+ +XEQ LP1620 card 



ENDJOB card 

PAUS (optional) 
+ 4= JOB 5 card 



(next Monitor job) 



These control cards (except PAUS) are required, regardless of whether the data is in 
the card reader or on disk. When the 1620-1311 LP System reads the ENDJOB card, it 
will return control to the 1620 Monitor System, which then types END OF JOB. If the 

PAUS card is used, there will be a halt; if not, 1620 Monitor will then try to read 
another card from the card reader. (See Bibliography Reference 12, for Monitor card 
format.) 



2. Set the console sense switches. 



SWITCH 


IF ON 


IF OFF 


1 


Each iteration will be logged on the 
typewriter. This switch setting may 
be changed as often as desired dur- 
ing problem run. Iterations will also 
be logged per ASSIGN control. 


Iterations will be logged 
per ASSIGN control. 


2 


Used to interrupt optimization pro- 
cedure. (See previous description of 
GETOFF procedure.) 


Normal operation. 


3 


Will not halt or type MINOR ERROR 
if minor errors occur. Will supress 
***** NOT IN COL FILE and ***** 
ROW NOT IN FILE messages which 
occur in INPUT when selecting 
columns and rows respectively 
from MATRK and RHS files. 


Will type MINOR ERROR 
and halt if minor error 
occurs. The position of 
switch 4 controls the 
succeeding action. 
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SWITCH 




IF ON 


IF OFF 


4 


Will 


type MAJOR ERROR and halt if 


Will type MAJOR ERROE 




major error occurs. If START is 


and halt if major error 




then 


pushed, will read in next con- 


occurs. If START is then 




trol 


card and attempt to continue. 


pushed, will execute 








END JOB. 



3. Set all console check switches to PROGRAM. 



4. Press RESET on console. 

5. Press LOAD on card read punch if LP deck is first in the stack of Monitor jobs. 

6. After a * + PAUS card, set console switches as desired and press START. 

Some fields on the Monitor control cards which are needed to describe the 1620 system 
used are listed below. (The Monitor control cards are fully described in Bibliography 
Reference 12.) 

+ + JOB 

Col. 1-5 

7 

8-11 



+ +XEQ 

Col. 1-5 
7-12 
31 - 32 



*DFINE 

Col. 1-6 

18 



+ + JOB 

5 must always be present 

blank if 1 disk system 

01 if 2 disk system 

012 if 3 disk system 

0123 if 4 disk system 

**XEQ 
LP1620 

03 if system has floating-point 

hardware 

blank otherwise 



*DFINE 

The number of disk storage drives on the system. 
Maybe one, two, three, or four. One drive is 
assumed if an *DFINE card is not present. See 
Reference 12 for a description of where the 
*DFINE card should be placed. 



108 



HALTS AND MESSAGES 



The following list contains all messages which will be typed out during a 1620-1311 LP System 
run, except for the following: 

As each agendum card is readmit is logged (on unit ASSIGNed for logging). 

All comments cards (* in column 1) in the agendum stream are logged. 

All disk I/O operations are under control of the Monitor I/O System, and any disk I/O 
messages are Monitor messages. If any disk I/O errors occur, a message is typed out 
for appropriate action. (See Reference 12.) 

The two program switches in the 1620-1311 LP System are the minor and major error 
switches. When any program finds an error condition, it types out an appropriate message, 
sets the program error switch, continues processing if possible, and then returns control to 
the executive program (LP1621). If the error is minor, the action of the executive program 
depends upon the setting of console switch 3. If ON, the error is ignored and the LP run 
continues without a halt. If OFF, the minor and major messages are typed and the program 
halts. If the error is major, the major error message is always typed and the program 
halts . 

The action that the user should take when an error halt occurs depends upon the cause of the 
error. On the following pages are listed all the error messages. The Error Restart pro- 
cedures referenced are described later under "Error and Interrupt Restart Procedures". 

There are several messages for which no restart procedure is given. These are due to pro- 
gram errors. If the user has modified the programs, he must correct his modifications. 

*MAX ERROR = xxxx. xxxxxxxx- CHECK. 

Information message only. Maximum row error generated during the preceding set of itera- 
tions. A CHECK and reinversion is automatically taken as specified by the inversion fre- 
quency. This message occurs at the end of CHECK, and yields the maximum error detected. 
This error is used to adjust the mantissa length and tolerances. 

****** COL NO LOWER BOUND REVMAT 

Minor error. ****** is the name of a column without a lower bound. It is not possible to 
revise the lower bound of a variable unless the variable is input with an explicit lower bound. 
The Revise Matrix element card contains a lower bound. If the card is mispunched use Error 
Restart B. If the lower bound is required use Error Restart C or D. 

****** COL NO UPPER BOUND REVMAT 

Minor error. Same as above except that this message is to the upper bound. 

****** COL NOT IN FILE REVCOL 

Minor error. ****** is a column name which is not in the COL. ED file. If the Revise 
column identification card is misspelled, use Error Restart B. If the column identification 
is to be added, use Error Restart C or D. 
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****** COL NOT IN MATRIX 



REVMAT 



Minor error. ****** is the name of a column which is not in the MATRIX file. The revise 
matrix element card identifies a column not in the matrix. If the column name is misspelled, 
use Error Restart B. If a column is to be added to the matrix, use Error Restart C or D. 



Minor error. ****** is the name of a variable (structural or logical) that has already been 
input. A second basis identification card names a previously read variable. The entry is 
correct if the first of the duplicate names specifies the correct basis type. In that case use 
Error Restart F. If the first of the duplicate basis cards specifies an incorrect basis type, 
use Error Restart G or A. 

****** DUP COL COL. ID 

Minor error. See above, except that reference is to column name. A second column identi- 
fication card names a previously read variable. The entry is correct if the first of the 
duplicate names specifies the correct basis type. In that case use Error Restart F. If the 
first of the duplicate basis cards specifies an incorrect column type, use Error Restart 
A or E. 

****** D UP COL NAME INPUT B 

Minor error. ****** is the name of two columns in the MATRIX file. Detected in the 
generation of a COL. ID file from the MATRIX file. This is usually due to: cards out of order 
in the MA.TRLX deck, a blank card in the MATRIX deck, two variables with the same name, 
two sets of coefficients for the same variable, or a misspelled variable name. Use Error 
Restart A. 

****** DUP NAME IN AJ FILE INPUT C 

Minor error. ****** is the name of two columns in the AJ VAR or MATRIX file. Detected 
in initialization of the MATRIX file from the COL. ID file. See preceding message for cause 
and action. 

****** DUP NAME IN AJ SLK FILE INPUT D 

Minor error. ****** is the name of two logical variables. Detected in initialization of 
AJ SLACK file from BASIS SLACK file. The program which generated the ROW. ID file, or 
the program which generated the logical variables from the ROW. ID file is in error. This 
error may occur due to user modification or replacement of input routines. The second 
logical variable will be classified as an omit (*) variable regardless of the row type. 

****** DUP NAME IN AJ VAR FILE INPUT D 

Minor error. ****** is the name of two columns in the MATRIX file. Detected in the 
initialization of the AJ VAR (MATRIX) file from the BASIS VARIABLE file. See ****** DUP 
COL NAME message for cause and action. 

****** DUP ROW ROW. ID 

Minor error. ****** is the name of a row that has already been input. A second row iden- 
tification card names a previously read row. The entry is correct if the first of the duplicate 
names specifies the correct row type. In that case use Error Restart F. If the first of the 
duplicates specifies an incorrect row type, use Error Restart G or A. 



****** DUP BASIS NAME 



BASIS. 
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****** G SLK NO UPPER BD 



OUTPUT 



Major error. The logical variable generated for row ****** has no upper bound. This 
variable is classified as a type G variable in the AJ SLACK file. The input or revise basis 
cards specified that the logical variable was a type G variable. This classification is valid 
only for equation rows and range constraints. Use Error Restart B if revised incorrectly. 
Use Error Restart A or G if input incorrectly. This error may also occur due to modifica- 
tion of optimization programs. 

****** G VAR NO UPPER BD OUTPUT 

Major error. The structural variable named ****** has been classified as a type G variable 
in the AJ VARIABLE or MATRIX file. This variable does not have an upper bound. The 
MATRIX file has been misused by an optimization program. This may be due to modification 
of an optimization program. 

nnnnnn* IN AJ SLK FILE INPUT D 

Minor error, nnnnnn is the name of a row classified as an OMIT (*) row. The logical 
variable for this row must not appear in the basis slack file . A basis identification card 
named the logical variable for an omit row. The row remains classified as an omit row and 
the basis classification is ignored. If that is the desired result, use Error Restart F. To 
enable this row to be used as a constraint, correct the row identification card and use Error 
Restart A or E . 

nnnnnn* IN A J VAR FILE INPUT D 

Minor error, nnnnnn is the name of a variable classified as an omit (*) in the MATRIX and 
COL. ID files. A basis identification card has been read for this variable. Omitted columns 
cannot be specified in the basis variable cards. The basis classification is ignored and the 
column remains omitted. If that is desired, use Error Restart F. If the variable is to be 
included in the problem, prepare a column identification card and use Error Restart A or E. 

****** INP ROW NOT IN MTRX ROW FILE MTXDSK 

Major error. ****** is the name of a row in the ROW. ED file for the problem which is not 
in the ROW. ID file of the matrix. A row identification card is mispunched. Use Error 
Restart A. 

****** IS not IN EQU TABLE LP1621 

Major error. The agendum name ****** col. 1-6 of the current agenda card does not cor- 
respond to any of the subroutines of LP1621 and this name is not in the IBM 1620 Monitor 
Equivalence (program name) table. The agendum card is misspelled (DObD/J instead of 
DO.D/J for example), or cards are out of order in the agenda deck, or there are duplicate 
EOF or ENDATA cards. Use Error Restart H. 

****** NO ENTRY IN RHS nnnnn REVFST 

Minor error. The input RHS nnnnn did not have an entry for row ******. A right-hand- side 
Revise element card has been read for this coordinate. If no entry is required, use Error 
Restart F. If this entry is required, use Error Restart C or D. li the card is misspelled, 
use Error Restart B. 
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****** NO ROW ENTRY IN COL nnnnnn 



REVMAT 



Minor error. The input column nnnnnn did not have an entry for row ******. A Revise 
matrix element card specifies this coordinate. If no entry is required, use Error Restart F. 
If the card is misspelled, use Error Restart B. If this entry is required, use Error Restart 
C or D. 

****** NOT IN BASIS YYY FILE REVBAS 

Minor error. ****** i s the name of a variable which is not in a basis file. YYY = SLK, 
basis slack file. YYY = VAR, basis variable file. A Revise basis identification card names 
a variable not in the basis file. If the card is misspelled, use Error Restart B. If this vari- 
able must be reclassified from the initial status (W - structural variables; F - logical vari- 
ables), use Error Restart C or D. 

****** NOT IN COL FILE INPUT C 

Minor error. ****** is the name of a variable which is in the MATRIX file but not in the 
column identification file. This message is suppressed if console switch 3 is on during 
INPUT. The input COL. ID file did not include this variable, which is in the MATRIX. This 
column is classed omit (*) . If that is desired, use Error Restart F. If the variable is to be 
included in the problem, use Error Restart A. 

***** RHS NOT IN B FILE REVFST 

Minor error. There is no RHS named ***** in the B (right-hand-side) file. The right-hand 
side specified by the FIRST. B or NEXT.B Revise indicator card does not exist. The right- 
hand-side element cards following this card are ignored. If the FIRST. B or NEXT.B card is 
misspelled, use Error Restart B. If a right-hand side is to be added, use Error Restart 
C or D. 

****** ROW NOT IN BETA 1 OUTPUT 

Major error. ****** is the name of a variable classified as a type F in the MATRIX or AJ 
VARIABLE FILE (structural), or the AJ SLACK FILE (logical). This name does not appear 
in the solution list BETA 1. This error is probably due to file destruction or misuse by an 
optimization or output program: 1 DUAL, INVERT, CHECK, OUTPUT, DO.D/J, or COST. R. 
It may be due to user program modification or replacement. 

****** ROW NOT IN FILE FIRST B 

MATRIX 

Minor error. Row named ****** is not in the ROW. ID file and is specified in a RHS row 
coordinate or on a MATRIX element card. This message is suppressed if console switch 3 
is on during INPUT. (Note: If switch 3 is off and there is a blank card in the MATRLX cards, 
the message will appear to be "bbbbbb ROW NOT IN FILE".) The input ROW. ID file did not 
include this variable. This card is bypassed and no entry occurs. If that is desired, use 
Error Restart F. If the row is desired, use Error Restart A. 

****** ROW NOT IN FILE REVROW 

Minor error. The input ROW. ID file did not include a row named ******. if the REVISE 
row identification card is misspelled, use Error Restart B. If a row is to be added, use 
Error Restart C or D. 
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****** ROW NOT IN ROW* ID FILE 



REVFST 
REVMAT 



Minor error. Row ****** j s no t m the ROW, ID file. The Revise matrix element or right- 
hand-side element card specifies a row name not in the ROW. ID file. If the card is mis- 
spelled, use Error Restart B. If a row is to be added to the MATRDX (and RHS), use Error 
Restart C or D. 

****** SLK NOT IN ROW. ID FILE SEQUENCE REVSLK 

Minor error. The logical variable generated for row ****** is not in the same order as the 
ROW. ID file. The AJ SLACK file has been generated incorrectly. This may be due to modi- 
fication of input programs. 

A AND/OR SLACK FILE NOT DEFINED CHECK. 

SAVE.p 

DO.D/J, NEWRHS 

Major error. The A J VAR (MATRIX) or the AJ SLACK file does not exist. This file(s) is 
required for processing the agendum. The agenda is out of order. Use Error Restart H. 

AGENDA NAME NOT EQUAL CNTL. DO.D/J 
FILE NAME 

Major error. The name on the DO.D/J agendum card does not match the specified control 
file name. The agendum card is mispunched, either the problem name and/or the control 
file location is incorrect. Use Error Restart H. 

BASIS FILE FULL BASIS. 

Minor error. The number of variables in the BASIS variable file exceeds the number of 
columns in the matrix, or the number of logical variables exceeds the number of rows in the 
ROW. ID file. Remove the unnecessary basis cards. Use Error Restart A. 

BETA FILE NOT DEFINED CHECK. 

Major error. The CHECK, agendum is called for before the Beta file (names and activity 
levels of intermediate-level basis variables) has been formed. The agendum cards are out 
of order. Use Error Restart H. 

CANNOT FIND MATCHING RHS CHECK. 

Major error. The RHS name specified in columns 8-12 of the INVERT, MIN, or MAX 
agendum card is not the name of an R-H-S in the RHS file. Correct the agendum card and use 
Error Restart H. 

CANT FIND NAMED RHS NEWRHS 

Major error. The RHS name specified in columns 8-12 of the INVERT, MIN, or MAX 
agendum card is not the name of an R-H-S in the RHS file. Correct the agendum card and use 
Error Restart H. 
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CANT CONTINUE. INV FILE GOING BEYOND 
USER AREA. 



INVERT 



Major error. The upper limit ASSIGNed by the user has been exceeded. 

1. Preferred method. Repuneh the ASSIGN card with a higher upper limit or a lower 
common computation area sector address. Use Error Restart H. 

2. Specify a lower starting sector address on the INPUT card. Use Error Restart A. 
Depending upon the size of the problem, this method may or may not yield a successful 
run. 

CANT CONTINUE. INV FILE GOING INTO INVERT 
DIM TABLE. 

Major error. The inverse file is about to exceed the disk storage area allowed. If the pro- 
gram had not halted, the Monitor DIM table required by the system would be destroyed. 

1. Preferred method. Prepare an ASSIGN card with 

a. Sector address — DIM 

b. Sector address — UL 

c. Common computation area sector address specifying an area above the DIM table 
large enough for the inverse file. 

Use Error Restart H. 

2. Specify a lower starting sector address on the INPUT card. Use Error Restart A. 
Depending upon the size of the problem, this method may or may not yield a successful 
run. 

CANT FIND OBJ FCT IN ROWID LP1621 

NEWRHS 

Major error. The objective function specified on the MIN or MAX card is not in the ROW. ID 
file. The agendum card is mispunched. Use Error Restart H. 

CDP ERR XXXXX 

Monitor error message; program halts. May be due to invalid disk sector addresses 
(multiple drives). The disks on drive n (as specified on the Monitor JOB card) may have 
invalid sector addresses. Usually occurs when a disk defined for use on drive (sector 
addresses 00000-19999) is placed on drive n, n ^ 0. May also be due to program error. 

CNTRL CARD NAME DOES NOT MATCH INPUT A 

PRB CTL FILE 

Major error. The input control card (ROW. ID, BASIS. , MATRIX, FIRST. D, COL. ID) 
problem name does not match the specified control file name. The control card is mis- 
punched. Use Error Restart A. 

COL FILE BEING MADE IS OVER LIMIT INPUT B 

Minor error. The number of columns in the MATRIX file exceeds capacity. The remaining 
columns are not entered in the COL. ID file and will subsequently be classed omit (*) in the 
MATRIX file. The problem is too large for the machine size. Use Error Restart A. 

COL FILE OVER LIMIT COL. ID 

Minor error. The number of column names in the column identification file exceeds capacity. 
The remaining column names are bypassed. The problem is too large for the machine size. 
Use Error Restart A. 



114 



DATA ERROR CANNOT OPERATE WITH 
MATRIX FILE 



1 DUAL, COST.R 



Major error. Vector product exponent exceeds 99. Program cannot operate with inverse 
and/or data files. The INVERSE file has been destroyed by a disk operation, or the MATRIX 
file has been incorrectly formed or has subsequently been destroyed. If the INVERSE is 
destroyed, use Error Restart I. If the MATRIX files have been destroyed, use Error Restart 
C or D. Incorrect MATRIX file formation may be due to modification or replacement of input 
routines. 

DATA ERROR, TYPE F NOT IN BETA COST.R 

Major error. A variable is classified as a type F variable in the AJ VAR (MATRIX) file. 
This variable is not in the Beta file (list of type F variables) . The error is probably due to 
file destruction or misuse by an optimization or output program: 1DUAL, INVERT, CHECK., 
OUTPUT, DO.D/J, or COST.R. It may be due to user modification or replacement. 

DSK ERR 0607l6l7363738 MONITOR 

Monitor error message. Machine halts. Error is probably due to incorrect sequence of 
agendum cards. Use Error Restart A. 

DUP B FILE INPUT A 

Major error. A FIRST. B input control card is read when a B or RHS file has already been 
input. The input deck is invalid. Use Error Restart A. 

DUP BASIS FILE INPUT A 

Major error. A BASIS, input control card is read when a BASIS, file has already been input. 
The input deck is invalid. Use Error Restart A. 

DUP COL FILE INPUT A 

Major error. A COL. ID input control card is read when a COL. ID file has already been 
input. The input deck is invalid. Use Error Restart A. 

DUP MATRIX FILE INPUT A 

Major error. A MATRIX input control card is read when a MATRIX file has already been 
input. The input deck is invalid. Use Error Restart A. 

DUP ROW FILE INPUT A 

Major error. A ROW. ID input control card is read when a ROW. ID file has already been 
input. The input deck is invalid. Use Error Restart A. 

DUPL ENTRY FIRSTB; MATRIX 

Minor error. A duplicate data entry for RHS or MATRIX element card is found. The 
duplicate card is punched out. The first element is stored in the disk MATRIX, the second 
is ignored. If the first element is correct, use Error Restart F. If the first element is not 
correct, use Error Restart A or E. 
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ENT ERROR 06071617363738 



MONITOR 



Monitor error message. Computation continues. A preceding I/O operation is erroneous 
and caused an error indicator to be turned on. Error may be due to machine malfunction or 
program error. If machine malfunction, use Error Restart K. 

ERROR IN INPUT CONTROL CARD INPUT. 

Major error. Card column 7 of the INPUT, agendum card is not C or D. Use Error 
Restart A . 

ETA 11 FILE NOT DEFINED CHECK. 

NEWRHS 

Major error. The ETA11 file does not exist. The agendum cards are out of order or a pro- 
gram error resulted in file destruction. Use Error Restart H if the agendum cards are out 
of order. 

FEASIBLE 1 DUAL 

Information message; no halt occurs. The current solution is now feasible. 

FILE MISSING xxxxxx COST.R 

1DUAL 

Major error, xxxxxx File is required for processing and has not been input or formed. (Can 
be Beta, MATRIX, AJ SLACK, ETA 11, or INVERSE file.) Agendum cards are out of order, 
no (or inadequate) input data exists, or COST.Ranging precedes optimization. Use Error 
Restart H. 

FILE OVER LIMIT BASIS. 

Minor error. The number of structural (logical) variables in the BASIS VARIABLES 
(SLACK) file exceeds machine-file capacity. The problem is too large for machine-size 
capacity. Use Error Restart A. 

G W/O UPBND INVERT 

Major error. A variable has been classified as a type G (in basis at upper bound), but the 
variable has no upper bound. Use a GETOFF agendum to determine the variables and types, 
correct the types and use Error Restart G. 

ILLEGAL MANTISSA LENGTH - xx INVERT 

Major error. The mantissa length specified on the ASSIGN card is not in the allowable 
range (5-18). Correct the ASSIGN card. Use Error Restart H. 

INFEASIBLE VARIABLE xxxxxx 1 DUAL 

Minor error/major error. Variable xxxxxx is infeasible and cannot be removed from the 
basis. This is a minor error the first time it occurs; system will automatically enter 
CHECK, and then INVERT to see whether this was caused by digital difficulties. If the 
situation recurs in 1 DUAL, it is a major error and the system halts because of infeasible 
problem formulation. Use Error Restart A. 
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INPUT CNTRL CARD NAME DOES NOT 



INPUT. 



MATCH PRB CTL FILE 

Major error. The INPUT, agendum card problem name does not match the specified control 
file name. The agendum card is mispunched. Use Error Restart A. 



Major error. The number of rows in the ROW. ID file for the problem is greater than the 
number of rows in the ROW. ID file of the matrix. If it is desirable to add rows to the matrix, 
use Error Restart C or D. If adding rows to the matrix is not desired, use Error Restart A. 



Major error. The inverse file is about to exceed the disk storage area allowed. If the pro- 
gram had not halted, the Monitor DIM table required by the system would have been destroyed. 

1. Preferred method. Prepare an ASSIGN card with 

a. Sector address — DIM 

b. Sector address — UL 

c. Common computation area sector address specifying an area above the DIM table 
large enough for the inverse file. 

Use Error Restart H. 

2. Specify a lower starting sector address on the INPUT card. Use Error Restart A. 
Depending upon the size of the problem, this method may or may not yield a successful 
run. 

INVERSE FILE EXCEEDS ULPRB. 1DUAL 

Major error. The upper limit ASSIGNed by the user has been exceeded. 

1. Preferred method. Repunch the ASSIGN card with a higher upper limit or a lower 
common computation area sector address. Use Error Restart H. 

2. Specify a lower starting sector address on the INPUT card. Use Error Restart A. 
Depending upon the size of the problem, this method may or may not yield a successful 
run. 

LP 1621 AUTO ASSIGN NEW MANSA IS XX LP 1621 

Information message; no halt occurs. The problem mantissa is now xx. 

LP 1621 NEW OBJ FCT LP1621 

Information message; no halt occurs. Optimization subroutine to change objective function 
row. 

LP 1621 TO RRRRRR LP 1621 

Information message; no halt occurs. The next routine in the optimization sequence is 
RRRRRR. 

RRRRRR = DUAL for iterations. 
RRRRRR = PRIMAL for iterations. 
RRRRRR = CHECK, to check accuracy. 
RRRRRR = INVERT to (re-)invert basis. 
RRRRRR = NEWRHS to change R-H-S. 

RRRRRR = INVERT CLEANUP when the basis to the INVERT program is singular. 



INPUT ROW FILE GREATER THAN MATRLX 
ROW FILE 



MTXDSK 



INVERSE FILE ENTERS DIM. 



ID UAL 
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MA. J. ERROR. SW 4 ON TO CONTINUE OR 
OFF TO ENDJOB 



LP 1621 



Information message; program halts. Major error has occurred. See preceding logical 
error messages to determine the cause of the major error. Consult message list to clarify 
the cause of the error and to determine the corrective action required. 

MAP REFERENCE NAME DOES NOT AGREE WITH LP 1621 

CONTROL FILE REFERENCE 

Major error. The MAP. . . agenda card problem name does not match the specified control 
file name. The agendum card is mispunched. Use Error Restart A. 

MATRIXD AND PRB CTL FILE NAMES DO MTXDSK 
NOT MATCH 

Major error. The name specified in columns 13-18 of the MATRIXD indicator card does not 
match the name in the sector given by columns 8-12 of the MATRIXD card. If the indicator 
card is mispunched, use Error Restart A. If the indicator card is not mispunched, the 
referenced file has probably been destroyed. A matrix must be entered on cards. Use 
Error Restart A. 

MINOR ERROR WITH SW 3 OFF MA J. ERROR. SW 4 LP 1621 

ON TO CONTINUE OR OFF TO ENDJOB 

Information message; program halts. A minor error has occurred. This message is sup- 
pressed and no halt occurs if console switch 3 is ON. See preceding logged error messages 
to determine cause of error. Consult message list to clarify the cause of the error and to 
determine the corrective action required. 

MOD ERR xxxxx MONITOR 

Monitor error message; machine halts. Invalid disk address due to failure to define, by a 
Monitor *DFINE card, the number of disk storage drives to MONITOR, or invalid JOB 
card, or invalid sector addresses on drive 1, 2, or 3. 

NO AREA SET UP FOR SAVE.p SAVE.p 

Major error. No area has been set aside for the basis file. Agendum cards are out of 
order, SAVE.p precedes optimization. Use Error Restart H. 

NO xxxxxx FILE REVISE 

Major error. No xxxxxx file exists but revise control card xxxxxx has been read. The data 
cards are bypassed. The file maybe BASIS., COL. ID, FIRST. B (R-H-S), MATRIX, ROW. ID. 
This is an attempt to revise a file which has not been INPUT or SAVEd. If the file is not 
required, use Error Restart F. If the file is required, use Error Restart C or D . 

NO DATA EXISTS FOR MIN, MAX, OR LP 1621 

INVERT 

Major error. A file or files required for processing have not been INPUT. Agenda deck or 
input deck out of order or inadequate. Use Error Restart A. 
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NO DATA FILES - CANNOT OUTPUT 



OUTPUT 



Major error. Data files are nonexistent or invalid. The agendum cards are probably out of 
order. Use Error Restart A. 

NO IND CARD BASIS. 

INPUT A 

Major error. The card which follows the INPUT, agendum or the BASIS, data indicator card 
is not a data indicator card. The data cards are bypassed until an indicator card is read. 
The INPUT, agendum must be followed by either a ROW. ID or a COL. ID data indicator card. 
The BASIS, data indicator card must be followed by a VARBLS data indicator card. Use 
Error Restart A. 

NO INPUT CONTROL CARD INPUT A 

Major error. The indicator card read (cols 1-3 not blank, col 1 not *) is not an INPUT, data 
indicator card. The card in question is punched out. Error is due to mispunched indicator 
card or missing ENDATA or EOF card. Use Error Restart A. 

NO INVERSE EXISTS. ONLY VARBLS AND SAVE.p 
SLKS SAVED 

Informative message; no halt occurs. The agendum card specified SAVE. I (save inverse), but 
no inverse file exists. The BASIS VARIABLE and BASIS SLACK files are formed and saved. 

NO MATRIX — CANNOT ACCEPT BASIS INPUTA 

Major error. The BASIS, indicator card has been read, but the MATRIX file has not been 
read. The data cards are bypassed. The input deck does not contain a MATRIX, or the 
MATRIX follows the BASIS. Use Error Restart A. 

NO MATRIX ~ CANNOT FORM B FILE INPUT A 

Major error. The FIRST. B indicator card has been read, but the MATRIX file has not been 
read. The data cards are bypassed. The input deck does not contain a MATRIX, or the 
MATRLX follows the BASIS. Use Error Restart A. 

NO MATRIX FILE TO REV B FILE REVISE 

Major error. The MATRIX file does not exist. A B or R-H-S file is valid only if a 
MATRIX file exists. The data cards are bypassed. The data file for this problem is invalid 
and must be corrected and re-INPUT. before use. Use Error Restart C or D. May also be 
due to modification or replacement of input programs. 

NO OBJ. FUNCTION ROW IN FILE INVERT 

Major error. There is no row classified as an objective function in the ROW. ID file. The 
objective function row identification card is either missing or mispunched. If it is missing, 
use Error Restart A. If it is misclassified, use Error Restart A or E. 
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NO REV IND CARD 



REVISE 



Major error. The card(s) which follows the REVISE agendum card is not a data indicator 
card (or a comments card * in col. 1). These cards are bypassed until an indicator card is 
read. A data indicator card must follow the REVISE agendum card. Cards are missing or 
out of order. Use Error Restart B. 

NO ROOM FOR PROB. ON DISK BASIS., COL. ID, 

FIRSTB, INPUT., 
INPUTB, MATRIX, 
MTXDSK, ROW. ID, 
SAVE.p 

Major error. The files created by the program at the sector address specified in the INPUT, 
or SAVE.p agendum card (col. 8-12) exceed the upper limit for the problem. An ASSIGN card 
can be used to modify this limit. A mispunched sector address in the agendum card or a 
missing ASSIGN card causes this error. Use Error Restart H. 

NO ROW FILE-CANNOT FORM MATRIX INPUT A 

Major error. The MATRIX data indicator is read before formation of the ROW. ID file. The 
ROW. ID file is missing or does not precede the MATRIX. Use Error Restart A. 

NO ROW FILE TO REV MATRIX REVISE. 

Major error. The current problem file does not contain a ROW. ID file. The revision is 
being made to the wrong file (incorrect INPUT agenda), the files have been destroyed by sub- 
sequent disk usage, or by trying to add a MATRIX to a file. 

NO SLACKS IND CARD BASIS. , REVBAS 

Major error. The BASIS, deck does not have a SLACKS indicator card. The SLACKS indi- 
cator card is misplaced or missing. Use Error Restart A if the error occurred during input 
of a BASIS. Use Error Restart B if the error occurred during RE Vision of a BASis. 

NO SPACE AVAILABLE ABOVE DIM LP1620 

Informative message only; not an error. The area above the DIMension table (cylinder 24) 
to the monitor program area (cylinder 80) is Monitor-protected. This area contains pro- 
grams and data under Monitor control and must not be used for LP1620 data storage or com- 
putation. The only area available for LP data use is sector 01805 to sector 04799, if the 
system has one IBM 1311 disk device. 

NO UPPER BOUND FOR TYPE G CHECK. 

Major error. A variable is classified as a type G variable in the AJ VAR (MATRIX) file. 
This variable has no upper bound. Error is probably due to file destruction or misuse by an 
optimization or output program: ID UAL, INVERT, CHECK., OUTPUT, DO.D/J, or 
COST.R. It may be due to user program modification or replacement. 

NO VARBLS IND CARD BASIS. , REVBAS 

Major error. The BASIS, deck does not have a VARBLS indicator card. The VARBLS data 
indicator card is misplaced or missing. Use Error Restart A if the error occurred during 
input of a basis. Use Error Restart B if the error occurred during revision of a basis. 
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NO VARBLS OR SLACKS CARD 



BASIS. 



Major error. The first data indicator card following the BASIS, card in the BASIS deck is 
not VARBLS and is not SLACKS. The input deck is out of order or indicator cards are 
missing. Use Error Restart A. 

NO VBLS LEFT IN BETA 1 CHECK. 

Major error. Program error in CHECK, program. Error is probably due to user 
modification. 

NOT A REVISE CONTROL CARD REVISE 

Major error. The indicator card (column 1-3 not blank and column 1 not *) is not a data 
indicator card: ROW. ID, COL.LD, BASIS., MATRLX, FIRST. B, EOF, or ENDATA. The 
indicator card is punched out. The REVISE deck is out of order or cards are missing. Use 
Error Restart B. 

NOT ENOUGH DISK FOR BETA LP1621 

Major error. The disk area remaining is inadequate to allow room for a BETA file. This 
area is part of the compute area. The starting sector address for the compute area can be 
specified on an ASSIGN card. Prepare an ASSIGN card and use Error Restart H. 

NOT ENOUGH ROOM FOR SAVED VBL AND LP1621 
SLACK NAMES NXTDSK = xxxxx 

Major error. The disk area remaining is inadequate to allow room for a BASIS file. The 
BASIS area is space for a BASIS SLACK file and a BASIS VARIABLE file. This area is part 
of the input area. It is set aside prior to optimization (MIN or MAX). The ASSIGN card 
(see "1620-1311 LP Agenda") can be used to increase the available area. A lower sector 
address on the INPUT agendum card, if possible, will increase the available disk area. Use 
Error Restart A if the INPUT agendum is changed. Use Error Restart H if the ASSIGN 
agendum is used. 

OBJ. FUNCTION OR OMIT ROW HAS CLASS G INVERT 
IN SLACK FILE 

Major error. The BASIS SLACK file classifies an objective function or an omit row as a 
class G variable. The ROW. ID file or the BASIS SLACK file should be revised to correct 
the discrepancy. Use Error Restart E or G respectively. 

OPTIMUM 1 DUAL 

Information message only. The current basis is optimum. 

OUTPUT AGDM CARD NAME DOES NOT OUTPUT 
MATCH PRB CTL FILE 

Major error. The name specified on columns 13-18 of the OUTPUT agendum does not match 
the name in sector given by columns 8-12 of the OUTPUT agendum. If the name or sector 
address is mispunched, use Error Restart H. If the name and sector address are correct, the 
referenced file has probably been destroyed. Prepare an OUTPUT agendum card with blanks 
in column 7-18. Use Error Restart H and manually select output. Examine INPUT and SAVE 
agenda and use MAP. . . agenda to determine which agenda caused destruction of the referenced 
file. 
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PROB. NAME ON CHECK CARD DOES NOT 
MATCH CTL. FILE 



CHECK. 



Major error. Same as preceding except that agenda is CHECK. 

PROBLEM COMMON CONFLICTS WITH DISK INPUT. 
LIMITS 

Major error. The sector address specified in columns 8-12 of the INPUT agendum card is 
not in the available disk area. An ASSIGN card should precede the INPUT card, making 
available the area specified by the sector address. Use Error Restart A. 

PROC. ERROR EXCEEDS TOL. FOR VARIABLE 1DUAL 

Information message only. Processing error for variable ****** selected to enter the basis 
in this iteration exceeds the error tolerance. The system will increase the mantissa length 
(if possible) and change tolerances, INVERT, and continue optimization. 

RECORD MARK ON AGENDUM CARD UNINTER PRE TABLE LP 1621 

AGENDUM CARD 

Major error. The agendum card has record marks in columns 1-6. Use Error Restart H. 

ROW. ID FILE NOT DEFINED NEWRHS 

Major error. The ROW. ID file does not exist. The agendum cards are out of order or a 
program error resulted in file destruction. Use Error Restart H if the agendum cards are 
out of order. 

RHS FILE NOT DEFINED CHECK., NEWRHS 

Major error. The R-H-S file does not exist. The agendum cards are out of order or a pro- 
gram error resulted in file destruction. Use Error Restart H if the agendum cards are out 
of order. 

RHS NAME NOT FOUND IN FILE INVERT 

Major error. The RHS name specified in columns 8-12 of the INVERT, MIN, or MAX 
agendum card is not the name of a R-H-S in the RHS file. Correct the agendum card and 
use Error Restart H. 

ROW FILE OVER LIMIT ROW. ID 

Minor error. The number of rows in the INPUT ROW. ID file exceeds machine capacity. 
The excess rows are bypassed. Use Error Restart A. 

ROW = * INVERT 

NEWRHS 

Major error. There are two entries in the right-hand side for a row defined as an equation 
in the ROW. ID file. Use Error Restart E to either change the row type or to delete one (or 
both) entries in the right-hand side. 
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ROW NOT IN FILE 



FIRSTB, MATRIX, 
REVROW 



Minor error . See error message ****** ROW NOT IN FILE . Probably a blank or mis- 
punched card in the R-H-S or MATRIX element cards. Use Error Restart A. 

S*****CONTROL FILE NOT YYYYYY COST.R 

Major error. The name on the COST.R agendum card does not match the specified control 
file name. The agendum card is mispunched. Either the problem name and/or the control 
file location are incorrect. Use Error Restart H. 

SLACK HAS INVALID CLASS INVERT 

Major error. A basis identification card contains an invalid slack class. Turn console 
switches 2 and 4 on, remove all cards from read hopper, place GETOFF agendum card in the 
read hopper, press START, turn console switches 2 and 4 off, list GETOFF output cards, 
correct invalid slack class(es). Use Error Restart G. 

THIS IS NOT A COST ROW INVERT 

Major error. The objective function name xxxxxx, specified on columns 13-18 of the INVERT, 
MAX, or MIN card, is not a row classified as an objective function in the ROW. ID file. If the 
row type is incorrect in the ROW. ID file, use Error Restart A or E. If the row name is 
incorrect in the INVERT, MAX, or MIN agendum card, use Error Restart H. 

THIS ROW NOT FOUND IN FILE INVERT 

Major error. The objective function name xxxxx, specified on columns 13-18 of the INVERT, 
MAX, or MIN card, is not in the ROW. ID file. The agendum card is probably mispunched. 
Use Error Restart H. 

TRIED TO FLOAT INVALID NUMBER FIRSTB, MATRIX, 

REVFST, REVMAT 

Minor error. The element card contains an invalid data field. The data is bypassed and the 
erroneous card is punched out. The card is probably mispunched. Use Error Restart A if 
the error occurred in data INPUT. Use Error Restart B if the error occurred in data 
REVISION. 

TWO CONSECUTIVE FAILURES IN OPTIMIZING LP1621 

Major error. 1 DUAL has found that the problem is infeasible or unbounded; or that the 
mantissa and tolerances are inadequate. The problem status (infeasible or unbounded) and the 
name of the (infeasible or unbounded) variable are given by the preceding error message: 
INFEASIBLE VARIABLE xxxxxx or UNBOUNDED VARIABLE xxxxxx. If the problem is in- 
correctly formulated, use Error Restart A or, if possible, E. If the problem is correctly 
formulated, turn console switches 2 and 4 on, remove cards from read hopper, place a GET- 
OFF agendum card in the read hopper, press START, turn console switches 2 and 4 off, list 
GETOFF output cards, adjust mantissa and/or tolerances by an ASSIGN card preceding the 
MIN or MAX agendum, and use Error Restart H. (See ASSIGN card format shown earlier.) 
The tolerances which may required adjustment are: element tolerance (set this to 9 unless the 
BASIS, card in the GETOFF output is larger than 9); pivot tolerance (increase this by 1 or 2, 
that is, from 4 to 5 or 6); the objective function tolerance — if the solution is unbounded 
only — (decrease this by 1 or 2, that is, from 6 to 5 or 4). Increase the mantissa to 18 (or 
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less). If the element tolerance is changed and the mantissa is not changed, the ERASE I 
(erase the inverse) agendum should precede the MIN or MAX agendum card. 

TYPE F YYYYY NOT IN BETA 1 CHECK. 

NEWRHS 

Major error. YYYYY = SLACK - AJ SLACK file. YYYYY = VBL. - AJ VAR (MATRIX) 
file. A variable is classified as a type F variable in the AJ VAR or AJ SLACK file. This 
variable does not appear in the solution list BETA 1. Error is probably due to file destruc- 
tion or misuse by an optimization or output program: 1DUAL, INVERT, CHECK., OUTPUT, 
DO.D/J, or COST.R. It may be due to user program modification or replacement. 

UNBOUNDED VARIABLE xxxxxx 1 DUAL 

Information message/major error. Variable xxxxxx has been chosen to enter the basis but 
has been found to be unbounded. The system will CHECK, the current basis for accuracy, 
increase the mantissa and tolerances if accuracy (max and min row error above maximum 
error tolerance) is inadequate, and re-INVERT the current basis to minimize possibility of 
digital difficulty. If the variable remains unbounded, the message repeats, this time as a 
major error. See TWO CONSECUTIVE FAILURES IN OPTIMIZING message for correction 
and restart procedure. 

UNKNOWN CODE, A VBL FILE INVERT 

Major error. A basis identification card for a structural variable contains an invalid class. 
Turn console switches 2 and 4 on, remove all cards from read hopper, place GETOFF 
agendum card in the read hopper, press START, turn console switches 2 and 4 off, list 
GETOFF output cards, correct invalid class(es), use Error Restart G. 

UPPER LIMIT DIM = NNNNN LP1620 

Information message. The last Monitor-protected sector on disk drive is NNNNN. Sectors 
04800 to NNNNN are Monitor-reserved for DIM, EQUIVALENCE tables and programs. The 
area from NNNNN to 15999 may be used for LP data if an ASSIGN specifies that the area can 
be used. This ASSIGN card must precede the use (INPUT. , SAVE) of the area above NNNNN. 



DISK STORAGE UTILIZATION 

AREAS UNAVAILABLE ON MONITOR DISK 



The following areas on the Monitor disk are not available for LP files: 



Sectors 


00000- 


-01804 


— Temporary Storage, Common File, LP System 


Tt 


04800- 


-04999 


— DIM Table 


tt 


05000- 


-05199 


— Equivalence Table 


It 


05200- 


-NNNNN ■ 


— Linear Programming and User Programs. The sector address 
NNNNN is given by the Message UPPER LIMIT DIM = NNNNN 


t! 


16000- 


-19999 


— Monitor 



Sectors 00000-04799 are used by Monitor as working storage; if the Monitor disk is being used 
for other jobs as well as LP runs, any files to be saved must be placed above the DIM table, 
Equivalence table, and programs. At the beginning of each LP run, a message is typed out 
giving the sector address of the first available sector above the program area. 
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If the Monitor disk is being used only for LP runs, any files to be saved may be stored in the 
following sectors: 

sectors 01805 through 04799 

first sector beyond programs through 15999 

On multidrive systems, the second and succeeding disks may be used in their entirety for 
saving LP files. 

PROTECTION OF DIM TABLE 

When files are created, they are placed on the disk beginning at the sector address specified 
on the INPUT, agendum card. If they require too much area to fit below the DIM table, 
Equivalence table, and programs, they will be continued beyond the program area, so that the 
DIM table, etc. , will be protected. This is true of all files except the Inverse File, where 
this technique is not possible. If the Inverse File is started below the DIM table and does not 
fit in that area, a message is typed out to that effect and a major error halt occurs. 

With a large problem, the Inverse File will utilize by far the greatest amount of disk storage; 
thus it is to be expected that the Inverse File will not always fit below the DIM table. As 
stated above, the files will be placed on the disk beginning at the sector address specified on 
the INPUT, agendum card; the ASSIGN agendum card may be used to place the Inverse and 
Beta Files on a different disk area. Thus, INPUT, may specify a sector address low enough 
to allow all other files to reside on the disk below the DIM table, while ASSIGN is used to 
place the Inverse and Beta Files beyond the program area. This technique may also be used 
with a multidrive system, where the Inverse and Beta Files are placed on the second disk, 
while the other files are put on the first disk. (An example of this usage is included in the 
paragraph on additional disk drives, later in this section.) 

LP DISK FILES 

The number and size of the generated files depends on the input data. The size of the Inverse 
and Beta Files (see below) depends also on the mantissa length. 

The files generated from the input data or formed in processing are: 

1. Program Control File, containing location, status, and current processing information. 

2. Identification Files 

a. ROW. ID file, formed from card input 

b. COL. ID file, formed from card input or generated from the Matrix file 

c. Basis slack file, formed from card input or SAVEd from Matrix (Aj Variable) file 

d. Basis Variable file, formed from card input or SAVEd from Aj Slack file 

3. Data Files 

a. Aj Variable (Matrix) file, formed from card input 

b. Aj Slack file, generated from ROW. ID file 

c. B (R-H-S) file, generated from card input. 
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4. Solution Level File: Beta file — solution levels and bounds of type F (intermediate-level) 
basis variables. 

5. Transformation Files 

a. Eta 11 Logical basis transformation. A single-digit-per-row coding to indicate the 
type of logical variable formed from the ROW. ID file. 

b. Inverse file — iteration transformation vectors, formed after inversion and by 
iteration process 

c. Alpha file — generated during inversion procedure. Partially transformed basis 
vectors 

LP FILE LOCATION AND SIZE 

The Problem Control File begins at the sector specified by the INPUT, or SAVE.p agendum. 
This file is five sectors in size. The location of all other files is given by the Problem 
Control File. The applicable fields are all in the first two sectors as follows: 

Sector Positions in Length Description 
sector 



01-12 12 Problem name 

15-23 9 Disk drive, sector address, amount of sectors for 

ROW. ID File 

29-37 9 Disk drive, sector address, amount of sectors for 

ETA11 File 

43-51 9 Disk drive, sector address, amount of sectors for 

COL. ID File 

57-65 9 Disk drive, sector address, amount of sectors for 

MATRIX File 

86-90 5 First available sector above DIM table, Equivalence 

table, and programs 

91-95 5 Upper limit of user's area on disk 

01-09 9 Disk drive, sector address, amount of sectors for 

Right-Hand-Side File 

15-23 9 Disk drive, sector address, amount of sectors for Basis 

Variable File 

29-37 9 Disk drive, sector address, amount of sectors for Basis 

Slack File 

43-51 9 Disk drive, sector address, amount of sectors for 

Inverse File 

57-65 9 Disk drive, sector address, amount of sectors for Beta 

File 
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Sector Positions in Length Description 
sector 



71-75 
76-80 

81-85 
86-90 

91-95 
96-100 



Sector address of first Aj group record 

Sector address of first sector beyond last Aj group 
record 



Sector address of first A slack group record 

Sector address of first sector beyond last A group 
record 

Sector address of first Right-Hand-Side group record 

Sector address of first sector beyond last Right-Hand- 
Side group record 



EXAMPLE OF LP FILE ARRANGEMENT 

A typical input deck is shown below, with the files which would be created for the problem. 
Note that the files start at the sector address specified on the INPUT, agendum card. 



Card Deck 
INPUT. C08000 
ROW. ID 
COL. ID 
MATRIX 
FIRST. B 

• 

NEXT.B 

BASIS 

ENDATA 

If no COL. ID data is input, a 
COL. ID File is created 
internally and put: 

1. directly after the RHS File if 
no BASIS cards are input, or 

2. directly after the BASIS File 
if BASIS cards are input. 



Disk Storage 



Sector 08000 
Sector 08005 



COMMON for this problem 



ROW. ID File 



ETA11 File 



COL. ID File 



Aj VARIABLE (MATRIX) File 



Aj SLACK File 



RIGHT-HAND-SIDE File 



BASIS VARIABLE File 



BASIS SLACK File 



BETA File 



INVERSE File (Alphas and Etas) 
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FILE SIZE 



With the exception of COMMON, all files vary in size according to the problem size. To find 
the number of sectors used for a file, take the results of the formula and use the next higher 
integer (for example, 2024/100 = 21 sectors). 



COMMON 



5 sectors 



ROW. ID File 



R x 14 
100 



R = number of rows 



ETA 11 File 



R 

100 



R = number of rows 



COL. ID File 



C x 14 
100 



number of columns 



Aj VARIABLE (MATRLX) FILE AND THE Aj SLACK FILE 

Each of these is composed of group records. A group record contains as many variables as 
possible within the maximum group record size to reduce disk usage. The maximum size of 
a group record depends on machine size. 



max. size of a 
group record: 



20K — 
40K ~ 
60K -- 



20 sectors 
50 sectors 
100 sectors 



MATRIX File 
Aj file 



MATRIX File 
Aj Slack file 



15 + (57C + 15N) 
100 



15 + 72R 
100 



C = number of columns 
(variables) 

N = average number of 
nonzero entries per 
column, including 
bounds 

R = number of rows in 
ROW. ID File 



RIGHT-HAND-SIDE FILE 

The RHS File also consists of group records, whose maximum size is shown below. However, 
when more than one group record is needed, one right-hand side is never split; that is, if a 
right-hand side will not all fit in a group record, it is put into the next group record. 



max. size of a 
group record: 



20K — 31 sectors 
40K — 80 sectors 
60K — 155 sectors 



RHS group 
record 



15 + (15H + 30N) 
100 



H = number of right-hand 
sides in this group 
record 



N = number of nonzero 
row entries in this 
group record 
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BASIS Variables 
File 



V x 14 



100 



V = number of variables 



BASIS Slack 
File 



S x 14 
100 



S 



number of slacks 



BETA File 



R x 44 
100 



R = number of rows 



INVERSE File maximum size 



[ (24 + (M + 2) * R)/100] * R + max (I, G * [F/E] ) 



where: 



E = [(G * 100 - 20)/((M + 5) * R + ll)] d 



[ ]y L is the smallest integer - [ ] 

[ ]^ is the largest integer < [ ] 

R = number of rows 

M = problem mantissa length 

F = inversion frequency 

E = min. number of Eta's per group record 

G = ETA group record size 

I and G vary according to machine size: 



size 


J_ 


G_ 


20K 


32 


24 


40K 


68 


50 


60K 


104 


100 



ADDITIONAL DISK DRIVES 

If more than one disk drive is used, the Monitor would be on the first drive (drive 0) ; the 
second (drive 1) and succeeding drives could then be used in their entirety for LP problems. 
With two disk drives, three different arrangements are possible for an LP problem: 

• Problem files entirely on first disk 

• Problem files entirely on second disk 

• Problem files on both disks 

To have the problem files entirely on one disk, use the appropriate sector address on the 
INPUT, agendum card. Thus, to begin the problem files on the first sector of the second 
disk, specify sector address 20000. 
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It is possible to have the files for one problem on both disks; this is recommended when large 
problems are run. The method for accomplishing this is to use the second disk for the BETA 
and Inverse Files, and the first disk for all others. This is done by using a sector address of 
the first disk on the INPUT, agendum card, and then a sector address of the second disk on an 
ASSIGN agendum card (columns 33-37) . For example: 



Card Deck 

ASSIGN 20000 
INPUT. C 090 20 
ROW. ID 
COL. ID 
MATRIX 
FIRST. B 
NEXT.B 
BASIS 
ENDATA 



Sector 09020 
Sector 09025" 



First Disk 
(Drive 0) 



Second Disk 
(Drive 1) 



COMMON for this problem 



ROW. ID File 



ETA11 File 



COL. ID File 



MATRIX File 



RIGHT-HAND-SIDE File 



BETA File 



INVERSE File 



Sector 20000 



BASIS File 



CORE STORAGE UTILIZATION 

Because each routine in the 1620-1311 LP System is loaded and executed individually and 
works with different files, there is no standard core storage arrangement. Communication 
from routine to routine is accomplished by the use of COMMON, which is stored in sectors 
01800-01804 on the disk and is called in by each routine. 

Monitor makes use of the memory area below 02402; above that area are located program and 
files, and floating-point subroutines in those programs which require them. 

Not only do routines differ from one another in the utilization of memory, but the organization 
in the memory of a specific routine may vary from one machine size to another in order to 
accommodate the files which are increased in size when core storage is increased. 



ERROR AND INTERRUPT RESTART PROCEDURES 



Procedures A through I describe input data and formulation restart procedures. These pro- 
cedures are referenced by error message action. In case of multiple input errors, the 
highest-priority restart procedure (lowest number) should be used. 

Procedure J is the standard GETOFF interruption resumption. 

Procedure K is the standard malfunction restart procedure. 
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ERROR RESTART PROCEDURES 

A. Correction of INPUT data. Priority 1. 

1. Correct card(s). 

2. Put corrected deck in card reader followed by interrupted agenda. 

3. Turn on console switch 4. 

4. Push START. 

5. Turn off console switch 4. 

B. Correction of REVISE data. Priority 4. 

1. Correct card(s). 

2. Put corrected deck in card reader. 

3. Put remaining agenda in card reader behind corrected deck. 

4. Turn on console switch 4. 

5. Push START. 

6. Turn off console switch 4. 

C. Additions to INPUT data, method 1. Priority 2. 

1. Correct original INPUT data and identification cards. 

2. Insert additional entries, identifications. 

3. Put corrected deck in card reader. 

4. Put remaining agenda in card reader. 

5. Turn on console switch 4. 

6. Push START. 

7. Turn off console switch 4. 

D. Additions to INPUT data, method 2. Priority 2. 

1. Insert additional entries, identifications in original input data deck. 

2. Consolidate all revisions: 

a. Order the revise decks from oldest (first) to newest (last). 

b. Remove all REVISE agendum cards except the first one. 

c. Remove all ENDATA and EOF cards except the last one. 

3. Put the INPUT deck, including INPUT, agendum, in the card reader. 

4. Put the consolidated REVISE deck in the card reader. 

5. Put interrupted agenda remaining in card reader. 

6. Turn off console switch 4. 

7. Push START. 

8. Turn on console switch 4. 

E. Corrections to input by a REVISE deck. Priority 3. 

1. Prepare a REVISE deck. 

a. REVISE agendum card. 

b. Data indicator card, ROW. ID, COL. ID, MATRIX or FIRST. B as required. 

c. Corrected identification card or element card. 

d. EOF card. 

2. Place REVISE deck in card reader followed by remaining agenda. 

3. Turn on console switch 4. 

4. Push START. 

5. Turn off console switch 4. 

F. Continuation if error does not result in incorrect disk storage of data. Priority 6. 

1. Turn on console switch 4. 

2. Push START. 

3. Turn off console switch 4. 
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G. Corrections to input basis by a REVISE deck. Priority 3. 

1. Prepare a REVISE deck. 

a. REVISE agendum card. 

b. BASIS, indicator card. 

c. VARBLS indicator card. 

d. Corrected basis identification card if the variable is a structural (column) vari- 
able; if it is a logical (slack) variable, no card. 

e. SLACKS indicator card. 

f . Corrected basis identification card if the variable is a logical (slack) variable; if 
it is a structural (column) variable, no card. 

g. EOF card. 

2. Place REVISE deck in card reader followed by remaining agenda. 

3. Turn on console switch 4. 

4. Push START. 

5. Turn off console switch 4. 



H. Corrections to agenda deck. Priority 1. 

1. Order and correct the agenda cards. 

2. Place remainder of agenda deck in card reader. 

3. Turn console switch 4 on. 

4. Press START. 

5 . Turn console switch 4 off. 



I. Inverse destruction. Priority 5. 

1. Prepare a restart agenda. 

a. Monitor cards; cold start #= 4= JOB, #= * XEQSLP1620. 

b. INPUT. Dsssssnnnnnn, where sssss is the sector address of the input control 
file and nnnnnn is the problem name. 

c. ERASEI to erase the inverse. 

d. Followed by the rest of the agenda. 

2. Place restart agenda in the card reader. 

3. Press RESET. 

4. Press LOAD. 



J. GETOFF restart procedure after standard interruption. When the 1620-1311 LP interrupt 
procedure is used (console switch 2 ON), GETOFF punches out the current basis. To re- 
start from the point of interruption: 

1. Prepare a new problem deck. 

a. Remove the original basis, if any. 

b. Insert the new basis just before the ENDATA or EOF card. 

2. Put the complete deck in the card reader. 

3. Press LOAD. 

K. Standard Machine Malfunction restart. 

1. Prepare a restart deck. 

a. LP1620 RESTORE card. 

Column 1 - 12 34 00032 00701 

13 - 24 36 00032 00702 

25 - 32 49 02402 

33-46 1 00100 100 00100 

b. Current agendum card followed by remaining data and agenda. 

2. Press RESET. 

3. Put the restart deck in the card reader. 

4. Press LOAD. 
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SAMPLE TIMINGS 



Problem ID 


m 


n 


Optimality 


Machine Model and 
Size 


Comments 


Iter 


Time 


SHARE 20 


17 


25 


21 


8 min. 


1620-1 40K 




SHARE 15 


66 


72 


58 


52 min. 


1620-1 40K 




SHARE 19 


78 


87 


131 


1 hr. , 10 min. 


1620-2 60K 




SHARE 1 


20 


31 


15 


3 min. 


1620-2 60K 




SHARE 22 


54 


44 


58 


36 min. 


1620-2 60K 




SHARE 16 


49 


62 


18 


1 hr. , 10 min. 


1620-1 40K 




SHARE 7 


36 


78 




1 hr. , 44 min. 


1620-2 60K 


Several options 
and solutions 


Sample Problem 


9 


7 




38 min. 


1620-1 40K 


For entire 
sequence 



MATHEMATICAL DESCRIPTION OF 1620-1311 LP 

This section describes the formulae and their application to the routines. The notation used 
is the program notation or symbology, wherever possible. For proofs, etc. , the reader 
should consult the mathematical references (5, 10, and 11) in the Bibliography. This section 
assumes a knowledge of vector algebra. 

INPUT. 

The input routine forms the structural and logical matrix files and classifies the variables, 
as specified by the user. The disk output of this routine is: 

aY = b, Y < b, Y > b 

where: 

a is the matrix of coefficients of the input structural variables and the generated logical 
variables. 

Y is the vector of variables. 

b is the right-hand-side vector (in force at optimization) . 
b are the upper bounds of Y. 
b are the lower bounds of Y. 
Note: The subscript-superscript notation will be used throughout this section. 
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For example: 

sli : column j of matrix a. 

b 1 : element (row) i of right-hand-side b. 

: variable j of the (column vector of) variables Y. 

a] : coefficient of i^h row. jth column, 
j 

b** , hf* : elements j of upper, lower bound vectors: b"' < Y"" ^ 

Dimensions . - The number of rows is M, the number of structural variables is N, the total 
number of variables is M + N. 

Tolerance and mantissa assignment . - The input routine assigns initial mantissa and toler- 
ances according to the following formulae, unless specified by a BASIS, input control card or 
an ASSIGN card. All tolerances are expressed as powers of 10: 10 e . Define N.D.R.a! 
as the number of digits to the right of the decimal point of input data element aj . J 



pivot 



■1 - max { 3, N.D.R.a 



3 f , 



1, . . . , M for nonfunctional rows. 



j=l,..., N 



functional 



-1 - max ^3, N.D.R.a! \ 

( J ) j=l N 



= 1, . . . , M for functional rows. 



feasible 



= -1 - max <3, N.D.R.b 



= 1, ... , M 



for j ranging over all right-hand sides. 



error 
Mantissa 



= e 



feasible 



min jl8, max j5, 2 + range - min je pivot , e functi(mal | [ (where range 
is the difference between largest and smallest exponent of the input data 
converted to floating-point form. 

€ element = " truncated j^antissa/4 | + min j€ pivot , € functiQnal } - 2. 
Mantissas and tolerances are modified during input by: 

Mantissa: add quantity (e error - e error ma x) ^ positive. e error max is the exponent 
of the maximum error found in a CHECK. 



element 



= € element " truncated (mantissa 1 st - mantissa old)/2. 



€. ^.e^. L . ,,6^ = truncated (mantissa 1 st - mantissa old)/3. 

pivot functional feasible old 
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Classification of variables , - This is determined according to basis status and activity level. 
W: class of variables at lower bound and not in basis. 



F : class of variables in the basis. 

G : class of variables at upper bound and not in basis. 

*: subclass of nonbasis variables, the omitted rows and columns. 

Reference to variables of the various classes is by superscript, , FJ , GJ ; and the 
column vectors by small letters, wj , fj , gj when differentiation is required. 

Initialization. - The generated logical variables constitute the initial basis, unless a basis is 
specified. The logical variables are classified as type W, unless a basis is specified. The 
omitted rows (class * logical variable in basis) and the omitted columns (class * structural 
variables) are not used in calculation. 

INVERSION 

The inversion procedure is similar to the revised simplex tableau procedure. The pivot 
selection is the major difference. 

Notation . A time index is required to clarify certain vectors, matrices, etc. , within 
parentheses. 

For example: 

F(0) : the basis at iteration 0; the logical basis. 
Right-hand-side initialization . 

/3 : the transformed right-hand side; the current basis activity levels. 

7r: the basis inverse. 

rj (1- 1) : the inverse of the logical basis, f(0). 

Note: In the program listings, 77 (1- 1) is used to reference V (0). For consistency, 
77 (1- 1) is used in this section. 

[f(O)]- 1 =77(1-1) 

Note : 77 (1 - 1) = f(0) . 

The first transformation effects a translation of axes to shift the lower bounds to zero. 
Define: X + b = Y 
Then: a (X + b) = b, 
aX = b - ab, 
Define: a = b - b 
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b < X + b < b 
< X < b - b 



F(0) = p (0) = t? (1- 1) [b - ab - gg] 
F(0) : Initial logical basis. 
/3 (0) : Logical basis activity. 

g : Matrix of vectors of variables at upper bound. 

g : The transformed (a) upper bounds of these variables. 

Logical variable bounds are formed for all range constraints and equations. 

S 1 = b 1 - b 1 For range constraints where: 

S 1 is the upper bound of the logical variable. 

b 1 is the upper bound row limit, b 1 is the lower. 

S 1 = for artificial logical constraints generated for equation rows . 

Alpha matrix . - A matrix of the basis (type F) structural variables is formed. This matrix 
is transformed by pivoting once on each column in the available (vacant logical variable for 
the row) rows until the matrix is completely transformed. Pivot choice is based on minimiz- 
ing the number of arithmetic operations required in the transformation. 

Initialization. - Transform structural variables by logical basis inverse. 

a. = V (l- 1) a. , j = 1, ... , y where aj's are the vectors of type F (Basis) structural 
"V 3 variables. 

Define . - Parameters used to select pivot row. 

r 1 + 1 = Number of nonzero entries in unpivoted a 1 . Row count. 

c. + 1 = Number of nonzero entries in unpivoted a, . Column count. 
J J 

Pivot candidate row: The logical variable for the row is: 

Type W; nonbasic at lower bound. 

Type G; nonbasic at upper bound. 

And row has not been pivoted during inversion. 
Pivot possible position: lo^l > 10 £ P iv °t, oil not pivoted, i pivot candidate row. 
Pivot procedure complete when there are no pivot possible positions. 
Normal completion if all oj- pivoted and no pivot candidate rows remain. 
Pivot Choice . - Minimize transformation operations. 
Choose R = pivot row ; s = pivot alpha 

c g = min j Cj min jr*| | for all pivot positions possible. 
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Pivot procedure . - Pivot selected column and row. Transform remaining alphas and beta. 

Define: I(R, a s ) as an identity matrix (dimension M) except that column R of that matrix 
is a s . 

Define: *7(T) = [I(R, a 



Note : [^(T)]^ ^ = 1^ where I is an identity matrix. 



[1(T)]^ R =-«V< 

Alpha transformation 

[a(T)]. =rj( T ) [a(T-l)]. 

Note: [a(T)]^ R = [a<T-l)]f + t£ [a(T-l)]f 

l i -tt l 

[o^T)]^* =^[«(T-l)] i R 

Beta transformation 

/S(T) = r7(T) ]8(T-1) 

Error conditions. If row candidates and/or alphas remain, but pivoting cannot continue — 
no possible pivot candidates remain — the remaining alphas, if any, are ignored and logical 
variable restored for any remaining pivot candidates. 

The inverse. The product of the etas (V) would form the inverse. The inverse is main- 
tained in the product form. 

tt(T) = *?(T) ?7(T-1) ...tj(I) 77(1-1) 
Inverse row . - During optimization a row of the inverse is required. 

tt 1 = I 1 * = I%(T) rj(T-l) ... t?(1) 77(1-1) 
Notation: Objective function row noted by <p 

Note: Only the vector [t?(T)] and the pivot row indication R are stored. 

R 

DUAL 

Row selection. - Choose R, the basis index or position of the variable to leave the basis such 
that: 

= min (/?*, f* - /?*)J for i = l, m excluding * and objective function artificials. 
K 1 = 1 if fi ^ 0, where f* is the upper bound of F 1 . 
K 1 = 100 if I 1 = 0. 
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The solution is said to be feasible if > - 10 *- eas ^ e m 
The variable F leaves the basis at its lower bound if /3 < 0. 
The variable F R leaves the basis at its upper bound if /3 1 > 0. 
Define: BETA TEST VALUE = (1/K R )/V R . 

Column selection, - Choose s, the index or position of the variable to enter the basis, such 
that: 

Z = min | D./C. I for C. < -10 pivot an( j is not an artificial logical variable 
S (•'"') (generated for an equation row) . 

Where: 

! + 1 if MAXimizing 1 { + 1 for type W variables 1 

- 1 if MINimizing / I - 1 for type G variables / 

This formula is used only for the first iteration following inversion. The objective 
function row is denoted by <p . 

And: 

( - 1 if /3 R > ) { + 1 for type W variables ) 

I + 1 if < J I - 1 for type G variables J 

Dj is computed during dual by Dj = Dj + Z g Cj except for the first iteration following inversion, 

A W variable increases from an explicit lower bound, or from the implicit zero lower bound, 
as it enters the basis. 

A G variable decreases from the explicit upper bound as it enters the basis. 
The problem is said to be infeasible if no Z g exists. 
Processing error test 

e 



1 - 



C /a R 
s s 



error 

< 10 Test failure results in increase of mantissa length, change 

in tolerance and reinversion. 



Bound test 

{ - 1 for type W"| _ 

New BETA TEST VALUE = old BETA TEST VALUE - J > a s • a 

( + 1 for type G ) s 

If the new BETA TEST VALUE is negative, the entering variable replaces the leaving vari- 
able; otherwise the variable enters at its bound. If the B.T.V. remains negative, the vari- 
able would violate its upper or lower bound if it replaced F R . 
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PRIMAL 

Column selection. Choose s, the index or position of the variable to enter the basis, such 
that: 



„ . , „ , , _ ^ - n " € functional 
D g = mm <D. } and D g < - 10 



Where: 



| + 1 If MAXimizing \ ( + 1 for W | 



D j = \ \ # \ ( a j 

U 1 if MINimizing ) I - 1 for G / 

The solution or basis is said to be optimum if no D c exists. 

fc> 

A W variable increases from its lower bound as it enters the basis. 
A G variable decreases from its upper bound as it enters the basis. 

In a degenerate iteration (no change in objective function), a W variable enters at its lower 
bound, and a G variable enters at its upper bound. 

Row selection. - Choose R, the index or position of the variable to leave the basis, such that: 

o 1 / 1 -f 1 ^ 1A f pivot "\ . - 

3/oi if a > 10 I i = 1 m 

s< s I 

yR _ J ( excluding * and 

(_ \ . e . ( objective function 

B - fM /ol if a < -10 and f 1 exists \ artificials. 

/ s s ; 

Indexing stops and row i is selected (i=R) if f * = and I a* I > 10 pivot # 

s 

The solution is said to be unbounded if is undefined and the variable entering the basis is 
type W with no upper bound. 

The variable leaves the basis at its lower bound if a* > 
The variable F^ leaves the basis at its upper bound if al < 
Processing error test 

e 



1 - 



V^s 



error 

< 10 . Test failure results in increase of mantissa length, change 

in tolerance and reinversion. 



Bound test 



If a < V the variable enters at bound. 

R — s 

If V is undefined and a exists, the variable enters at bound. 
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UPDATE 

a Q UPDATE TO CURRENT BASIS 
a s = 7ra s 

= [r?(T-l) [*7(T-2)... [1(1) [7,(1-1) [a g ] ]]...]] 
Define: [ce(0)] g = [1(1-1)] [a g ] 

Where 1(1-1) is the inverse of the logical (starting) basis f(0) . 

Note: [f(O)]" 1 = 77(1-1) = [f(0)] 
Define: [&(i)] = [^(i)] [^(i-l)] » anc * is the pivot row for iteration i. 

o S 

Then: [«(i)]^ R(i) = [a(l-l)]* + [«(i-l>]* (1) 
r ...tKH) k r ,. ,,,R(i) 

[o(1 >V "W^v 



and 

a(T-l) = a s 

1(T) Formation . - Required for vector types W and G entering the basis at intermediate 
level (becoming basic or type F), where T is the iteration index. 

7](T) = [HR.ag)]- 1 

th 

where I (R, o> s ) is an identity matrix (dimension M) except for the R— * column which is a £ 

Then: [,(T)]^ R = ^ 

and 

emir 



and 



J R 

Only the index, R, and the vector [tj(T)]-^ are stored. 

Beta update to current basis 

/3(T) = jS(T-l) - <j> + 6 + y 
Where 

/3(T-1) is the vector of basis activity levels at iteration T-l . 
j3(T) is the vector of basis activity levels at iteration T . 
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A zero column vector except for <p^, i.e. , = 

/ T? s 
= < v xv = OifX enters at upper or lower bound. 

b = /3(T-1) if the entering variable enters at intermediate level. 



-a s a„ if the variable is W entering at upper bound. 

o s 

+ a s o> 3 if the variable is G g entering at lower bound. 
R 

[jS(T-l)] V (T) if the variable is W s or G g entering at intermediate level and 
Fj^ is leaving at lower bound. 

[ j/3(T-l)| ^-^f 11 ] V(T) if the variable is W g or G g entering at intermediate level 

and F^ is leaving at upper bound. 

A zero column vector except for Y^, i.e. , y* 5 *^ = 
y <^ y& = a s jf f- ne entering variable is G g entering at intermediate level. 
y = otherwise. 



NEW RIGHT-HAND SIDE 

Initialization. Same as inversion except for alpha formation. 
Beta transformation by inverse 

PIT) = 7T/?(0) 

= r?(T) t?(T-1) . . . 0(0) 

Same as a s update to current basis except that the vector is ]S(0) and is not transformed 
by 77(1-1). 

CHECK. 

Row level. The product-sum of the row coefficients multiplied by lower bounds (all vari- 
ables), basis activity level (type F structural variables), and upper bounds (type G structural 
variables) . 

L 1 = a 1 b + fp* = g X g* 

* structural variables only. 

Selective output includes product-sums only on the selected structural variables (for each 
selected row) . 

Row error . The difference between input value, adjusted by slack level and row level. 

E i _ B i _ for row type and + Is*- L 1 

\ - 1 : for row type blank and - J 
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For selective output this lumps the row error with the product-sum of the nonselected 
variables . 

B 1 is the right-hand-side entry. 

S 1 is the activity level of the logical variable for row i . 
Row limit values . - The current right-hand-side coefficients for the row. 
COST. RANGE 

High-limit variable selection . Choose H, the index or position of the variable which limits 
the highest cost of basis variable F^, such that: 

Z T = max |d. /C. X for C. > 10 pivot and X J is not an artificial logical variable. 
H I J J ) J 

Where: 

_ ( +1 if MAXimizing\ $ ( + 1 for type W variables \ 
j \-l if MINimizing J 71 (-1 for type G variables j j 

The Dj's are not computed during COST. Ranging, but are assumed to be left over from the 
final primal iteration. 

And: 

_ R ( + 1 for type W variables \ 
j I - 1 for type G variables f j 

If Zj[ is undefined, there is no high-limit variable. 

Lo-limit variable selection. Choose L, the index or position of the variable which limits 
the lowest cost of basis variable F^, such that: 

( 1 6 * t 

Z T , min <D./C. > for C. < -10 and X 1 is not an artificial logical variable. 

See high-limit variable definition for Cj and Dj . If Z is undefined, there is no lo-limit 
variable. 

Cost-range calculation 

, . , , . f-1 if MAXimizing) r (-1 if MAXimizing \ 01 ~ 
The highest cost = < , „^ T . . . > Z R + < , ... . . } a. if Zu exists 

& \+l if MINimizing j L H \+ 1 if MINimizing j k J n 

i j. A C-l if MAXimizing) r (-1 if MAXimizing) 1 „ 

The lowest cost = < - . . > Z T + < , . . 6 > a, if Z T exists 

(+1 if MINimizing j L L 1+1 if MINimizing j k J ^ 

Where aj^ is the input objective row coefficient. 

For undefined Zjj , Z ^ "infinity" is placed in the data field. 

A minus sign follows this word when the sense is negative. 
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Current cost. Input (or revised) objective row coefficient. 

Reduced cost . Product of inverse objective function row by variable, reversed in sign for 
variables at upper bound, type G. 



Basis value . Difference of current cost and reduced cost. 

v _ a <P / +1 if MINimizing \ 
j j " \ -1 if MAXimizing ) j 

where a^ is the objective row coefficient. 

Note: Basis value is output only if there is a current cost entry. 

Increase B value, decrease B value: Simplex multipliers or marginal values. Entries only 
for nonzero simplex multipliers 
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